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

GTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508516‑3508811
                                                                                                                                                  |                                                                                                                                                      
gttgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATca                                                                                                                                                     >  2:156796/1‑150 (MQ=255)
     ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCacacac                                                                                                                                                <  2:189231/150‑1 (MQ=255)
       ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACgg                                                                                                                                              <  1:1472/150‑1 (MQ=255)
         aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGCGCGGGGGGAAAATCacacacag                                                                                                                                              >  1:29278/1‑146 (MQ=255)
         aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAATTCacacaccg                                                                                                                                              >  2:242817/1‑146 (MQ=255)
              cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAACAAGGGAG‑ccccc                                                                                                                                     >  1:566330/1‑151 (MQ=255)
                   agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAG‑CCCCCACATc                                                                                                                                >  2:650042/1‑151 (MQ=255)
                     aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATAAAACACGGTAA‑ccccc                                                                                                                                     >  2:824461/1‑144 (MQ=255)
                       cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAG‑CCCCcac                                                                                                                                   >  2:418771/1‑144 (MQ=255)
                        aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCCCGGGAAAAATCAAAAAAGGGAG‑ACCCCACACCAGaa                                                                                                                            >  1:399532/1‑150 (MQ=255)
                           atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAG‑CCCCCCCCTCAGaaa                                                                                                                           >  2:876858/1‑148 (MQ=255)
                             gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCgtg                                                                                                                      <  1:198283/151‑1 (MQ=255)
                               aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCgtggt                                                                                                                    <  1:467949/151‑1 (MQ=255)
                                     tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGTGCGGGGGAAAAATAACACACGGGAG‑CCCCCCCATCCGAAAACGTGGGGGGcac                                                                                                              >  2:632882/1‑151 (MQ=255)
                                      cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAG‑CCCCCACATCAGAAAGCGGGGGgcg                                                                                                                 >  2:871436/1‑147 (MQ=255)
                                              aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCCCATAAAAAAGCGGGGGGGGCACAcc                                                                                                           >  2:243222/1‑145 (MQ=255)
                                                   gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAg                                                                                                <  2:250054/151‑1 (MQ=255)
                                                        ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                                                                           <  2:796362/151‑1 (MQ=255)
                                                           tCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAATTACACACGGGAA‑CCCCCCCACCCAAAAAAGCGGTGCGCGCCCACAGCCGGGTAGCGaccat                                                                                         >  2:281116/1‑150 (MQ=255)
                                                               aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                                                                    <  1:487706/151‑1 (MQ=255)
                                                                   aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCCCGGGAAAAATTAAACACCGGAGACCCCCCCACCA‑AAAACGGGGGGGGCACAACGAGCGGCGGAGCGATCCTcccct                                                                                    >  1:1012113/1‑145 (MQ=255)
                                                                        cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATa                                                                             <  1:671134/149‑1 (MQ=255)
                                                                           aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAcc                                                                         <  2:159054/150‑1 (MQ=255)
                                                                                gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAG‑CCCCCACATCAGAAAAAGGGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGATAAACCaac                                                                      >  2:328714/1‑146 (MQ=255)
                                                                                    ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGGGGGAAAAATTACACACAGGAG‑CCCCCACATCAAAAAGAGGGGGGCGGAAACCCAGCTGCGGAGCGGACCTCCCCTTCTGGTAAAACAGCAcctg                                                                 >  1:470871/1‑146 (MQ=255)
                                                                                          aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACtt                                                         <  2:807188/151‑1 (MQ=255)
                                                                                           aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                        <  2:635653/151‑1 (MQ=255)
                                                                                               gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAG‑CCCCCACATCAGAAAACGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTTCTGATAAAACAGCACCGGCACCACTTCGCgg                                                    >  2:26542/1‑151 (MQ=255)
                                                                                               gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTACACACGGGAG‑CCCCCCCATCAGAAAACGTGGGGGGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACTGGCACCACCTCGCgg                                                    >  2:463124/1‑151 (MQ=255)
                                                                                                 gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGt                                                  <  1:869017/151‑1 (MQ=255)
                                                                                                        aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAG‑CCCCCACAACAGAAAGCGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAACACCGGCACCACCTCGCGGGGAtggggg                                           >  1:914531/1‑151 (MQ=255)
                                                                                                         agaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTACACACGGGAG‑CCCCCACATCAGAAAACGTGGTGCGCACACCCAGCTGCGGAGCGATCATCCCTTTCTGATAAAACAGCACTGGCACCACCTCGCGGGGAtggggg                                           >  2:564235/1‑150 (MQ=255)
                                                                                                          gaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCt                                         <  1:418766/151‑1 (MQ=255)
                                                                                                           aaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtg                                        <  1:689848/151‑1 (MQ=255)
                                                                                                            aaagaGCGACCATTTTGTGGCGGTGGGCCGCCGAAAAATTCAAACACCGGAA‑TCCCCCCATCTAAAAGGGTGGGGCGCACCCCCAGGCGCGTATCGGTCATCCCTTTCTGGTAAACCAACACCCGCCCCCCCTCCCGGGGAtggtgg                                           >  1:577420/1‑147 (MQ=255)
                                                                                                            aaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGTAG‑ACCCCACATCAGAAAGAGGGGTGGGCCCACCCCGCGGCGTCGCGGTCATCCCCTTCTGATAAAACAACACTGGCGCCACTTCGCGGGGATGGGGGGtg                                        >  1:249862/1‑150 (MQ=255)
                                                                                                                    aCCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATg                               <  1:200948/151‑1 (MQ=255)
                                                                                                                       aTTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc                            <  2:797243/151‑1 (MQ=255)
                                                                                                                       aTTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc                            <  2:899866/151‑1 (MQ=255)
                                                                                                                              ggcggtggcgCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                      <  1:143732/150‑1 (MQ=255)
                                                                                                                                  gtggcgCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGcc                 <  2:39989/151‑1 (MQ=255)
                                                                                                                                  gtggcgCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGcc                 <  1:204187/151‑1 (MQ=255)
                                                                                                                                          cgGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg         <  2:249021/151‑1 (MQ=255)
                                                                                                                                              aaaaTTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAt     >  1:198504/1‑151 (MQ=255)
                                                                                                                                                 aTTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAc  <  1:839324/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508516‑3508811

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: ‑