Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6620TA42.0% 71.9 / 22.3 50E384D (GAA→GAT) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (8/21);  new base A (21/0);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 6.50e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATC  >  NC_000913/3508526‑3508809
                                                                                                                                        |                                                                                                                                                   
aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGc                                                                                                                                       <  1:769461/151‑1 (MQ=255)
 aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCAAACACCGTAGcc                                                                                                                                      >  1:544768/1‑151 (MQ=255)
    cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAATTAAACAAGGTAGccccc                                                                                                                                   >  2:122709/1‑151 (MQ=255)
    cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAATTAAAAACGGGAGccccc                                                                                                                                   >  2:571053/1‑151 (MQ=255)
    cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGgcccc                                                                                                                                   >  2:58297/1‑151 (MQ=255)
       agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCcaca                                                                                                                                <  2:217748/151‑1 (MQ=255)
       agagaTCAATATGCAGATGCTCATGGGCGAGGGCGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAATTCACACACGGGAGCCCCcac                                                                                                                                 >  1:217735/1‑150 (MQ=255)
             cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAATCACACACGGGAGCCCCCacaacaaa                                                                                                                           >  2:674793/1‑148 (MQ=255)
                   gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtg                                                                                                                    <  2:448212/151‑1 (MQ=255)
                     aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtgg                                                                                                                   <  1:487320/150‑1 (MQ=255)
                          cTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGTAGCCCCCACCTCAGAAAACgtgggg                                                                                                                 >  1:531497/1‑145 (MQ=255)
                             aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCacac                                                                                                          <  2:104837/151‑1 (MQ=255)
                             aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGAGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTAGAGCGGCGAAAAATTCAAAAACCGCAACCCCCCCAACAAAAACCgaggggggc                                                                                                              >  1:298954/1‑142 (MQ=255)
                                 gCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAACACGGGAGGCCCCCCCTCAGAAAACgcgggggggcccccc                                                                                                        >  2:39564/1‑138 (MQ=255)
                                  cGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTAAAAACCGGAGCCCCCCCATCAGAAAGAGGGGGGCGGCCAcc                                                                                                         >  2:408484/1‑147 (MQ=255)
                                    aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTg                                                                                                   <  1:185815/151‑1 (MQ=255)
                                     gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGCGCCGGGGAAAAATCACAAACGGAAGCCCCCACATCACAAAACGCGGGGGGCCCAcc                                                                                                         >  1:61938/1‑144 (MQ=255)
                                        tgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAGCGTGGGGCGCACACCGAgcggcg                                                                                                 >  1:303220/1‑149 (MQ=255)
                                        tgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATCAAAAACCGGAGACCCCCCAACCCAAAGAGGGGGGCGcaca                                                                                                           >  2:296260/1‑139 (MQ=255)
                                         gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa                                                                                               <  2:762432/150‑1 (MQ=255)
                                           aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCg                                                                                            <  2:311300/151‑1 (MQ=255)
                                             cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGc                                                                                             <  1:545189/148‑1 (MQ=255)
                                                  cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACAACGGAGGCCCCACATCAGAAAAAGTGGGGCGCACACCCAGCTGCGTAGCGATCATcc                                                                                     >  1:374033/1‑151 (MQ=255)
                                                    gATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGTGGGCGCCGAAAAAATCAAAAACAGTAACCCCCCCATCAGAAAGCCTGGGGCGGACCCCCGGCCACGGAGCGAACAccccc                                                                                   >  1:630454/1‑150 (MQ=255)
                                                     aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                                                                  <  2:81994/151‑1 (MQ=255)
                                                            cTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAGCCCCCACATCAGAAAACGTGGTGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGAt                                                                            >  1:627669/1‑150 (MQ=255)
                                                              cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaaa                                                                         <  2:293395/151‑1 (MQ=255)
                                                               aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCAAACACGGTAGCCCCCACATCAGAAAACGTGGTGGGCACAACCAGCAGCGTAGCGATCATCCCTTGCTGATaaa                                                                         >  2:492362/1‑150 (MQ=255)
                                                                 aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAATTACACACGGTAGCCCCCCCATCAGAAAGCGGGGGGGGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAcc                                                                       >  1:240761/1‑150 (MQ=255)
                                                                 aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                                                                      <  1:35678/151‑1 (MQ=255)
                                                                       tctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACCGTAGCCCCCCCCTCAAAAAGCGGGGGGCGCACACCCAGCTGCGGAGCGATCATCCCCTTCTGGTAAAACAGCa                                                                   >  1:733052/1‑148 (MQ=255)
                                                                          ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGAAAAAATTACACACAGTAGCCCCCACCTCACAAAGCGGGGGGCGCCCACCCAGCTGCGTAGCGAGCATCCCTTGCTGGTAAAACAACACTgg                                                               >  2:728190/1‑149 (MQ=255)
                                                                          ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCACATCAGAAAACGTGGGGGGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCAc                                                                  >  2:369223/1‑146 (MQ=255)
                                                                          ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTgg                                                               <  1:78024/149‑1 (MQ=255)
                                                                                aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACAAACGGGAGCCCCCCCATCAGAAAACGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGGTAAAACAACACTGGcacc                                                           >  1:130886/1‑147 (MQ=255)
                                                                                 aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                      <  1:645063/151‑1 (MQ=255)
                                                                                     gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg                                                  <  2:260204/151‑1 (MQ=255)
                                                                                       gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGt                                                <  2:720782/151‑1 (MQ=255)
                                                                                                gaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGTAGCCCCCACATCAGAAAGCGGGGGGGGCACACCCAGCCGCGTAGCGGTCATCCCTTGCTGGTAAACCAACACTGGCACCACCTCGCGGGGAtggtgg                                         >  1:237507/1‑149 (MQ=255)
                                                                                                       gCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCa                                <  2:672105/151‑1 (MQ=255)
                                                                                                        cGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACCGGAGCCCCCACATCAGAAAGCGGGGGGGGGACACCGCGCTGCGTCGCGAACATCCCTTGCTGATAAAACAGCAATGGCACCCCTTCGCGGGGATGGTGGGTGTGAcc                                 >  2:410994/1‑149 (MQ=255)
                                                                                                                ttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCgg                        <  1:145909/150‑1 (MQ=255)
                                                                                                                  gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGAt                      <  1:523710/150‑1 (MQ=255)
                                                                                                                             gcgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAt          <  2:698693/151‑1 (MQ=255)
                                                                                                                              cgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc         >  2:715475/1‑151 (MQ=255)
                                                                                                                              cgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc         >  2:149361/1‑151 (MQ=255)
                                                                                                                                cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg       <  1:107008/151‑1 (MQ=255)
                                                                                                                                cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg       >  2:392224/1‑151 (MQ=255)
                                                                                                                                  gAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc     >  2:435023/1‑151 (MQ=255)
                                                                                                                                      aaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATc  >  1:360235/1‑150 (MQ=255)
                                                                                                                                        |                                                                                                                                                   
AAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATC  >  NC_000913/3508526‑3508809

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