Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,513,846 C→T 5.8% S5S (AGC→AGT rsmH → ribosomal RNA small subunit methyltransferase H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,513,8460CT 5.8% 208.1 / 2.4 69S5S (AGC→AGTrsmHribosomal RNA small subunit methyltransferase H
Reads supporting (aligned to +/- strand):  ref base C (41/24);  new base T (2/2);  total (43/26)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

TGCTGGTGGGGCAACTGAACAAATTCCAGCTGTGGGATGAAGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGG  >  NC_002947/1513705‑1513994
                                                                                                                                              |                                                                                                                                                     
tGCTGGTGGGGCAACTGAACAAATTCCAGCTGTGGGATGAAGATGCCTGGAACGCGGTTTCGGCAGCTGACATTGAAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTc                                                                                                                                                 >  7:243794/1‑148 (MQ=255)
   tggtggGGCAACTGAACAAATTCCAGCTGTGGGATGAAGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAAcc                                                                                                                                             <  4:324373/149‑1 (MQ=255)
   tggtggGGCAACTGAACAAATTCCAGCTGTGGGATGAAGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAAc                                                                                                                                              >  6:84726/1‑148 (MQ=255)
     gtggGGCAACTGAACAAATTCCAGCTGTGGGATGAAGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACca                                                                                                                                            <  2:40786/148‑1 (MQ=255)
          gCAACTGAACAAATTCCAGCTGTGGGATGAAGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACcac                                                                                                                                           >  1:287513/1‑144 (MQ=255)
          gCAACTGAACAAATTCCAGCTGTGGGATGAAGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACcac                                                                                                                                           <  2:287513/144‑1 (MQ=255)
          gCAACTGAACAAATTCCAGCTGTGGGATGAAGATGCCTGGAACGCGGTATCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTa                                                                                                                                       >  8:364386/1‑148 (MQ=255)
            aaCTGAACAAATTCCAGCTGTGGGATGAAGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCg                                                                                                                                    >  5:188662/1‑149 (MQ=255)
            aaCTGAACAAATTCCAGCTGTGGGATGAAGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCg                                                                                                                                    <  7:271400/149‑1 (MQ=255)
                 aaaaaaTTACAGCTGTGGGATGAAGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGTCGAAATGCACGACGAATTGCGTGACATCATACTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTg                                                                                                                               >  8:27164/4‑149 (MQ=255)
                 aaCAAATTCCAGCTGTGGGATGAAGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTg                                                                                                                               >  8:174525/1‑149 (MQ=255)
                  aCAAATTCCAGCTGTGGGATGAAGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGTGGCTTCAACCACATTCCCGTCCTGc                                                                                                                              >  6:20975/1‑149 (MQ=255)
                  aCAAATTCCAGCTGTGGGATGAAGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGTGGCTTCAACCACATTACCGTCCTGc                                                                                                                              >  6:15670/1‑149 (MQ=255)
                            gCTGTGGGATGAAGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGc                                                                                                                    >  7:47137/1‑149 (MQ=255)
                               gtcggaggaAGATGCCTGGCACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGt                                                                                                                 <  8:243794/146‑1 (MQ=255)
                                    aTGAAGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAgg                                                                                                            >  2:180048/1‑149 (MQ=255)
                                    aTGAAGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAgg                                                                                                            >  5:164453/1‑149 (MQ=255)
                                       aaGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGTGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCAt                                                                                                         <  5:20975/149‑1 (MQ=255)
                                       aaGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCAt                                                                                                         >  6:162984/1‑149 (MQ=255)
                                        aGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGTGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCAt                                                                                                         <  5:15670/148‑1 (MQ=255)
                                             ccTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCt                                                                                                   >  7:195026/1‑149 (MQ=255)
                                              cTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCtt                                                                                                  >  4:81032/1‑149 (MQ=255)
                                              cTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCtt                                                                                                  <  7:27164/149‑1 (MQ=255)
                                               tGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTg                                                                                                 >  8:318808/1‑149 (MQ=255)
                                                ggAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTgc                                                                                                <  2:136235/149‑1 (MQ=255)
                                                  aaCGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTgcgc                                                                                             >  2:66804/1‑149 (MQ=255)
                                                   aCGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTgcgcg                                                                                            >  4:11105/1‑149 (MQ=255)
                                                             ggCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTg                                                                                  >  2:264795/1‑149 (MQ=255)
                                                             ggCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTg                                                                                  >  5:238031/1‑149 (MQ=255)
                                                             ggCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTg                                                                                  >  3:30002/1‑149 (MQ=255)
                                                                aGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTa                                                                               >  4:320342/1‑149 (MQ=255)
                                                                 gCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTAt                                                                              >  4:95143/1‑149 (MQ=255)
                                                                       cTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACTCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGa                                                                         >  1:91623/1‑149 (MQ=255)
                                                                       cTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGAc                                                                        >  6:381470/1‑149 (MQ=255)
                                                                          gCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACgg                                                                      <  2:38483/148‑1 (MQ=255)
                                                                            aGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCAc                                                                   >  3:42554/1‑149 (MQ=255)
                                                                               tATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACttt                                                                >  7:360306/1‑149 (MQ=255)
                                                                                 ttCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTc                                                               <  1:236568/148‑1 (MQ=255)
                                                                                   caacaaCCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCgg                                                             <  7:364386/148‑1 (MQ=255)
                                                                                    aacaacCGGGCGCCCTGCCCGACGAATTGCGTGACGTGATCCTGTCACCATAGCTAGCGGCTTCAACCCCATTACCGTCCTGCTCGACGCAGCTGTCGAGGCATGGGCCTTGCGCGCCGCCGGTTGCGATCTGGACGGCACTTTCGTcc                                                           >  5:373992/1‑149 (MQ=255)
                                                                                    aacaacCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGcc                                                           >  5:348247/1‑149 (MQ=255)
                                                                                      caacCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCgg                                                         <  3:81032/149‑1 (MQ=255)
                                                                                       aacCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCgg                                                         <  1:264795/149‑1 (MQ=255)
                                                                                       aacCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCggg                                                        >  1:321926/1‑149 (MQ=255)
                                                                                                  aTGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGc                                             <  3:11105/149‑1 (MQ=255)
                                                                                                   tGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGcc                                            <  8:47137/149‑1 (MQ=255)
                                                                                                    gCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAg                                              <  3:292005/146‑1 (MQ=255)
                                                                                                    gCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAg                                              >  4:292005/1‑146 (MQ=255)
                                                                                                      ccGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTc                                         >  5:80790/1‑149 (MQ=255)
                                                                                                      ccGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCCGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTc                                         <  1:66804/149‑1 (MQ=255)
                                                                                                       cgacgaATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCt                                        >  6:197838/1‑149 (MQ=255)
                                                                                                        gacgaATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGACGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCt                                        >  4:325094/1‑148 (MQ=255)
                                                                                                         acgaCTTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGa                                      <  4:223586/149‑1 (MQ=255)
                                                                                                           gaATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATa                                    <  6:188662/149‑1 (MQ=255)
                                                                                                           gaATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATa                                    <  3:325094/149‑1 (MQ=255)
                                                                                                                gCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAg                               >  6:375499/1‑149 (MQ=255)
                                                                                                                  gTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGc                                                      <  6:355430/124‑1 (MQ=255)
                                                                                                                  gTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGc                                                      >  5:355430/1‑124 (MQ=255)
                                                                                                                  gTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGCGGCGGa                                                   >  3:47630/1‑127 (MQ=255)
                                                                                                                  gTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAACCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGa                                                   <  4:47630/127‑1 (MQ=255)
                                                                                                                     aCCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGc                          >  6:381504/1‑149 (MQ=255)
                                                                                                                           aTCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCggg                     <  4:42554/149‑1 (MQ=255)
                                                                                                                               tgtgACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCg                  >  2:305857/1‑148 (MQ=255)
                                                                                                                                   aCCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAg              <  5:197838/148‑1 (MQ=255)
                                                                                                                                    ccATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAggg            >  7:39074/1‑149 (MQ=255)
                                                                                                                                    ccATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAggg            >  1:139582/1‑149 (MQ=255)
                                                                                                                                     cATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGggc           >  3:71632/1‑149 (MQ=255)
                                                                                                                                     cATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGggc           <  5:381470/149‑1 (MQ=255)
                                                                                                                                        agatagCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGggcggc        >  4:310434/1‑149 (MQ=255)
                                                                                                                                              cGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTgg  >  7:246406/1‑149 (MQ=255)
                                                                                                                                              cGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTgg  >  4:183320/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
TGCTGGTGGGGCAACTGAACAAATTCCAGCTGTGGGATGAAGATGCCTGGAACGCGGTTTCGGCAGCTGACCTTGCAGCTATTCAACAACCGGGCGCCATGCCCGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGG  >  NC_002947/1513705‑1513994

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