Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 226,027 C→T noncoding (269/2904 nt) rrlH → 23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913226,0270CT92.9% 40.3 / ‑4.4 28noncoding (269/2904 nt)rrlH23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/14);  minor base G (0/2);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

GCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGG  >  NC_000913/225826‑226232
                                                                                                                                                                                                         |                                                                                                                                                                                                             
gCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGgcgcgc                                                                                                                                                              <  1:144762/251‑1 (MQ=14)
 cGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGCAAACCCAGTGTGTTTAGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATAAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGCAAGgcgcg                                                                                                                                                               >  1:57640/1‑249 (MQ=9)
                       atGAACCGATATGACCGTGGATTTCCGAATGGGGAAACCCGGTGTGTTTCGACACACTATGATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTGCCCCGAAGAAAAGAAAGCAACCGAGGTTCCCCCAGTAGCTGCGAGCGAACGGGGAGGAGCCGAGAGCCTGAATCAGTGTGTGGGTTAGTGGGATCGTCTGGAAAGGAGCGCGATACAGGGTGAGAGCCCCGTa                                                                                                                                        <  1:167526/250‑1 (MQ=11)
                                     ccGGCGATTTCCGAATGGGGAAACCCAGTTTGTTTCGACACACTATCATTAACTGAATCCAGAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCATTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTGAAtn                                                                                                                                                                                                   <  1:49779/177‑2 (MQ=255)
                                              tCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGAGTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGTCGAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAggg                                                                                                                                                     <  1:219510/214‑1 (MQ=9)
                                                ctaaTGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTc                                                                                                              >  1:408864/3‑251 (MQ=14)
                                                                    gtTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCtg                                                                                                                     <  1:290833/224‑1 (MQ=14)
                                                                                 aTCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTc                                                                              <  1:192616/250‑1 (MQ=14)
                                                                                         cTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAgg                                                                                                                                                      <  1:350734/170‑1 (MQ=17)
                                                                                            aaTCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTGAATCAgtgtgtgtg                                                                                                                                                                                         <  1:413133/132‑1 (MQ=18)
                                                                                                  nnggNNNNNGAGNCGAACCGGGGGAACTGAAACAGCTAAGTAACCCGAGGAAAAGAAATCATCCGAGATTCCTCCAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTGAAGCAGTGTGTGTGGTAGGGGAAGCGTATGGaaa                                                                                                                                                                     <  2:283034/144‑1 (MQ=11)
                                                                                                        aaTGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAACTCAACCGAGATTCCCACAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTGAATCAGTGTGTGTATTAGTGGAAGCGTCTGGAAAgg                                                                                                                                                                   >  1:403871/1‑142 (MQ=11)
                                                                                                             gctggaCCGGNGGAACTGAAANAGCNNAGGANCNCNAGGAAAAGAAAGCAACCGAGATTCCCCCANTAGCGCCGAGCGAACGGGGAGGANNCTNNNGCNNGAATc                                                                                                                                                                                                   <  2:339057/100‑1 (MQ=2)
                                                                                                             gctggaCCGGGGGAACTGAAACCTCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTgcctc                                                                                                                                                                                                   >  1:339057/6‑101 (MQ=11)
                                                                                                                gAACCGGGGGAACTGAAACATCTAAGTACAACGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCgc                                                                                                                                         >  1:354129/1‑159 (MQ=255)
                                                                                                                                aaCATCTAAGGACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCTAGCGAACGGGGAGGAGCCTAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGGATCCTGTCTGAATATGGGGGGACCATCCTCAAAGGc                                                  <  1:147591/231‑1 (MQ=255)
                                                                                                                                  cATCTAAGTACCCCGAGGAAAATAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGTGGAGGAGCCTAGAGCCTGAATCAGTTTGGGTGGTAGTGGAAGCGGCTTGAAAGGCGCGCGATACACGGTGACGGCCCCGTACACAAAAATGCACATGCTGTGa                                                                                                                  <  1:382539/165‑1 (MQ=11)
                                                                                                                                            ncccGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAgg                                                                                                                                                                   >  1:88099/2‑106 (MQ=18)
                                                                                                                                             nccGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGcccc                                                                                                                                           >  1:278546/2‑129 (MQ=18)
                                                                                                                                             cccGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGcccc                                                                                                                                           <  2:278546/129‑1 (MQ=18)
                                                                                                                                                    naaaGAAATCAACCGAGCTTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTg           >  1:109672/2‑250 (MQ=14)
                                                                                                                                                     aaagaaagCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCt                                                                                                               <  1:401518/149‑1 (MQ=12)
                                                                                                                                                     aaaGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATg                                                                                                                        <  1:94803/140‑1 (MQ=18)
                                                                                                                                                     aaaGAAATCAACCGAGATTCCCCCAGNAGCGGCNNGCNAACGGGGAGGAGCCTAGAGNCTNNNTNNNTNTGTGTNTTNNTGGAAGCGTNTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATg                                                                                                                        >  2:94803/1‑140 (MQ=11)
                                                                                                                                                                       tCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCTAGCGCCTGAATCAGAGTGTGTGTTAGTGGAAGCGTCTGGAAAgga                                                                                                                                                                  >  1:402971/1‑79 (MQ=11)
                                                                                                                                                                                         cgaACGGGGAGGAGCCTAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACCTGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTg                                                                            >  1:400714/1‑148 (MQ=12)
                                                                                                                                                                                           aaCGGGGAGGAGCCTAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACATGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTATGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGTACCATCCTCCAAGGCTAAATGCTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGa      <  1:329751/216‑1 (MQ=255)
                                                                                                                                                                                                ggaggaGCCTAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAACAATGCACATGCTGTGAGCTCGCTGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAgg  >  1:224542/1‑215 (MQ=9)
                                                                                                                                                                                                         |                                                                                                                                                                                                             
GCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGG  >  NC_000913/225826‑226232

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