Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913908,7940AG74.5% 87.5 / 29.0 51G167G (GGT→GGCltaElow‑specificity L‑threonine aldolase
Reads supporting (aligned to +/- strand):  ref base A (6/7);  new base G (20/18);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCG  >  NC_000913/908656‑908926
                                                                                                                                          |                                                                                                                                    
gTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCAc                                                                                                                                     <  2:437103/140‑1 (MQ=255)
gTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCAc                                                                                                                                     <  1:942376/140‑1 (MQ=255)
   aTCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGt                                                                                                                                  <  1:842315/140‑1 (MQ=255)
    tCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTc                                                                                                                                 >  1:1006246/1‑140 (MQ=255)
    tCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTc                                                                                                                                 >  1:274831/1‑140 (MQ=255)
          ttACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATg                                                                                                                           <  2:696674/140‑1 (MQ=255)
          ttACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATg                                                                                                                           <  2:274831/140‑1 (MQ=255)
          ttACCGAAGAGAAATGCACCGACAGGCGTGCCAAGACCTTTCGACAGGGAAATGGTGAACGAATCACAATATTGTGTGATCTCTTTAAGTTCGCAGCCGTAAGCCGCCACGGCATTAAAGATGCGGGCGCCATCAACATg                                                                                                                           <  2:753511/140‑1 (MQ=255)
                          aaCCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTgcgc                                                                                                           >  1:286764/1‑140 (MQ=255)
                               aCTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTcgcg                                                                                                      >  2:804763/1‑140 (MQ=255)
                               aCTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTcgcg                                                                                                      >  2:557284/1‑140 (MQ=255)
                               aCTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTcgcg                                                                                                      >  2:339703/1‑140 (MQ=255)
                               aCTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTcgcg                                                                                                      >  1:709342/1‑140 (MQ=255)
                               aCTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTcgcg                                                                                                      >  1:802434/1‑140 (MQ=255)
                                         cAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTcc                                                                                            >  2:509611/1‑140 (MQ=255)
                                               cTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGct                                                                                      <  1:1118393/140‑1 (MQ=255)
                                                tttCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGctt                                                                                     <  2:802434/140‑1 (MQ=255)
                                                 ttCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGcttc                                                                                    >  1:1005949/1‑140 (MQ=255)
                                                  tCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGcttct                                                                                   <  2:1090808/140‑1 (MQ=255)
                                                   cGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGcttctt                                                                                  <  1:509611/140‑1 (MQ=255)
                                                        ggCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAgg                                                                             <  1:557284/140‑1 (MQ=255)
                                                         gCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGt                                                                            >  1:688970/1‑140 (MQ=255)
                                                              tGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTcc                                                                       <  2:1244197/140‑1 (MQ=255)
                                                                  gAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCcgcg                                                                   >  1:281783/1‑140 (MQ=255)
                                                                    aCGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGc                                                                 <  1:886932/140‑1 (MQ=255)
                                                                    aCGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGc                                                                 <  2:524557/140‑1 (MQ=255)
                                                                        aTCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAaca                                                             >  1:1150115/1‑140 (MQ=255)
                                                                                 ttGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGt                                                    >  1:664229/1‑140 (MQ=255)
                                                                                  tGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGtt                                                   >  1:248829/1‑140 (MQ=255)
                                                                                  tGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGtt                                                   >  1:281755/1‑140 (MQ=255)
                                                                                  tGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGtt                                                   >  1:331163/1‑140 (MQ=255)
                                                                                    cGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTtgt                                                 >  1:142376/1‑140 (MQ=255)
                                                                                        aTCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTGGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACCCTTTGCCGTTGTGGgt                                             <  1:804763/140‑1 (MQ=255)
                                                                                        aTCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGgt                                             <  2:1005949/140‑1 (MQ=255)
                                                                                        aTCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGgt                                             <  1:1150662/140‑1 (MQ=255)
                                                                                        aTCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGgt                                             <  2:248829/140‑1 (MQ=255)
                                                                                        aTCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGgt                                             <  2:142376/140‑1 (MQ=255)
                                                                                        aTCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGgt                                             <  2:1006246/140‑1 (MQ=255)
                                                                                           tctTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGtt                                          >  1:155411/1‑140 (MQ=255)
                                                                                             tttCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGtttt                                        >  1:1393356/1‑140 (MQ=255)
                                                                                               tCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGgt                                             >  2:389534/1‑133 (MQ=255)
                                                                                               tCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGgt                                             <  1:389534/133‑1 (MQ=255)
                                                                                                  gTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCAGCAACACTTTGCCGTTGTGGGTGTTTTCCAGa                                   <  2:281783/140‑1 (MQ=255)
                                                                                                   ttCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGAc                                  >  1:1145981/1‑140 (MQ=255)
                                                                                                    tCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACGCTTTGCCGTTGTGGGTGTTTTCCAGACt                                 <  2:274174/140‑1 (MQ=255)
                                                                                                      gCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGa                               <  1:416665/140‑1 (MQ=255)
                                                                                                       cAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAg                              >  1:1227268/1‑140 (MQ=255)
                                                                                                       cAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAg                              >  2:533615/1‑140 (MQ=255)
                                                                                                                    ccacGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCgg                 <  2:709342/140‑1 (MQ=255)
                                                                                                                         gCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGaa            <  1:533615/140‑1 (MQ=255)
                                                                                                                                  aTGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATAtc   >  1:755182/1‑140 (MQ=255)
                                                                                                                                   tGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATAtcg  >  1:372132/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                    
GTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCG  >  NC_000913/908656‑908926

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: