Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,530,1730AC33.3% 20.1 / 15.7 18*399S (TAA→TCA) 
noncoding (62/244 nt)
ftsZ
PP_mr19
Cell division protein FtsZ
Xoo8 ncRNA
Reads supporting (aligned to +/- strand):  ref base A (4/8);  new base C (3/3);  total (7/11)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCC  >  AE015451/1529963‑1530375
                                                                                                                                                                                                                  |                                                                                                                                                                                                          
agaCCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGt                                                                                                                                                                                  <  1:226124/235‑1 (MQ=255)
         aGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTAt                                                                                                                                                                         >  1:167806/1‑237 (MQ=255)
           aCTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCTGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATAAGGGGTATAAAGGt                                                                                                                                                                                  >  2:403829/1‑226 (MQ=255)
                                      aaTCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTa                                                                                                                                                                                                            >  1:238649/1‑173 (MQ=255)
                                      aaTCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTa                                                                                                                                                                                                            <  2:238648/173‑1 (MQ=255)
                                           tGCGCCGGTGCGCCAGGGGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCGCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCa                                                                                                                                          <  1:163398/234‑1 (MQ=255)
                                                     cgcCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGt                                                                                                                                                                                  <  1:139697/184‑1 (MQ=255)
                                                                                     ccTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTCCCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATACAGGTGATTGGTATTCATCACAGGTTGGTGCTGtntt                                                                                                                                                  >  2:53393/1‑181 (MQ=255)
                                                                                      cTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTTCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGc                                                                                                                    <  2:167805/213‑1 (MQ=255)
                                                                                        ggAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTatcatc                                                                                                                                              <  1:93577/185‑1 (MQ=255)
                                                                                                       cATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCACTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGc                                                                            <  1:206897/236‑1 (MQ=255)
                                                                                                                ccAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTAt                                                                                                                                                                                                 >  2:65425/1‑110 (MQ=255)
                                                                                                                ccAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTAt                                                                                                                                                                                                 <  1:65426/110‑1 (MQ=255)
                                                                                                                                          ccGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATAAAAGGTTg                                                                                                                                                            <  1:403831/121‑1 (MQ=255)
                                                                                                                                                                               cGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGt                                                                                                                                                                                  <  1:48611/62‑1 (MQ=255)
                                                                                                                                                                               cGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGt                                                                                                                                                                                  >  2:48611/1‑62 (MQ=255)
                                                                                                                                                                               cGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGc   >  1:325271/1‑237 (MQ=255)
                                                                                                                                                                                 atatCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGcc  >  2:215194/1‑236 (MQ=255)
                                                                                                                                                                                          cATTCCTGCGTCGTCAGTCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGc   <  1:365083/226‑1 (MQ=255)
                                                                                                                                                                                                  cgtcgtCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCa                                   <  1:24409/186‑1 (MQ=255)
                                                                                                                                                                                                                  |                                                                                                                                                                                                          
ATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCC  >  AE015451/1529963‑1530375

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