Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,604,784 G→T R131L (CGC→CTC)  rsmD → 16S rRNA m(2)G966 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,604,7840GT98.4% 213.3 / ‑3.3 61R131L (CGC→CTC) rsmD16S rRNA m(2)G966 methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (29/31);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CGCGGTTTCTCAGCAGTTAATTAAGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGA  >  NC_000913/3604647‑3604919
                                                                                                                                         |                                                                                                                                       
cgcgGTTTCTCAGCAGTTAATTAAGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTtcctcc                                                                                                                                       >  2:833428/1‑140 (MQ=255)
cgcgGTTTCTCAGCAGTTAATTAAGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTtcctcc                                                                                                                                       >  2:1106572/1‑140 (MQ=255)
cgcgGTTTCTCAGCAGTTAATTAAGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTtcctcc                                                                                                                                       >  2:280508/1‑140 (MQ=255)
    gTTTCTCAGCAGTTAATTAAGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTgg                                                                                                                                   <  2:1063588/140‑1 (MQ=255)
    gTTTCTCAGCAGTTAATTAAGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTgg                                                                                                                                   <  2:1094641/140‑1 (MQ=255)
    gTTTCTCAGCAGTTAATTAAGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTgg                                                                                                                                   <  2:235695/140‑1 (MQ=255)
    gTTTCTCAGCAGTTAATTAAGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTgg                                                                                                                                   <  1:433589/140‑1 (MQ=255)
    gTTTCTCAGCAGTTAATTAAGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTga                                                                                                                                   <  1:780185/140‑2 (MQ=255)
     tttCTCAGCAGTTAATTAAGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGc                                                                                                                                  >  2:326957/1‑140 (MQ=255)
     tttCTCAGCAGTTAATTAAGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGc                                                                                                                                  <  2:449137/140‑1 (MQ=255)
      ttCTCAGCAGTTAATTAAGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCt                                                                                                                                 <  1:95905/140‑1 (MQ=255)
      ttCTCAGCAGTTAATTAAGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCt                                                                                                                                 <  1:684890/140‑1 (MQ=255)
         tcAGCAGTTAATTAAGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGt                                                                                                                              <  1:75539/140‑1 (MQ=255)
               gTTAATTAAGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTagaag                                                                                                                        >  1:370755/1‑140 (MQ=255)
                      aaGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGAt                                                                                                                 <  1:1060312/140‑1 (MQ=255)
                       aGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTagaa                                                                                                                         >  2:1012104/1‑131 (MQ=255)
                       aGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTagaa                                                                                                                         <  1:1012104/131‑1 (MQ=255)
                          aTCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTTAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAAt                                                                                                             >  1:848376/1‑140 (MQ=255)
                          aTCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAAt                                                                                                             >  1:869326/1‑140 (MQ=255)
                           tCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAAtt                                                                                                            <  2:370755/140‑1 (MQ=255)
                            cTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAAttt                                                                                                           >  1:450336/1‑140 (MQ=255)
                             tGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTa                                                                                                          >  2:700649/1‑140 (MQ=255)
                              ggCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTAc                                                                                                         <  1:280508/140‑1 (MQ=255)
                                 gACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTgg                                                                                                      <  2:869326/140‑1 (MQ=255)
                                 gACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTgg                                                                                                      <  2:450336/140‑1 (MQ=255)
                                      tAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGAt                                                                                                 <  1:833428/140‑1 (MQ=255)
                                              ggCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTgg                                                                                                      >  2:424183/1‑127 (MQ=255)
                                              ggCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTgg                                                                                                      <  1:424183/127‑1 (MQ=255)
                                              ggCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGAt                                                                                                 <  1:182593/132‑1 (MQ=255)
                                              ggCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGAt                                                                                                 >  2:182593/1‑132 (MQ=255)
                                              ggCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACggctg                                                                                         >  1:771127/1‑140 (MQ=255)
                                                        gcgTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGAc                                                                               >  2:152313/1‑140 (MQ=255)
                                                         cgTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACg                                                                              >  1:207423/1‑140 (MQ=255)
                                                             gtgAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGc                                                                          <  1:102694/140‑1 (MQ=255)
                                                             gtgAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGc                                                                          >  1:662238/1‑140 (MQ=255)
                                                              tgAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGcc                                                                         <  1:1106572/140‑1 (MQ=255)
                                                                aaCAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCt                                                                       >  1:315132/1‑140 (MQ=255)
                                                                        cgcgATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATg                                                               >  2:901284/1‑140 (MQ=255)
                                                                          cgATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTc                                                             >  1:749585/1‑140 (MQ=255)
                                                                            aTGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGa                                                           >  2:484875/1‑140 (MQ=255)
                                                                               tCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAg                                                        >  1:643945/1‑140 (MQ=255)
                                                                                 aTTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCg                                                      <  1:853185/140‑1 (MQ=255)
                                                                                       ggCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCg                                                >  1:599319/1‑140 (MQ=255)
                                                                                            aaaaaGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAAc                                           <  2:207423/140‑1 (MQ=255)
                                                                                             aaaaGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACg                                          <  2:749585/140‑1 (MQ=255)
                                                                                                    acacCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGcc                                   <  2:848376/140‑1 (MQ=255)
                                                                                                      acCGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCa                                 <  1:152313/140‑1 (MQ=255)
                                                                                                        cGCATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACt                               >  2:500856/1‑140 (MQ=255)
                                                                                                          cATAATATCGTGTTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGt                             >  2:69048/1‑140 (MQ=255)
                                                                                                                    tgtTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCaaa                    <  2:319910/139‑1 (MQ=255)
                                                                                                                    tgtTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCaaa                    >  1:319910/1‑139 (MQ=255)
                                                                                                                    tgtTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAAc                   <  1:987783/140‑1 (MQ=255)
                                                                                                                     gtTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACt                  >  2:380172/1‑140 (MQ=255)
                                                                                                                     gtTTGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACt                  <  1:700649/140‑1 (MQ=255)
                                                                                                                        tGTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGt               >  1:621172/1‑140 (MQ=255)
                                                                                                                         gTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTc              >  1:497385/1‑140 (MQ=255)
                                                                                                                         gTCGATCCACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTc              <  1:500856/140‑1 (MQ=255)
                                                                                                                                cACCGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACAt       <  1:326957/140‑1 (MQ=255)
                                                                                                                                   cGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCa                      >  1:1107445/1‑122 (MQ=255)
                                                                                                                                   cGTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCa                      <  2:1107445/122‑1 (MQ=255)
                                                                                                                                    gTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCggg   >  1:398029/1‑140 (MQ=255)
                                                                                                                                    gTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCggg   >  2:167993/1‑140 (MQ=255)
                                                                                                                                    gTTCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCggg   >  2:855752/1‑140 (MQ=255)
                                                                                                                                     ttCCTCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGa  >  2:94727/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                       
CGCGGTTTCTCAGCAGTTAATTAAGAATCTGGCGACACTAAAAGCAGGCAATGCACGCGTGGTGAACAGCAACGCGATGTCATTCCTGGCGCAAAAAGGTACACCGCATAATATCGTGTTTGTCGATCCACCGTTCCGCCGTGGCTTGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGA  >  NC_000913/3604647‑3604919

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