Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*U000962,947,5910A.23.7% 104.9 / 29.1 39intergenic (+18/‑16)metW/metVtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (21/8);  new base . (7/2);  total (28/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGC  >  U00096/2947491‑2947737
                                                                                                    |                                                                                                                                                  
cacggACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGa                                                                                                                                     <  2:1950940/112‑1 (MQ=255)
    gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaa                                                                                                                                                  <  1:1950974/99‑1 (MQ=39)
    gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAg                                                                                                                                             <  1:1950983/104‑1 (MQ=255)
    gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcgcg                                                                                                                                 <  1:1950968/116‑1 (MQ=255)
    gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcgc                                                                                                                                  <  1:1950976/115‑1 (MQ=255)
    gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAc                                                                                                                                    <  1:1950985/113‑1 (MQ=255)
                                 tcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGc                                                                                                             <  1:1950987/107‑1 (MQ=255)
                                                                     ggCACCCGCTATCAATTAAATTTGTAGTATGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtt                               <  2:1950999/149‑1 (MQ=14)
                                                                         cccGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtttttt                           >  2:1951102/1‑149 (MQ=255)
                                                                         cccGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtttttt                           >  2:1951038/1‑149 (MQ=255)
                                                                         cccGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtttttt                           >  1:1951113/1‑149 (MQ=255)
                                                                         cccGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtttttt                           >  2:1951118/1‑149 (MQ=255)
                                                                          ccGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttttt                          >  2:1951066/1‑149 (MQ=255)
                                                                            gCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCCAATTATTGAACACCCTAACGGGTGTTTTTTTGt                        >  2:1951088/1‑149 (MQ=255)
                                                                               aCCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTc                     >  1:1951110/1‑149 (MQ=255)
                                                                                ccAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCgg                                                                                                                                <  1:1950991/41‑1 (MQ=255)
                                                                                ccAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTc                     >  1:1951064/1‑148 (MQ=255)
                                                                                 cAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTg                   >  2:1951085/1‑149 (MQ=255)
                                                                                  aaTCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTgg                  >  1:1951086/1‑149 (MQ=255)
                                                                                   aTCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGt                 >  1:1951108/1‑149 (MQ=255)
                                                                                   aTCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGGGTTTTTTTGTTTCTGGt                 >  1:1951145/1‑149 (MQ=255)
                                                                                    taaaaTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCCTAAAGGGTGTTTTATTGTTTCTGGtct               >  2:1951025/3‑149 (MQ=17)
                                                                                    tCAAATTTGATGAAGTAAAAGCAGTACGGTGACGTGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtc                >  2:1950992/1‑149 (MQ=255)
                                                                                    tCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtc                >  2:1951134/1‑149 (MQ=255)
                                                                                     cAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtct               >  2:1951049/1‑149 (MQ=255)
                                                                                      aaaTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAATGAACACCCTCACGAGTGTTTTGTGGTATCTGGgctac             >  1:1951053/1‑147 (MQ=11)
                                                                                       aaTTTGATGAAGTAAAAGCAGTACGGTGACGCGTGGTGGAGCAGCATGGTAGCTCGTCGGGCTCATAAACCGAAGGTCGTCGGTTCAAATACGGCCACAGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTcc             >  2:1951084/1‑149 (MQ=21)
                                                                                        aTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAATGAACACCCAAACCGGTGTGTTTTTGTTTCGGGGCTCCCa           >  2:1951001/1‑149 (MQ=14)
                                                                                         tttGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCa           >  1:1951105/1‑149 (MQ=255)
                                                                                         tttGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTAGGGCTAATAACCCGAAGGTCTTCGGTTCAAATCCGGCCCCCGCACACAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCa           >  2:1951133/1‑149 (MQ=21)
                                                                                          ttGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGa                                                >  2:1951000/1‑110 (MQ=11)
                                                                                           tgatgaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtttttt                           <  1:1950996/130‑1 (MQ=255)
                                                                                           tgatgaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtttttt                           >  2:1950996/1‑130 (MQ=255)
                                                                                            gatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaa        >  1:1951109/1‑149 (MQ=255)
                                                                                            gatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaa        <  2:2434574/149‑1 (MQ=255)
                                                                                             atgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaaa       >  1:1951083/1‑149 (MQ=255)
                                                                                             atgaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAaca                                             <  2:1950993/110‑1 (MQ=18)
                                                                                             atgaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCATTACGGGTGTGTTTTTGTTTCAGGTCTCCAATaaca      >  2:1951052/1‑147 (MQ=17)
                                                                                              tgaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAATTAACACCCTAACAGGGGATTATTTGTATCGGGTCTCCCAaaaaaa     >  2:1951062/1‑149 (MQ=14)
                                                                                                  gTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgc  >  2:1951031/1‑149 (MQ=255)
                                                                                                    aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATaa                                                                                              <  1:1950997/55‑1 (MQ=31)
                                                                                                    |                                                                                                                                                  
CGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGC  >  U00096/2947491‑2947737

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