Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,130,1910CA62.2% ‑2.6 / 35.2 46intergenic (‑240/+64)gcvH‑I/PP_0990glycine cleavage system protein H/DoxX family protein
Reads supporting (aligned to +/- strand):  ref base C (6/11);  new base A (28/0);  total (35/11)
Fisher's exact test for biased strand distribution p-value = 1.22e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGA‑CGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCC  >  NC_002947/1130057‑1130325
                                                                                                                                      |                                                                                                                                       
cAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCg                                                                                                                                        <  1:2310829/136‑1 (MQ=255)
         aCGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTTCGA‑CGt                                                                                                                              >  2:162670/1‑136 (MQ=255)
           gTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGA‑CGTAg                                                                                                                            >  2:1841168/1‑136 (MQ=255)
            tAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGA‑CGttgg                                                                                                                           >  1:92974/1‑133 (MQ=255)
                       agaAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGA‑CGTAGCGAAAACGATaa                                                                                                                >  2:1633227/1‑135 (MQ=255)
                                    attGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGA‑CGTAGCGAAAACGATACATTTGGTGCCTgg                                                                                                   >  2:2488606/1‑136 (MQ=255)
                                         gtgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGG‑CGTAGGGGAAAAGATACATTTGGTGgctgggctgg                                                                                              >  1:2984627/1‑136 (MQ=255)
                                          tgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACACTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGGTAAGA‑CCATGCGGAAACGATAAAATTGGTGACAGGACTGGc                                                                                             >  1:2832296/1‑136 (MQ=255)
                                          gggAATTATACCGTACCCGGTGCCCCTGACAGCAAACGGGTTTCCGCACCAAACCCTTGTGCCCCGACGGCGAGGCGGCGGTGTATCGTTTCCCCTACCA‑CGTAGCTAAAACGATACATGTGGTGCCTGGACTGGc                                                                                             <  2:1716881/135‑1 (MQ=255)
                                           gtaaTTATATCGTTCCCGGGCCCTCTGAGAGCAAACGGGTTGCCCCACCAACCCGTTGTCCCCCGACGGCGTGGCGGCTGTGTATCGTTTCCCCACCCC‑CGTAGCTAAAACGATACATGTGGTGCCTGGACTGGCg                                                                                            <  2:275338/134‑1 (MQ=255)
                                             aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGC‑CGTAGCGAAAAAGATACATGTGGTGCCTGGACTGGCGAg                                                                                          >  1:104378/1‑136 (MQ=255)
                                             aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGA‑CGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAg                                                                                          >  1:2121017/1‑136 (MQ=255)
                                             aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCAAACA‑CGTATCGAAAACCAAAAATTTGGTGGCATGACAGGCGAg                                                                                          >  1:2232149/1‑136 (MQ=255)
                                                 atctCGTTACCGGCGCCTCTGACAGGAATAGGGTTGCACCATCAAACCCTTGTGCACCGACGGCGCGTTGGCGGGGTATCGTTTTCGCTCCGC‑CGTAGCTAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGc                                                                                      <  2:2894274/133‑1 (MQ=255)
                                                    tCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGA‑CGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCaa                                                                                   >  2:2387987/1‑136 (MQ=255)
                                                     cGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTTCCA‑CGTAGCGAAAACGATAAATGTGGTGCCTGGACTGGCGAGGTGCCaaa                                                                                  >  1:1236039/1‑135 (MQ=255)
                                                       ttACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTAAGA‑CGTAGGGAAAACGATACATGTGGTGCCTGGACTGGCGTGGTGGCAATgg                                                                                >  2:783100/1‑136 (MQ=255)
                                                       ttACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGA‑CGTAGCGAAAACGATACAAGTGGTGGCTGGACTGGCGAGGTGCCAATgg                                                                                >  2:2403934/1‑136 (MQ=255)
                                                         aCAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACAGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGA‑AGGATCGAAAACCATACAATTGGGGCCTGGAATGGCGAGGTGACAATgggc                                                                              >  1:746651/1‑134 (MQ=255)
                                                          cAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGA‑CGTAGCGAAAACGAGACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGccg                                                                             >  2:690102/1‑135 (MQ=255)
                                                           cggTGCCTCTGACAGGAAAAGGGTTGCCCCAACAACCCCTTGTCCACCTACGGCCAGCCGGCGGTGAAGCTTATTCGCAACCA‑CCTACCTAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCtt                                                                            <  2:190083/135‑1 (MQ=255)
                                                            ggTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTACGTACCACCGTA‑AGAAAACCATAAAATTGGTTCCAGGGCTGGGCACGTGCCAATGGcccgg                                                                           >  1:1532429/1‑133 (MQ=255)
                                                              tGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTAACG‑CGTAGCGAAAACCATAAATGTGGTTCCAGGGCTGGCGAGGTGCCAATTGCCtaggg                                                                         >  1:1973160/1‑132 (MQ=255)
                                                                   ctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTTGGTTCGG‑CCGAGCGAAAACCAAAAAAGTGGTGGCGGGCCGGGCCAGGGGCCAAAGGCCGTAGGTCAta                                                                    >  1:233487/1‑135 (MQ=255)
                                                                   ctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGA‑CGTAGCGAAAACGATACATGTGGTGCCTGGTCTGGCGAGGTTCCAATGGCCTTAGGtactg                                                                    >  1:1281441/1‑132 (MQ=255)
                                                                   ctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCGTAGA‑CCGTGCGAAAAAGATAAATGTGGTGGCAGGACAGGCGCGGTGGCAATGGCCATCGGGCAtg                                                                    >  1:1581388/1‑135 (MQ=255)
                                                                    tGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTATGTACGA‑CGTAGCGAAAACGATAAAATTGGTGGCAGGACAGGCGGGGTGGCAAAGGCCGTGGGTGAtgt                                                                   >  1:3158097/1‑134 (MQ=255)
                                                                          gAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGA‑CGGAGCGAAAACGATACAAGTGGTGGCGGGACTTGCGAGGTGCCAAAGGCCTTAGGTCATCTCACatt                                                             >  2:1879599/1‑136 (MQ=255)
                                                                                ggTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCGACGG‑CGGAGCGAAAACGATACATGTGGTGGCGGGACTGGCGGGGTGGCAATGGGCTTCGGTGATGTGACTTTcagcag                                                       >  2:2156967/1‑136 (MQ=255)
                                                                                 gTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCTTTTTCGCTACGA‑CGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAga                                                      <  2:1191222/136‑1 (MQ=255)
                                                                                     gAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCGTCGA‑CGGAGCGAAAACGATTCATGTGGTGGCTGGGCTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAgatgct                                                  >  2:2661509/1‑134 (MQ=255)
                                                                                             aaCCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTGCGA‑CGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGaa                                          >  1:1871636/1‑136 (MQ=255)
                                                                                                    gtgCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGA‑CGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACgcgc                                   <  2:2121017/136‑1 (MQ=255)
                                                                                                    gtgCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGA‑CATAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACgcgc                                   <  2:335088/136‑1 (MQ=255)
                                                                                                     tgCACCGACGGCGAGGAGGCTGTGTATCGTTTTACCTACCA‑CGTATCGAAAACGAAACAAGTGGTGCCAGGACTGGCGAGGTGGCAATGGCCTTAGGTAATCACACATTAAGCAGAAGAAAACGCGAAAACGCGcc                                  >  1:1269408/1‑136 (MQ=255)
                                                                                                      gCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCGACGA‑CGTAGCGCAAACGATACAAGTGGTGCCTGGACCGGCGGGGTGCCAATGGCCTTCGGTGATCGGACTTTCAGCAGATGATCCAGCGTAAACGCGCCg                                 >  1:221976/1‑136 (MQ=255)
                                                                                                             cGGCGAGGAGGCTGTGTATCGTTTTAGCTACGA‑CGTAGCGAAAACGATACATGTGGTGCCGGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTc                          >  1:3036468/1‑136 (MQ=255)
                                                                                                               gCGAGGAGGCTGTGTATCGTTTTAGCTACGC‑CGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGGCTTCGGTCATCTCACTTTGAGCAGATGATCCAGCGAAAACGCGCCGGCGCCACgg                        >  1:94257/1‑136 (MQ=255)
                                                                                                                  aggaggCTGTGTATCGGTTTAGGGAAGA‑CGGTGCGGAAACGCGAAAAGTGGTGGCGGGGCAGGGGGGGTGGCAAATGGCGTGGGTGATCTGACCTTGAGCAGATTTTACACCGAAAACGCGCCGGCGCCGCggggc                     >  2:1473241/1‑134 (MQ=255)
                                                                                                                     aggCTGTGTATCGTTTTAGCTACGA‑CGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAAcc                  >  1:2400920/1‑136 (MQ=255)
                                                                                                                        gtgtgTATCGTTTTCGCTACCC‑CGTAGCTAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAgg               <  2:3036468/135‑1 (MQ=255)
                                                                                                                        gtgtgTAGCTTTTCCGCTACGA‑CGTAGCTAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAgg               <  2:104378/135‑1 (MQ=255)
                                                                                                                          gtgtATCGTTTTCGCTACGA‑CGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGCGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGag             >  2:2087077/1‑136 (MQ=255)
                                                                                                                               tCGTTTTAGCTACGG‑CGTAGCGCAAACGAGACAAGTGGTGCCTGGTCTGGGGCGGTGGCAATGGCCGTGGGTCATCGCCCTTTCAGCAGCTGATCCAGCGAAAACGCGCCGGCGCCGCGGGCCCCCAGGagcagcc        >  2:632470/1‑135 (MQ=255)
                                                                                                                                  ttttCGCTTCGC‑GTTAGCTAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAaga     <  2:746651/136‑1 (MQ=255)
                                                                                                                                     tCGCTACGA‑CGTAGCGAAAACGATTCATGTGGTGCCTGGGCTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGcc  >  2:1606928/1‑136 (MQ=255)
                                                                                                                                      |                                                                                                                                       
CAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGA‑CGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCC  >  NC_002947/1130057‑1130325

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