Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,478,3110AG41.0% 43.8 / 14.9 41E29G (GAG→GGG) ynbBputative CDP‑diglyceride synthase
Reads supporting (aligned to +/- strand):  ref base A (6/17);  new base G (16/0);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 2.06e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAAT  >  NC_000913/1478166‑1478448
                                                                                                                                                 |                                                                                                                                         
gCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGgagaga                                                                                                                                       <  1:414288/150‑1 (MQ=255)
     tctGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGGGAAAAggg                                                                                                                                 >  2:558380/1‑151 (MQ=255)
        ggACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAggggg                                                                                                                               <  2:44213/150‑1 (MQ=255)
            tGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGGGGAAAGGGGGGGAgg                                                                                                                          >  2:545020/1‑151 (MQ=255)
               gATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTa                                                                                                                       <  1:534130/151‑1 (MQ=255)
                       gTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGt                                                                                                               <  1:710520/151‑1 (MQ=255)
                                 tgtgCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATCCGAGAGAAAAGGGGGGGGGGGTTAATTACgcctttgttcc                                                                                                      >  2:359683/1‑140 (MQ=255)
                                  gtgCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTTCCt                                                                                                     <  1:576793/150‑1 (MQ=255)
                                  gtgCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCtg                                                                                                    <  2:705744/151‑1 (MQ=255)
                                  gtgCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCtg                                                                                                    <  2:383451/151‑1 (MQ=255)
                                      ttcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGAGGGAAAAGAGGGGGGGGGTTCATTAAGTATTTGTAGCttgtg                                                                                                 >  1:625297/1‑150 (MQ=255)
                                        cttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTTCGGGGGGAAAAGGGGGGGGGGGTAAAATAAGTAATTGGACCCGgggggg                                                                                              >  2:289729/1‑151 (MQ=255)
                                         ttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGGGAAAAGGGGGGGGGGGTACATTTCGTATTTGTACCTGgggggg                                                                                              >  2:787009/1‑150 (MQ=255)
                                                 gcTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGGGGGGAAAGGGGGGGGGGGTTCATTACGGATTTCTACCTGGTgggg                                                                                              >  2:103392/1‑142 (MQ=255)
                                                   tGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGGAAAGGGGGGGGGGGTACATTTAGTATTTGTACCTGGGGGGGGATCGTCAttt                                                                                   >  2:217966/1‑151 (MQ=255)
                                                    gAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCAttt                                                                                   <  1:391562/150‑1 (MQ=255)
                                                    gAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGACGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTg                                                                                  >  2:698483/1‑151 (MQ=255)
                                                           aaTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCa                                                                           <  2:420550/151‑1 (MQ=255)
                                                               cTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGATATTGGGGGGAAAAGAGGGGGGGGTTACATTACATATTCTTTGCTCGTGGGGGAGCGTCTGTTGttttt                                                                             >  1:178198/1‑145 (MQ=255)
                                                                   aaaaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGGGAGAAAAAGGGGGGAGGGTACATTACGTATTCGTACCTGGGGGGGGATCGTCCTTTGTTTTTTAATgggg                                                                     >  2:729047/1‑147 (MQ=255)
                                                                    aaaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTACGGGGAGAAAAAGGGGGGAGGGTACATTACGTATTCGTACCTGGTGGGGGGTCGGCATTTTTTTTTTTATGGGGAtt                                                                  >  2:89042/1‑151 (MQ=255)
                                                                       aTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGCGAGAAAAGGGGGGGGGGGTACATTTCGTAATTTTTCCTGGGGGGGGAAAGTTATTTTTTTTTTATTTGGGAttt                                                                 >  1:775019/1‑149 (MQ=255)
                                                                         ctctGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTGCCGGAGAGAAAAGGGGTTGAGGTTACATTACGTATTCGTCCCTGGTGGGGCATCGTCATTTGTTTTTCATTGGTGAt                                                                   <  1:101728/145‑1 (MQ=255)
                                                                           ctGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAGGGGGTGAGGGTACATTTCGTATTTGTACCTGGGGGGGCGTCGGCATTTTTTTTTTAATTGTGATTTTAgg                                                             >  1:674932/1‑149 (MQ=255)
                                                                           ctGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAg                                                              <  1:275465/148‑1 (MQ=255)
                                                                                 acacTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCt                                                     <  2:141797/151‑1 (MQ=255)
                                                                                  cacTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTg                                                    <  2:324516/151‑1 (MQ=255)
                                                                                   acTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCCGAGAGAAAAGGGGGGGGGGGTACATTATGTATTCGTACCTGGGGGGGGGTCGTCATTTTTTTTTTTATGGTGATTTTAGGTCCACGCTgg                                                   >  2:773030/1‑151 (MQ=255)
                                                                                              tgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgttga                                         <  2:374720/150‑1 (MQ=255)
                                                                                                  aaTTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGGGAGAAAAGGGGGGGGGGGTACATTACGTATTTGTTCCTGGTGGGGGATCGTCATTTTTTTTTTTATTGTGATTTTAGGGCCCCGCTGGGTGACGTTGACGttt                                    >  2:746690/1‑151 (MQ=255)
                                                                                                    ttttAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGGGGAAAAGGGGGGGGGGTATATTAAGTATTTGTACCTGGGGGGGGATCGTCATTTTTTTTTTTATGGTGATTTTAGGTCCCCGCTGGGTGACGTTGACGttttt                                  >  2:37966/1‑151 (MQ=255)
                                                                                                    ttttAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGGGAAAAAGGGGGGAGGGTATATTACGTTTTTGTACCTGGGGGGGGATCGTCATTTTTTTTTTTATTGTGATTTTTGGGCCACGCTGGATGACGGTGGCGgtttt                                  >  1:466331/1‑151 (MQ=255)
                                                                                                       tAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGGGAGGAAAGGGGGGGGGGGTTCAATAAGTATTTGTACCCGGGGGGGGATCGTCATTTTTTTTTTTATTGTGATTTTAAGGCCACGCTGGGTGACGTTGGCGtttttt                                 >  2:230855/1‑149 (MQ=255)
                                                                                                       tAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAGGGGGTGAGGGTACATTAAGTATTCGTACCTGGGGGGGGATCGTCATTTTTTTTTTTATGGTGATTTCAGGTCCACGCTGGGTGACGTTGACGtttttt                                 >  1:250215/1‑149 (MQ=255)
                                                                                                                tGATAAACCGTTTTCTGCTATGGCGGTTTCCCGCGAGAAAAAGGGGGGAGGGTTCATTTCGTATTTTTACCTGGTGGGGGGTCGTCCTTTGTTTTTTTATGGTGGTTTTAGGTCCCCGCTGGGTGAAGGTGTCGttttttttgttggt                         >  1:59453/1‑140 (MQ=255)
                                                                                                                  aTAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAAGGGGGGGGGGTACATTACGTATTTGTAACTGGGGGGGGATCGTCCTTTGTTTTTTAATGGGGATTTTAGGGCCACGCTGGGTGACGGTGACGttttttttggtggt                         >  2:114048/1‑138 (MQ=255)
                                                                                                                    aaaCCGTTTTCTGCTATGGCGGTTAACGGGGGGAAAAGGGGGGGGGGGTACATTAAGTATTCGTACCTGGTGGGGGATCGTCATTTGTTTTTTAATGGGGATTTTAGGGCCACGCTGGATGACGGTTGCGTTTTTTTCGCTGATAAGtttt                  >  2:460944/1‑151 (MQ=255)
                                                                                                                      aCCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCt                <  1:24053/151‑1 (MQ=255)
                                                                                                                                   aTGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAg     <  2:580496/149‑1 (MQ=255)
                                                                                                                                   aTGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGa    <  1:302220/150‑1 (MQ=255)
                                                                                                                                     ggCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAat  >  2:306625/1‑150 (MQ=255)
                                                                                                                                       cGGTTACCGGAGAGAAAAGGGGGTGAGGTTCCATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTGGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAat  <  1:570290/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                         
GCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAAT  >  NC_000913/1478166‑1478448

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