Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1351,4920GA50.0% 4.0 / 29.6 30intergenic (‑6569/+60384)unknown/unknowntRNA‑Ser‑GCT/LSU rRNA ## 23S rRNA, large subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (4/11);  new base A (7/8);  total (11/19)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
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.

GTGTCTAAACTTTATAAGATTCCTCATTCTCAAGTCATTTATTGGTTAAAAGCCTTTCGTGAAAGAGGGCTTGATGGCATAAAATCGCCTTATACTAAACCTAAAACCTCCAAAATTCTGAAACCCAAAATGAAAAAAACGCAGATTGAAGTGGCTGAAACAACCGACTTTTCACCACAAGCCTTTAAAAAACTGAAACGAGAACCCGCCCTTGAAAAAGCGGAGATTGATTATCTAAAAAAGTTAGTGGAATTGGGGCGATCAAAAACAAATTAAACTCGATTTAAAAG  >  1/351448‑351737
                                            |                                                                                                                                                                                                                                                     
gtggctaaaCTTCATCAGATTCCTCATTCTCAAGTCATTTATTGATTAAAAGCCTTTCGTGAAAGAGGACTTGATGGCGTAAAATCACCTTATACTAAACCTAAAAtacccaaaatgaaaaaaatgcagattgaagtgactgaaacaaccgacttttcaccacaagccttgaaaaaactcaaacgagaagtcataggcgataatttcgccattgaagagatctttaattggtgataaatagagcttgcc                                           <  1:295425‑M1/245‑144 (MQ=255)
  ggctaAACTTCATCAGATTCCTCATTCTCAAGTCATTTATTGATTAAAAGCCTTTCGTGAAAGAGGACTTGATGGCGTAAAATCACCTTATACTAAACCTAAAAtacccaaaatgaaaaaaatgcagattgaagtgactgaaacaaccgacttttcaccacaagccttgaaaaaactcaaacgagaagtcataggcgataatttcgccattgaagagatctttaattggtgataaatagagc                                                <  1:313133‑M1/240‑139 (MQ=255)
    cTAAACTTCATCAGATTCCTCATTCTCAAGTCATTTATTGATTAAAAGCCTTTCGTGAAAGAGGACTTGATGGCGTAAAATCACCTTATACTAAACCTAAAAtacccaaaatgaaaaaaatgcagattgaagtgactgaaacaaccgacttttcaccacaagccttgaaaaaactcaaacgagaagtcataggcgataatttcgccattgaagagatctttaattggtgataaatagagcttgccttcat                                      <  1:85033‑M1/250‑149 (MQ=255)
    cTAAACTTCATCAGATTCCTCATTCTCAAGTCATTTATTGATTAAAAGCCTTTCGTGAAAGAGGACTTGATGGCGTAAAATCACCTTATACTAAACCTAAAAtacccaaaatgaaaaaaatgcagattgaagtgactgaaacaaccgacttttcaccacaagccttgaaaaaactcaaacgagaagtcataggcgataatttcgccattgaagagatctttaattggtgataaatagagcttgccttcat                                      <  1:327219‑M1/250‑149 (MQ=255)
       aaCTTCATCAGATTCCTCATTCTCAAGTCATTTATTGATTAAAAGCCTTTCGTGAAAGAGGACTTGATGGCGTAAAATCACCTTATACTAAACCTAAAAtacccaaaatgaaaaaaatgcagattgaagtgactgaaacaaccgacttttcaccacaagccttgaaaaaactcaaacgagaagtcataggcgataatttcgccattgaagagatctttaattggtgataaatagagcttgccttcatcgca                                  >  1:48938‑M1/1‑99 (MQ=255)
           tcatcagattcCTCATTCTCAAGTCATTTATTGATTAAAAGCCTTTCGTGAAAGAGGACTTGATGGCGTAAAATCACCTTATACTAAACCTAAAAtacccaaaatgaaaaaaatgcagattgaagtgactgaaacaaccgacttttcaccacaagccttgaaaaaactcaaacgagaagtcataggcgataatttcgccattgaagagatctttaattggtgataaatagagcttgccttcatcgcattta                              >  1:124268‑M1/6‑95 (MQ=255)
           tcatcagattcCTCATTCTCAAGTCATTTATTGATTAAAAGCCTTTCGTGAAAGAGGACTTGATGGCGTAAAATCACCTTATACTAAACCTAAAAtacccaaaatgaaaaaaatgcagattgaagtgactgaaacaaccgacttttcaccacaagccttgaaaaaactcaaacgagaagtcataggcgataattt                                                                                      >  1:241022‑M1/6‑95 (MQ=255)
              ctaaatttctcgttCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTGATGGCGTAAAATCACCTTACAC‑AAACCTAAAACCTCCAAAATTCTGAAACCCAAAATGAAAAAAACGCAGATTGAAGTGGCTGAAACAACCGACTTTTCACCACAAGCCTTTAAAAAACTGAAACGAGAACCCGCCCTTGAAAAAGCGGAGATTGATTATCTAAAAAAGTTAGTGGAATTGGGGCGATCaaaaa                        >  1:660182/13‑251 (MQ=11)
               taaatttctcgttCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTGATGGCGTAAAATCACCTTACAC‑AAACCTAAAACCTCCAAAATTCTGAAACCCAAAATGAAAAAAACGCAGATTGAAGTGGCTGAAACAACCGACTTTTCACCACAAGCCTTTAAAAAACTGAAACGAGAACCCGCCCTTGAAAAAGCGGAGATTGATTATCTAAAAAAGTTAGTGGAATTGGGGCGATCAAAAAc                       <  1:500863/240‑1 (MQ=11)
               taaatttctcgttCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTGATGGCGTAAAATCACCTTACAC‑AAACCTAAAACCTCCAAAATTCTGAAACCCAAAATGAAAAAAACGCAGATTGAAGTGGCTGAAACAACCGACTTTTCACCACAAGCCTTTAAAAAACTGAAACGAGAACCCGCCCTTGAAAAAGCGGAGATTGATTATCTAAAAAAGTTAGTGGAATTGGGGCGATCAAAAAc                       <  1:10064/240‑1 (MQ=11)
               taaatttctcgttCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTGATGGCGTAAAATCACCTTACAC‑AAACCTAAAACCTCCAAAATTCTGAAACCCAAAATGAAAAAAACGCAGATTGAAGTGGCTGAAACAACCGACTTTTCACCACAAGCCTTTAAAAAACTGAAACGAGAACCCGCCCTTGAAAAAGCGGAGATTGATTATCTAAAAAAGTTAGTGGAATTGGGGCGATCAAAAAc                       <  1:777479/240‑1 (MQ=11)
               taaatttctcgttCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTGATGGCGTAAAATCACCTTACAC‑AAACCTAAAACCTCCAAAATTCTGAAACCCAAAATGAAAAAAACGCAGATTGAAGTGGCTGAAACAACCGACTTTTCACCACAAGCCTTTAAAAAACTGAAACGAGAACCCGCCCTTGAAAAAGCGGAGATTGATTATCTAAAAAAGTTAGTGGAATTGGGGCGATCAAAAAc                       <  1:749024/240‑1 (MQ=11)
                aaatttctcgttCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTGATGGCGTAAAATCACCTTACAC‑AAACCTAAAACCTCCAAAATTCTGAAACCCAAAATGAAAAAAACGCAGATTGAAGTGGCTGAAACAACCGACTTTTCACCACAAGCCTTTAAAAAACTGAAACGAGAACCCGCCCTTGAAAAAGCGGAGATTGATTATCTAAAAAAGTTAGTGGAATTgggg                                  >  1:270846/11‑239 (MQ=11)
                    tCCTCATTCTCAAGTCATTTATTGATTAAAAGCCTTTCGTGAAAGAGGACTTGATGGCGTAAAATCACCTTATACTAAACCTAAAAtacccaaaatgaaaaaaatgcagattgaagtgactgaaacaaccgacttttcaccacaagccttgaaaaaactcaaacgagaagtcataggcgataatttcgccattgaagagatctttaattggtgataaatagagcttgccttcatcgcatttaaattcaat                      >  1:381160‑M1/1‑86 (MQ=255)
                     ccTCATTCTCAAGTCATTTATTGATTAAAAGCCTTTCGTGAAAGAGGACTTGATGGCGTAAAATCACCTTATACTAAACCTAAAAtacccaaaatgaaaaaaatgcagattgaagtgactgaaacaaccgacttttcaccacaagccttgaaaaaactcaaacgagaagtcataggcgataatttcgccattgaagagatctttaattggtgatagatagagcttgccttcatcgcatttaaattcaatca                    >  1:377559‑M1/1‑85 (MQ=255)
                     ccTCATTCTCAAGTCATTTATTGATTAAAAGCCTTTCGTGAAAGAGGACTTGATGGCGTAAAATCACCTTATACTAAACCTAAAAtacccaaaatgaaaaaaatgcagattgaagtgactgaaacaaccgacttttcaccacaagccttgaaaaaactcaaacgagaagtcataggcgataatttcgccattgaagag                                                                         <  1:263412‑M1/198‑114 (MQ=255)
                        cgttCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTGATGGCGTAAAATCACCTTACAC‑AAACCTAAAACCTCCAAAATTCTGAAACCCAAAATGAAAAAAACGCAGATTGAAGTGGCTGAAACAACCGACTTTTCACCACAAGCCTTTAAAAAACTGAAACGAGAACCCGCCCTTGAAAAAGCGGAGATTGATTATCTAAAAAAGTTAGTGGAATTGGGGCGATCAAAAACAAATTAAAc              <  1:349475/249‑1 (MQ=14)
                           tctcAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTGATGGCGTAAAATCACCTTACAC‑AAACCTAAAACCTCCAAAATTCTGAAACCCAAAATGAAAAAAACGCAGATTGAAGTGGCTGAAACAACCGACTTTTCACCACAAGCCTTTAAAAAACTGAAACGAGAACCCGCCCTTGAAAAAGCGGAGATTGATTATCTAAAAAAGTTAGTGGAATTGGGGCGATCAAAAACAAATTAAACTCg           >  1:687968/1‑251 (MQ=255)
                            ctcAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTGATGGCGTAAAATCACCTTACAC‑AAACCTAAAACCTCCAAAATTCTGAAACCCAAAATGAAAAAAACGCAGATTGAAGTGGCTGAAACAACCGACTTTTCACCACAAGCCTTTAAAAAACTGAAACGAGAACCCGCCCTTGAAAAAGCGGAGATTGATTATCTAAAAAAGTTAGTGGAATTGGGGCGATCAAAAAc                       >  1:736955/1‑238 (MQ=14)
                            ctcAAGTCATTTATTGATTAAAAGCCTTTCGTGAAAGAGGACTTGATGGCGTAAAATCACCTTATACTAAACCTAAAAtacccaaaatgaaaaaaatgcagattgaagtgactgaaacaaccgacttttcaccacaagccttgaaaaaactcaaacgagaagtcataggcgataatttcgccattgaagagatctttaattggtgataaatag                                                   <  1:508647‑M1/213‑136 (MQ=255)
                             tcAAGTCATTTATTGATTAAAAGCCTTTCGTGAAAGAGGACTTGATGGCGTAAAATCACCTTATACTAAACCTAAAAtacccaaaatgaaaaaaatgcagattgaagtgactgaaacaaccgacttttcaccacaagccttgaaaaaactcaaacgagaagtcataggcgataatttcgccattgaagagatctttaattggtgataaatagagcttgccttcatcgcatttaaattcaatcacatcggtg            <  1:501201‑M1/251‑175 (MQ=255)
                              cAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTGATGGCGTAAAATCACCTTACAC‑AAACCTAAAACCTCCAAAATTCTGAAACCCAAAATGAAAAAAACGCAGATTGAAGTGGCTGAAACAACCGACTTTTCACCACAAGCCTTTAAAAAACTGAAACGAGAACCCGCCCTTGAAAAAGCGGAGATTGATTATCTAAAAAAGTTAGTGGAATTgggg                                  <  1:647501/225‑1 (MQ=11)
                                 gTCATTTATTGATTAAAAGCCTTTCGTGAAAGAGGACTTGATGGCGTAAAATCACCTTATACTAAACCTAAAAtacccaaaatgaaaaaaatgcagattgaagtgactgaaacaaccgacttttcaccacaagccttgaaaaaactcaaacgagaagtcataggcgataatttcgccattgaagagatctttaattggtgataaatag                                                   <  1:714722‑M1/208‑136 (MQ=255)
                                   cATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTGATGGCGTAAAATCACCTTACAC‑AAACCTAAAACCTCCAAAATTCTGAAACCCAAAATGAAAAAAACGCAGATTGAAGTGGCTGAAACAACCGACTTTTCACCACAAGCCTTTAAAAAACTGAAACGAGAACCCGCCCTTGAAAAAGCGGAGATTGATTATCTAAAAAAGTTAGTGGAAt                                       <  1:621030/215‑1 (MQ=11)
                                   cATTTATTGATTAAAAGCCTTTCGTGAAAGAGGACTTGATGGCGTAAAATCACCTTATACTAAACCTAAAAtacccaaaatgaaaaaaatgcagattgaagtgactgaaacaaccgacttttcaccacaagccttgaaaaaactcaaacgagaagtcataggcgataatttcgccattgaagagatctttaattggtgataaatagagcttgccttcatcgcatttaaattcaatcacatcggtgacccat      >  1:131326‑M1/1‑71 (MQ=255)
                                   cATTTATTGATTAAAAGCCTTTCGTGAAAGAGGACTTGATGGCGTAAAATCACCTTATACTAAACCTAAAAtacccaaaatgaaaaaaatgcagattgaagtgactgaaacaaccgacttttcaccacaagccttgaaaaaactcaaacgagaagtcataggcgataatttcgccattgaagagatctttaattggtgataaatagagcttgccttcatcgcatttaaattcaatcacatcggtgacccat      >  1:26145‑M1/1‑71 (MQ=255)
                                    aTTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTGATGGCGTAAAATCACCTTACAC‑AAACCTAAAACCTCCAAAATTCTGAAACCCAAAATGAAAAAAACGCAGATTGAAGTGGCTGAAACAACCGACTTTTCACCACAAGCCTTTAAAAAACTGAAACGAGAACCCGCCCTTGAAAAAGCGGAGATTGATTATCTAAAAAAGTTAGTGGAATTGGGGCGATCAAAAACAAATTAAACTCGATTTAAAAg  <  1:472045/251‑1 (MQ=255)
                                    aTTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTGATGGCGTAAAATCACCTTACAC‑AAACCTAAAACCTCCAAAATTCTGAAACCCAAAATGAAAAAAACGCAGATTGAAGTGGCTGAAACAACCGACTTTTCACCACAAGCCTTTAAAAAACTGAAACGAGAACCCGCCCTTGAAAAAGCGGAGATTGATTATCTAAAAAAGTTAGTGGAATTGGGGCGATCAAAAACAAATTAAACTCGATTTAAAAg  <  1:240872/251‑1 (MQ=255)
                                    aTTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTGATGGCGTAAAATCACCTTACAC‑AAACCTAAAACCTCCAAAATTCTGAAACCCAAAATGAAAAAAACGCAGATTGAAGTGGCTGAAACAACCGACTTTTCACCACAAGCCTTTAAAAAACTGAAACGAGAACCCGCCCTTGAAAAAGCGGAGATTGATTATCTAAAAAAGTTAGTGGAATTGGGGCGATCAAAAACAAATTAAACTCGATTTAAAAg  <  1:395212/251‑1 (MQ=255)
                                    aTTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTGATGGCGTAAAATCACCTTACAC‑AAACCTAAAACCTCCAAAATTCTGAAACCCAAAATGAAAAAAACGCAGATTGAAGTGGCTGAAACAACCGACTTTTCACCACAAGCCTTTAAAAAACTGAAACGAAAACCCGCCCTTGAAAAAGCGGAGATTGATTATCTAAAAAAGTTAGTGGAATTGGGGCGATCAAAAACAAATTAAACTCGATTTAAAAg  <  1:155611/251‑1 (MQ=17)
                                            |                                                                                                                                                                                                                                                     
GTGTCTAAACTTTATAAGATTCCTCATTCTCAAGTCATTTATTGGTTAAAAGCCTTTCGTGAAAGAGGGCTTGATGGCATAAAATCGCCTTATACTAAACCTAAAACCTCCAAAATTCTGAAACCCAAAATGAAAAAAACGCAGATTGAAGTGGCTGAAACAACCGACTTTTCACCACAAGCCTTTAAAAAACTGAAACGAGAACCCGCCCTTGAAAAAGCGGAGATTGATTATCTAAAAAAGTTAGTGGAATTGGGGCGATCAAAAACAAATTAAACTCGATTTAAAAG  >  1/351448‑351737

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑