Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,093,686 T→C 100% V130A (GTA→GCA)  ycdT → predicted diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,093,6860TC98.1% 175.8 / ‑4.5 53V130A (GTA→GCA) ycdTpredicted diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (25/27);  minor base G (1/0);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTG  >  AP009048/1093619‑1093750
                                                                   |                                                                
gtCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAt                                                                <  1:116524/70‑1 (MQ=255)
 tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTg                                                              <  1:4407875/71‑1 (MQ=255)
   cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                            >  1:4339955/1‑71 (MQ=255)
       aTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAAt                                                        <  1:1523834/71‑1 (MQ=255)
        tGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAAt                                                        >  1:2549892/1‑70 (MQ=255)
            aTTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                   <  1:4601475/71‑1 (MQ=255)
             ttGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                            <  1:3689056/61‑1 (MQ=255)
              tGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                   >  1:3638681/1‑69 (MQ=255)
                cAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATa                                                       <  1:314457/63‑1 (MQ=255)
                   tCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAt                                                                <  1:513793/51‑1 (MQ=255)
                   tCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGa                                             >  1:3764622/1‑70 (MQ=255)
                     tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                            >  1:2339493/1‑53 (MQ=255)
                     tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                          >  1:1306780/1‑71 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaa                                           <  1:1696195/68‑1 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                            <  1:847321/67‑1 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAac                                         <  1:2541861/70‑1 (MQ=255)
                         attTGTTTCGTCAGATGAGTTTGTGCATATTAATTTTTCTGGCATTg                                                              <  1:3572279/47‑1 (MQ=39)
                         attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                   >  1:3852549/1‑58 (MQ=255)
                         attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                       <  1:3003305/70‑1 (MQ=255)
                          tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAg                                                >  1:4188294/1‑60 (MQ=255)
                          tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGTGAATAAAGTTAGTGAAAACACaa                                     <  1:3525125/70‑1 (MQ=255)
                            tGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                   >  1:2135016/1‑71 (MQ=255)
                             gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                   <  1:4020141/70‑1 (MQ=255)
                             gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                   <  1:3477996/70‑1 (MQ=255)
                              tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                                 <  1:957090/71‑1 (MQ=255)
                                  gTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGt                                               >  1:1375121/1‑53 (MQ=255)
                                   tCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                          <  1:3015482/57‑1 (MQ=255)
                                     aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                          >  1:1324241/1‑55 (MQ=255)
                                     aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAat                          >  1:3040917/1‑71 (MQ=255)
                                     aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAat                          >  1:1530458/1‑71 (MQ=255)
                                      gATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaaa                                                     <  1:2222134/43‑1 (MQ=255)
                                      gATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatt                         >  1:725424/1‑71 (MQ=255)
                                       aTGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGt                                               >  1:352643/1‑48 (MQ=255)
                                       aTGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGt                                               >  1:561966/1‑48 (MQ=255)
                                         gAGTTTGTGCATATTAATATTTCTGGCATTGGTGaa                                                         >  1:11285/1‑36 (MQ=255)
                                          aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttat                      <  1:4043144/70‑1 (MQ=255)
                                            tttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattt                        <  1:515050/66‑1 (MQ=255)
                                              tgtgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttat                      >  1:1309709/1‑66 (MQ=255)
                                              tgtgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTc                  <  1:3771172/70‑1 (MQ=255)
                                               gtgCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                   <  1:1770770/36‑1 (MQ=255)
                                                  cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaa              >  1:4067273/1‑70 (MQ=255)
                                                    tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaa           <  1:1743274/71‑1 (MQ=255)
                                                     atTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaaa          >  1:2094690/1‑71 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                                 >  1:620147/1‑47 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                                 >  1:614728/1‑47 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTcc                 <  1:3470366/63‑1 (MQ=255)
                                                       tAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaa            >  1:37782/1‑67 (MQ=255)
                                                        aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATg        <  1:3039354/70‑1 (MQ=255)
                                                          tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                   <  1:2793995/41‑1 (MQ=255)
                                                          tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACt     >  1:971042/1‑71 (MQ=255)
                                                            tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt   <  1:1487509/71‑1 (MQ=255)
                                                             ttCTGGGATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtg  >  1:4445436/1‑71 (MQ=255)
                                                             ttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtg  >  1:1711690/1‑71 (MQ=255)
                                                                   |                                                                
GTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTG  >  AP009048/1093619‑1093750

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: