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

GCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTT  >  NC_000913/1478166‑1478459
                                                                                                                                                 |                                                                                                                                                     
gCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGgaga‑g                                                                                                                                                   <  2:601752/149‑1 (MQ=255)
gCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGA‑Gaa                                                                                                                                                 <  1:280812/151‑1 (MQ=255)
 cTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGggga‑gga                                                                                                                                                 >  1:20428/1‑148 (MQ=255)
      ctGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGA‑GAAAAgg                                                                                                                                             <  2:581987/149‑1 (MQ=255)
           cTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGG‑GAAAAggggggggg                                                                                                                                      >  2:348978/1‑147 (MQ=255)
           cTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCCGGGA‑GAAAAagggggg                                                                                                                                        >  1:302510/1‑147 (MQ=255)
            tGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTTCGGAGA‑GGAAAAGGGGGGAgg                                                                                                                                     >  2:472666/1‑151 (MQ=255)
             gCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGA‑GAAAAGGGGGTGAGGt                                                                                                                                    <  1:474434/151‑1 (MQ=255)
               gATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGA‑GAAAAGGGGGTGAGGTTa                                                                                                                                  <  2:489605/151‑1 (MQ=255)
                aTTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGA‑GAAAAgggggggggggt                                                                                                                                   >  2:372064/1‑142 (MQ=255)
                    aTCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATCGGAGG‑GAAAAGGGGGGGGGGGTACa                                                                                                                                >  1:446042/1‑148 (MQ=255)
                             ggAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGA‑GAAAAAGGGGGGGGGGTACATTTCGTAtt                                                                                                                       >  1:314941/1‑148 (MQ=255)
                              gAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGGGA‑GGAAAGGGGGGGGGGTGACATTTAGTAtt                                                                                                                       >  2:271437/1‑147 (MQ=255)
                                                cgcTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCa                                                                                                 <  1:664147/151‑1 (MQ=255)
                                                     aaaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCCGAGA‑GGAAAGTGGGGGAGGTTACATTACGTATTTGTAACTGGGGGGGGATCGTCAtt                                                                                               >  2:129313/1‑148 (MQ=255)
                                                      aaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGGGG‑GAAAAAGGGGGGAGGGTACATTTCGTATTCGTACCTGGTGGGGGCTCGTCATTTGtt                                                                                           >  1:387377/1‑151 (MQ=255)
                                                        tataATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGA‑GGAAAGGGGGGGGGGGTAAATTTAGTAATCGTAACTTGTGGGGGATAGTCCtttttttt                                                                                         >  1:403295/1‑151 (MQ=255)
                                                        tataATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGGGG‑GAAAAAGGGGGGAGGGTTGATTATCTCTTTCGAACTGGTGGGGGATCGGCATTTGgttt                                                                                         >  1:228791/1‑151 (MQ=255)
                                                           aaTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTc                                                                                       <  2:20427/150‑1 (MQ=255)
                                                                   aaaaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGAGG‑GAAAAAGGGGGGGGGTTTCATTATGTATTTGTACCTGGGGGGGGATCGTGATTTGtttttttat                                                                                    >  1:442868/1‑141 (MQ=255)
                                                                            tGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTACGGGGA‑GAAAAAGGGGGGAGGTTAAATTTCGTATTTGTACCTGGGGGGGGATCGTCCTTTTTTTTTTTATGGGGGTTTTAGGGcc                                                                     >  1:145951/1‑151 (MQ=255)
                                                                             ggCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGGGA‑GAAAAAGGGGTGGGGGTACATTACGGATTTGGACCTGGTGGGGGCTCGTCATTTTTTTTTCTATGGGGGTTTTAGGTccc                                                                    >  2:582714/1‑150 (MQ=255)
                                                                                      tttCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAt                                                            <  2:256493/150‑1 (MQ=255)
                                                                                        tCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtga                                                          <  1:489873/150‑1 (MQ=255)
                                                                                         cgcgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGA‑GAAAAGGGGGGGGGGATACATTTACTATTTGTACCTTGGGGGGGATTGGCAATTTTTTTTTTATTGTGGttttcaggccaccccgg                                                              >  1:290805/1‑131 (MQ=255)
                                                                                          gcgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGGGG‑GGAAACGGGGGGGGGGTAAAATACTTATTTGTAACTCGGGGGGGGACGTCATTTTTTTTTTTATGGTGATTTTATGTCCCCCCGGGAt                                                            >  1:351056/1‑146 (MQ=255)
                                                                                           cgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACCGGGA‑GAAAAAGGGGGGGGGGTACAATACGTATTTTTACCTCGTGGGGGATCGTCCTTTTTTTTTTTATGGTGATTTTAGGTCCACCCTTGAtggcggt                                                      >  1:156823/1‑149 (MQ=255)
                                                                                            gttgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACCGTGG‑GAAAAGGGGGGGAGGGTATATTAAGTATTTGTACCTTGTGGGGGGTCGTCCTTTTTTTTTTTATTGTGGTTTTTGGGCCACGCTGGGtgacggt                                                      >  1:93195/1‑148 (MQ=255)
                                                                                              tgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgttga                                                    <  1:514569/150‑1 (MQ=255)
                                                                                                    ttttAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGA‑GAAAAGGGGGGGAGGGTACATTACGTATTCGTACCTGGTGGGGGATCGTCATTTGTTTTTTAATGGTGATTTTAGGTCCACGCTGGGTGACGTTGTCGttttt                                             >  2:152380/1‑151 (MQ=255)
                                                                                                        aGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGA‑GAAAAGGGCGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTgcg                                         <  1:116380/151‑1 (MQ=255)
                                                                                                        aGCAACGCTGATAAACCGTTTTCTGCTATGGCGGGTTCCGGGGAGGAAAGGGGGGGGAGGAT‑CATTACGTTATCGTACCCGGGGGGGGAACGGCATTTTTTTTTTAATGGTGAGTTTTGGACCCCGCAGGGTGGCGTTTACGgttttt                                            >  2:593020/1‑148 (MQ=255)
                                                                                                           aaCGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCCGGGA‑GGAAAGGGGGGGGGGGTACATTACGTATTTGTAACTGGTGGGGGGACGTCATTTTTTTTTTTATGGGGGTTTTAGGGCCACGCTGGGTGCAGTTGACGtttttttt                                          >  1:523326/1‑145 (MQ=255)
                                                                                                             cGCTGATAAACCGTTTTCTGCTATGGCGGTTATTGGGGA‑GAAAAAGGGGGGAGGTTTTATTACGCATTTGGGACTCGGGGGGGAACGTCATTTTTTTTTCCATTGTGATTTCATGGCCACGCCTGGGGACGTTGTCGGTTTTTtcg                                         >  1:507225/1‑146 (MQ=255)
                                                                                                             cGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGAt                                    <  2:397147/151‑1 (MQ=255)
                                                                                                              gCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGGGG‑GAAAAAGGGGGGGGGGTACATTAAGTATTTGTTACTGGTGGGGGATCGGCATTTGTTTTTTAATGGGGGTTTTAGGGCCACGCTGGGTGACGTTGACGgttttt                                            >  2:68235/1‑142 (MQ=255)
                                                                                                                 gATAAACCGTTTTCTGCTATGGCGGTTACCGGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATaa                                  <  1:348975/149‑1 (MQ=255)
                                                                                                                    aaaCCGTTTTCTGCTATGGCGGTTACCGGGGA‑GAAAAAGGGGGGGGGTTACATTTCGTATTCGTACCCGGGGGGGCATCGTCATTTTTTTTTTAATGGTGGTTTTAGGTCCACGCTGGGTGACGTTGACGTTTTTTGCGCTGATAAGtttt                             >  2:359643/1‑151 (MQ=255)
                                                                                                                    aaaCCGTTTTCTGCTATGGCGGTTACCGGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGtttt                             <  2:292400/151‑1 (MQ=255)
                                                                                                                       ccGTTTTCTGCTATGGCGGTTACCGGGGA‑GAAAAGGGGGGGGGGGTACATTTCGTATTTGTACCTGGGGGGGGATCGTCATTTTTTTTTTAATGGGGATTTTAGGGCCACGCTGGGTGACGTTGGCGTTTTTTTTGCTGGTAAGtttttt                           >  1:238935/1‑148 (MQ=255)
                                                                                                                                gCTATGGCGGTTACCGGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTg                    <  2:257584/148‑1 (MQ=255)
                                                                                                                                   aTGGCGGTTACCGGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAg                <  2:479249/149‑1 (MQ=255)
                                                                                                                                      gCGGTTACCGGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAatat           <  1:83457/151‑1 (MQ=255)
                                                                                                                                           tACCGGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTAc      >  1:173371/1‑151 (MQ=255)
                                                                                                                                             ccGGTGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGAGATTTCAGGTCCATGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGc    >  1:384769/1‑151 (MQ=255)
                                                                                                                                               ggAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCtt  >  2:665826/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCTTCTCTGGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTT  >  NC_000913/1478166‑1478459

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