Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,037,4960AT41.9% 5.3 / 33.9 32intergenic (+161/‑23)alaT/rrlAtRNA‑Ala/23S ribosomal RNA of rrnA operon
Reads supporting (aligned to +/- strand):  ref base A (11/7);  new base T (6/7);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTAC  >  NC_000913/4037355‑4037538
                                                                                                                                             |                                          
aaaaaTACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaa                                           <  1:446116/143‑1 (MQ=38)
aaaaaTACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaa                                           <  1:447537/143‑1 (MQ=38)
 aaaaTACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgtaa                                          >  2:296230/1‑143 (MQ=25)
    aTACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaa                                       <  1:606276/143‑1 (MQ=39)
    aTACGTCAGAGTGTACCTGCAAAGGTTAACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCAGAACAAGTAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaa                                       <  1:102316/143‑1 (MQ=14)
     tACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAACTTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGCATCgaaagaac                                      >  2:53798/1‑142 (MQ=21)
      aCTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAAAAATGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAAc                                     <  2:673569/143‑1 (MQ=255)
        ttCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACAt                                   <  2:447458/143‑1 (MQ=25)
        ttCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAAGCGAAAGAAACAt                                   >  2:663612/1‑143 (MQ=38)
         tCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTATCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATc                                  >  2:705776/1‑143 (MQ=25)
              gtgtACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAGATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTGCGCAACACGATGATGAATCGTAAGAAACATCTTCgg                             <  2:684794/143‑1 (MQ=255)
               tgtACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGCGGCTGTCAAAGTGTCGCAACACGATGATGAATCGAAAGAAACGTGTTCggg                            >  1:197958/1‑143 (MQ=11)
               tgtACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCggg                            <  1:477901/143‑1 (MQ=39)
                   ccTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTtgt                        >  1:89338/1‑143 (MQ=37)
                    cTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGCACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGGAACATCTTCGGGTtgtg                       >  1:271823/1‑143 (MQ=255)
                      gCAAAGGTTCACTGCGAAGTTTTGGTCTGTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAGGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGCAAGAACCATCTTCGGGTTGTGAg                     >  2:442932/1‑143 (MQ=11)
                      gCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGACAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCACATTTTCGCAACACGATGATGAATCGTAGGAAACATCTTCGGGTTGTGAg                     >  1:8811/1‑143 (MQ=255)
                      gCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAGATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAg                     <  2:237436/143‑1 (MQ=255)
                      gCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAg                     >  2:544325/1‑143 (MQ=21)
                      gCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAg                     >  1:681918/1‑143 (MQ=39)
                      gCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTGTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAg                     <  2:480030/143‑1 (MQ=25)
                         aaGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGtt                  >  1:624307/1‑143 (MQ=25)
                         aaGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGtt                  <  2:404601/143‑1 (MQ=21)
                           ggTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTaa                <  2:384838/143‑1 (MQ=38)
                           ggTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTaa                <  1:270566/143‑1 (MQ=38)
                               cACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGa            <  1:566027/143‑1 (MQ=38)
                               cACTGCGAAGTGTTGCTCTTTAAAAATCTGGATCAAGCTGAAAAATGACACTCCATACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGa            <  1:722749/143‑1 (MQ=11)
                                  tGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTGAGCGACTa         >  2:659591/1‑143 (MQ=25)
                                   gCGAAGTTTTGCTCTTTAAAAATCTGGTTCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTaa        <  1:703706/143‑1 (MQ=17)
                                   gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGGGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTaa        <  2:540224/143‑1 (MQ=25)
                                    cGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAg       >  2:188348/1‑143 (MQ=17)
                                    cGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAg       >  1:721521/1‑143 (MQ=34)
                                    cGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTAGACCGATGAGGTTAAGCGACTAAg       >  2:666017/1‑143 (MQ=11)
                                     gAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAGACATCTTCGGGTTGTGAGGTTAAGCGACTAAGc      >  1:6379/1‑143 (MQ=18)
                                      aaGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAACTTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg     <  2:219629/143‑1 (MQ=17)
                                      aaGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg     >  2:91682/1‑143 (MQ=32)
                                         ttttGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACACCGAAAGTTGTTCGTGAGTCTCTCCAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGAGAATAAGCgggc  >  2:538468/1‑140 (MQ=11)
                                         ttttGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTac  >  1:580469/1‑143 (MQ=11)
                                                                                                                                             |                                          
AAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTAC  >  NC_000913/4037355‑4037538

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