Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,138,3500TG50.0% 8.2 / 20.7 36F240V (TTT→GTT) rbsC_2Ribose import permease protein RbsC
Reads supporting (aligned to +/- strand):  ref base T (8/10);  new base G (18/0);  total (26/10)
Fisher's exact test for biased strand distribution p-value = 3.44e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATTCGGCCTGCGTCTGTACGCCATCGGAGAATACCCGGAGGCGGCTAATGCCGCTGGCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGTTTTGTTCGGCGGCCTGGTTCAGCGGCAGCACCACCGGTTCCGGCGAGATGCTGCTGTCCGTGGTGGCCATCGCTTTTCTTGGGGTGATCTTCTCCCAGCGTCTGCAGCCGAGCATCGGCGGCACGCTGCT  >  ACB143/3138226‑3138480
                                                                                                                            |                                                                                                                                   
aTTCGGCCTGCGTCTGTACGCCATCGGAGAATACCCGGAGGCGGCTAATGCCGCTGGCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGGTTTGTTgggcg                                                                                                                          >  5:145892/1‑136 (MQ=255)
aTTCGGCCTGCGTCTGTACGCCATCGGAGAATACCCGGAGGCGGCTAATGCCGCTGGCCTGTCGGAGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGTgttgttatggg                                                                                                                          >  3:225656/1‑131 (MQ=255)
   cGGCCTGCGTCTGTACGCCATCGGAGAATACCCGGAGGCGGCTAATGCCGCTGGCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGGTTTGTTCGGCGGcc                                                                                                                       >  2:45278/1‑136 (MQ=255)
   cGGCCTGCGTATGTACGCCATCGGAGAATACCCGGAGGCGGCTAATGCCGCTGGCCTGTCGGTGCGTGGCTACGTGTTCGCCCGCTACCTGATCAGCGGTCTGTGCGGAAGCCTGGCCGCGGGTTGTGTGGCGGcc                                                                                                                       >  6:152873/1‑136 (MQ=255)
     gCCTGCGTCTGTACGCCATCGGAGAATACCCGGAGGCGGCTAATGCCGCTGGCCTGTCGGTGCGTGGCTACGTGTTCGCCCGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGGTTTTTTAggcggcgag                                                                                                                     >  5:190273/1‑133 (MQ=255)
     gCCTGCGTCTGTACGCCATCGGAGAATACCCGGAGGCGGCTAATGCCGCTGGCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGTTTTTTTCGGCGGCCTg                                                                                                                     >  8:157860/1‑136 (MQ=255)
     gCCTGCGTCTGTACGCCATCGGAGAATACCCGGAGGCGGCTAATGCCGCTGGCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGGTTTTTTGGGGGGccgg                                                                                                                     >  4:198992/1‑134 (MQ=255)
            tCTGTACGCCATCGGAGAATACCCGGAGGCGGCTAATGCCGCTGGCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGGTTTGTTGGGGGGCCGGGGTCAgg                                                                                                              >  3:18311/1‑135 (MQ=255)
            tCTGTACGCCATCGGAGAATACCCGGAGGCGGCTAATGCCGCTGGCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGGTTTGTTAGGCGGGCGGGTTTAGc                                                                                                              >  3:323497/1‑136 (MQ=255)
                   gCCATCGGAGAATACCCGGAGGCGGCTAATGCCGCTGGCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGGTTTGTTGGGCGGCCGGGTTGAGCGgcagca                                                                                                      >  1:155996/1‑136 (MQ=255)
                                  ccGGAGGCGGCTAATGCCGCTGGCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGGGTTGGTCGGCGGCCTGGGTCAGCGGCAGCACCCCCCGGTgcggcg                                                                                       >  4:85859/1‑136 (MQ=255)
                                          ggcTAATGCCGCTGGCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGTTTTGTTCGGCGGCCTGGTTCAGCGGCAGCACCACCGGTTCCGGCGAGAtgctg                                                                               >  6:179235/1‑136 (MQ=255)
                                          ggcTAATGCCGCTGGCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGTTGTTTTAGGGGGCCGGGGTCAGCGGCAGCACCCCCCGCTCCGGCgagaggaag                                                                               >  3:256279/1‑131 (MQ=255)
                                            cTAATGCCGCTGGCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGGTTTGTTGGGCGGCCGGGTTGAGCGGCCGCCCCCCCCGGTGCGGCGCGAtgctggt                                                                             >  7:130911/1‑134 (MQ=255)
                                                     cTGGCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGTTTTGTTCGGCGGCCTGGTTGAGCGGCAGCACCACCGGTTCCGGCGAGATGCTGCTGTCCgtggt                                                                    <  1:324589/136‑1 (MQ=255)
                                                        gCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGTTTTGTTCGGCGGCCTGGTTCAGCGGCAGCCCCCCCGGTTCCGGCGAGATGCTGCTGTGCGGGGTGGc                                                                 >  1:199252/1‑136 (MQ=255)
                                                        gCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGTTTTGTTAGGCGGCCGGGTTCAGCGGCAGCCACCCCCGGTCGGGCGAGATGCAGCTTTCAGTGGTGGc                                                                 >  6:255771/1‑136 (MQ=255)
                                                        gCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGGTTTGTTCGGCGGGCTGGTTCAGCGGCAGCACCACCGGTTCCGGCGAGATGCTGCTGTCCGTGGTGGc                                                                 >  1:43610/1‑136 (MQ=255)
                                                        gCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGGTTTGTTCGGCGGCCTGGGTCAGCGGCAGCACCACCCGTTCCCGCGAGATGCTGCAGTCCGGGGTGGc                                                                 >  6:99319/1‑136 (MQ=255)
                                                             tCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGTTTTGTTCGGCGGCCTGGTTCAGCGGCAGCACCACCGGTTCCGGCGAGATGCTGCTGTCCGTGGTGGCCATc                                                             <  4:370950/135‑1 (MQ=255)
                                                             tCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGGCTGTGCGGCAGCCTGGCCGCGTTTTGTTCGGCGGCCTGGTTCAGCGGCAGCACCACCGGTTCCGGCGAGATGCTGCTGTCCGTGGTGGCCAt                                                              <  3:38350/134‑1 (MQ=255)
                                                                gTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGTTTTGTTCGGCGGCCTGGTTCAGCGGCAGCACCACCGGTTCCGGCGAGATGCTGCTGTCCGTGGTGGCCATCGCtt                                                         <  5:88576/136‑1 (MQ=255)
                                                                     tGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGTTTTGTTCGGCGGCCTGGTTCAGCGGCAGCACCACCGGTTCCGGCGAGATGCTGCTGTCCGTGGTGGCCATCGCTTTTCtt                                                    <  3:39034/136‑1 (MQ=255)
                                                                     tGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGGTTTTTTAGGGGGGCGGGGTGAGCGGCAGCACCACCGGGTACGGGGAGGTGCTGGGGTGCGTGGTGGGCGTAGCTtttttt                                                    >  5:266795/1‑136 (MQ=255)
                                                                      ggCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGGTTTTTTGGGCGGGCGGGTTCAGCGGCCGCCCCACCGGTTGCGGCGAGCTGCGGCTGTCCGGGGTGGGCATCGCTTTTTTTg                                                   >  6:189660/1‑136 (MQ=255)
                                                                      ggCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGGTTTGGTCGGCGGCCAGGGTGAGCGGCCTCACACCCCGGTGCGGGGGGATGCTGGTGGGCGGGGGGGGCCTCGGGTTTCTTg                                                   >  3:104048/1‑136 (MQ=255)
                                                                       gCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGTTTTGTTCGGCGGCCTGGTTCAGCGGCAGCACCACCGGTTCCGGCGAGATGCTGCTGTCCGTGGTGGCCATCGCTTTTCTTgg                                                  <  2:305136/136‑1 (MQ=255)
                                                                          aCGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGGTTTGTTGGGCGGGCGGGGTGAGCGGCCGCACCCCCGGGTGCGGCGGGATGCTGGTGTCCGTGGTGGGCCTCGCTTTTTGTggggg                                               >  1:229125/1‑135 (MQ=255)
                                                                          aCGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGGTTTGTTCGGCGGCCTGGGTGAGCGGCAGCACCCCCGGTTCCGGCGAGATGGTGGTGTCCGTGGTGGGCATCGCGTTTGTTGGGGt                                               >  5:178463/1‑136 (MQ=255)
                                                                           cGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGTTTTGTTCGGCGGCCTGGTTCAGCGGCAGCACCACCGGTTCCGGCGAGATGCTGCTGTCCGTGGTGGCCATCGCTTTTCTTGGGGTg                                              <  3:86283/136‑1 (MQ=255)
                                                                              gtTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGGTTTTTTAGGCCGCCAGGTTAAGCCGCACCACCACCCGGTACGGCGAGATGCGTCTGTCCGGGGTTGCGCTAGCTATTATTTGGGTGTtc                                           >  2:329161/1‑136 (MQ=255)
                                                                                          ccTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGTTTTGTTCGGCGGCCTGGTTCAGCGGCAGCACCACCGGTTCCGGCGAGATGCTGCTGTCCGTGGTGGCCATCGCTTTTCTTGGGGTGATCTTCTCCCAGc                                 <  8:278937/134‑1 (MQ=255)
                                                                                          ccTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGTTTTGTTCGGCGGCCTGGTTCAGCGGCAGCACCACCGGTTCCGGCGAGATGCTGCTGTCCGTGGTGGCCATCGCTTTTCTTGGGGTGATCTTCTCCCAGc                                 <  4:225656/134‑1 (MQ=255)
                                                                                                      tCTGTGCGGCAGCCTGGCCGCGTTTTGTTGGGCGGGCTGGTTAAGCGGCAGCACCACCGGTTGCGGCGAGATGCTGCTGTACGTGGTGGCCATCGCTTTTGTTGGGGTGTTCTTCTCCCCGCGTCTGCAGCCGAg                    >  1:111854/1‑135 (MQ=255)
                                                                                                              gCAGCCTGGCCGCGTTTTGTTCGGCGGCCTGGTTCAGCGGCAGCACCACCGGTTCCGGCGAGATGCTGCTGTCCGTGGTGGCCATCGCTTTTCTTGGGGTGATCTTCTCCCAGCGTCTGCAGCCGAGCATcggcgg           <  7:138639/136‑1 (MQ=255)
                                                                                                                       ccGCGTTTTGTTCGGCGGCCGGGGTGAGCGGCCGCACCACCGGTTGCGGCGAGATTCTGGTGTGCGTGGTGGGCCTCGGGGTTGTTGGGGGGGTGTTATACCAGCGTCTGCCGCCCAGCATAGGCGGCCCgcagag  >  6:250676/1‑132 (MQ=255)
                                                                                                                            |                                                                                                                                   
ATTCGGCCTGCGTCTGTACGCCATCGGAGAATACCCGGAGGCGGCTAATGCCGCTGGCCTGTCGGTGCGTGGCTACGTGTTCGCCAGCTACCTGATCAGCGGTCTGTGCGGCAGCCTGGCCGCGTTTTGTTCGGCGGCCTGGTTCAGCGGCAGCACCACCGGTTCCGGCGAGATGCTGCTGTCCGTGGTGGCCATCGCTTTTCTTGGGGTGATCTTCTCCCAGCGTCTGCAGCCGAGCATCGGCGGCACGCTGCT  >  ACB143/3138226‑3138480

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