Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,644,169 G→T P39T (CCC→ACC)  rlmN ← dual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,644,1690GT93.1% 450.0 / 15.5 147P39T (CCC→ACC) rlmNdual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (6/4);  new base T (63/71);  total (71/76)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

AGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGT  >  NC_000913/2644072‑2644263
                                                                                                 |                                                                                              
aGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGttt                                                                                              >  2:2543265/1‑100 (MQ=255)
  tttGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTtttt                                                                                            <  2:433659/100‑1 (MQ=255)
   ttGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGtttttt                                                                                           >  1:2720003/1‑100 (MQ=255)
   ttGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGtttttt                                                                                           <  2:2133507/100‑1 (MQ=255)
    tGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTc                                                                                          <  2:336224/100‑1 (MQ=255)
    tGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTc                                                                                          <  2:2513084/100‑1 (MQ=255)
     gCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCa                                                                                         <  1:794876/100‑1 (MQ=255)
     gCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCa                                                                                         <  1:1911262/100‑1 (MQ=255)
     gCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCa                                                                                         <  1:1986011/100‑1 (MQ=255)
      ccGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCAc                                                                                        <  1:986129/100‑1 (MQ=255)
      ccGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCAc                                                                                        >  2:2689560/1‑100 (MQ=255)
         cgcAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCAcc                                                                                       <  1:2968662/98‑1 (MQ=255)
           cAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATTGTGATACATCCACTTCATCTCCTGATCGGCGCGGAAGGTTTTTTCACCTaaa                                                                                   <  2:2437358/100‑1 (MQ=255)
               acTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCtt                                                                               >  1:3019813/1‑100 (MQ=255)
                cTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACACCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCttt                                                                              >  1:1955640/1‑100 (MQ=255)
                 tttGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTa                                                                             <  2:2439827/100‑1 (MQ=255)
                 tttGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTa                                                                             <  1:239992/100‑1 (MQ=255)
                  ttgttgATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTaa                                                                            <  1:970331/100‑1 (MQ=255)
                   tgttgATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTaaa                                                                           >  1:277741/1‑100 (MQ=255)
                     ttgATGTCTGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTaaaaa                                                                         <  2:1857469/100‑1 (MQ=255)
                       gATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAAc                                                                       >  2:648048/1‑100 (MQ=255)
                       gATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAAc                                                                       <  1:3300702/100‑1 (MQ=255)
                       gATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAAc                                                                       >  2:88968/1‑100 (MQ=255)
                         tGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTc                                                                     >  2:2543362/1‑100 (MQ=255)
                           tCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGTGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTcac                                                                   >  1:641152/1‑98 (MQ=255)
                           tCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTccc                                                                   >  2:1268077/1‑100 (MQ=255)
                           tCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTccc                                                                   >  2:2721554/1‑100 (MQ=255)
                            cGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCg                                                                  >  1:377584/1‑100 (MQ=255)
                            cGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCg                                                                  <  2:1086454/100‑1 (MQ=255)
                            cGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCg                                                                  <  2:222978/100‑1 (MQ=255)
                             ggTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGc                                                                 >  2:1904990/1‑100 (MQ=255)
                             ggTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGc                                                                 >  2:2659536/1‑100 (MQ=255)
                             ggTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGc                                                                 <  1:2256839/100‑1 (MQ=255)
                             ggTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGc                                                                 >  1:1121516/1‑100 (MQ=255)
                              gTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCa                                                                <  1:1159250/100‑1 (MQ=255)
                                 atctcatcAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATct                                                             >  2:2017916/1‑100 (MQ=255)
                                  tctcatcAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATctg                                                            >  1:2578872/1‑100 (MQ=255)
                                  tctcatcAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATctg                                                            >  1:1042754/1‑100 (MQ=255)
                                   ntcatcaAAGTTGTCGCAGCAATAGTCATACATCCCCTTCATCACCAGATCCCCACGGACGGCTTTTCCACCTAAATCCTTAACCAACTCCCGCATctcc                                                           >  1:480693/2‑98 (MQ=255)
                                   ctcatcAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATctgc                                                           >  1:881655/1‑100 (MQ=255)
                                   ctcatcAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATctgc                                                           >  1:1563037/1‑100 (MQ=255)
                                     catcaAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATctgctg                                                         >  2:986203/1‑100 (MQ=255)
                                     catcaAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATctgctg                                                         >  2:776299/1‑100 (MQ=255)
                                     catcaAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATctgctg                                                         <  2:3242546/100‑1 (MQ=255)
                                      atcaAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGa                                                        >  1:1444080/1‑100 (MQ=255)
                                       tcaAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGAc                                                       >  1:2171077/1‑100 (MQ=255)
                                       tcaAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGAc                                                       >  1:871815/1‑100 (MQ=255)
                                        caAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACg                                                      >  1:1258224/1‑100 (MQ=255)
                                         aaaGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACgg                                                     >  2:1481475/1‑100 (MQ=255)
                                           agctctCGCAGCAATAGTGTTACATCCACTTCATCACCTGATCGGCGCCGAAGGTTTTTTCACCTAATTCTTTAAAAAACTCCCGCATCTGCTGACgg                                                     <  1:2750776/93‑1 (MQ=255)
                                              tGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTgag                                                <  1:684349/100‑1 (MQ=255)
                                              tGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTgag                                                <  1:986191/100‑1 (MQ=255)
                                               gTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTgaga                                               >  1:2139337/1‑100 (MQ=255)
                                                tCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTgag                                                <  1:562632/98‑1 (MQ=255)
                                                tCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGAt                                              >  1:1991550/1‑100 (MQ=255)
                                                  gcagcaATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATcc                                            <  1:776288/100‑1 (MQ=255)
                                                  gcagcaATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATcc                                            >  1:388153/1‑100 (MQ=255)
                                                     gcaATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCcagc                                         <  1:2789690/100‑1 (MQ=255)
                                                      caATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCcagca                                        >  1:1829252/1‑100 (MQ=255)
                                                      caATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCcagca                                        >  1:3087232/1‑100 (MQ=255)
                                                      caAGAGTCATACATCCCCTTCATCACCTGATCGGCGCGGAAGGCTTTTTCACCTAAATCTGTCAAAAACTCCCGCATCTGCTGACGGTTGAGATCcagca                                        <  1:9345/100‑1 (MQ=255)
                                                         tAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAgg                                      >  1:2363132/1‑99 (MQ=255)
                                                           gTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTg                                   <  1:2086188/100‑1 (MQ=255)
                                                           gTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTg                                   >  2:538960/1‑100 (MQ=255)
                                                             gATACATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGAt                                 >  1:3233330/1‑100 (MQ=255)
                                                             gATACATCCACTTCATCACCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGAt                                 >  1:2525929/1‑100 (MQ=255)
                                                                aCATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGAtttt                              <  2:2795301/100‑1 (MQ=255)
                                                                aCATCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGAtttt                              <  2:305278/100‑1 (MQ=255)
                                                                  aTCCACTTCATCACCTGATCGGCGCTGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTc                            <  1:3462140/100‑1 (MQ=255)
                                                                   tCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTcc                           <  2:448869/100‑1 (MQ=255)
                                                                   tCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTcc                           <  2:2910601/100‑1 (MQ=255)
                                                                   tCCACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTcc                           <  2:2872701/100‑1 (MQ=255)
                                                                    ccACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCa                          >  1:3441948/1‑100 (MQ=255)
                                                                    ccACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCa                          <  2:1829277/100‑1 (MQ=255)
                                                                     cACTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCAt                         <  2:2574899/100‑1 (MQ=255)
                                                                      aCTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATc                        >  1:1080219/1‑100 (MQ=255)
                                                                       cTTCATCACCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCt                       <  1:1477406/100‑1 (MQ=255)
                                                                          catcaCCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATGTGCTGACGGTTGAGATCCACCAGGTTTATTTTTCCATCttt                     >  1:1374371/1‑99 (MQ=255)
                                                                          catcaCCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTc                    >  2:887699/1‑100 (MQ=255)
                                                                          catcaCCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTc                    <  1:2420837/100‑1 (MQ=255)
                                                                           atcaCCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCg                   >  1:1847083/1‑100 (MQ=255)
                                                                            tcaCCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCgt                  <  2:3352443/100‑1 (MQ=255)
                                                                             caCCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCgtg                 >  1:153725/1‑100 (MQ=255)
                                                                             caCCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCgtg                 <  2:377589/100‑1 (MQ=255)
                                                                             caCCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCgtg                 <  2:138354/100‑1 (MQ=255)
                                                                              aCCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCgtgg                <  1:1290247/100‑1 (MQ=255)
                                                                              aCCTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCgtgg                >  2:1068989/1‑100 (MQ=255)
                                                                              aCCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCgtgg                <  1:138120/100‑1 (MQ=255)
                                                                               ccTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCgtggt               <  2:984404/100‑1 (MQ=255)
                                                                               ccTGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCgtgg                >  2:2648810/1‑99 (MQ=255)
                                                                                 tGATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGa             <  1:2543231/100‑1 (MQ=255)
                                                                                  gATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGAc            <  2:2510308/100‑1 (MQ=255)
                                                                                  gATCGGCGCGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGAc            >  1:452913/1‑100 (MQ=255)
                                                                                      ngcgcGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGttt        >  2:3151483/2‑100 (MQ=255)
                                                                                       gcgcGGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGtttt       >  1:693699/1‑100 (MQ=255)
                                                                                         gcgGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCa     >  2:3001558/1‑100 (MQ=255)
                                                                                          cgGAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAg    <  1:2647323/100‑1 (MQ=255)
                                                                                           ggAAGGTTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAgg   <  2:1444102/100‑1 (MQ=255)
                                                                                            gAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGgt  >  1:209422/1‑100 (MQ=255)
                                                                                             aaGGTTTTTTCTCCTAAATCTTTAAAACACTCCCGCATCTGCTTACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAg    <  1:2724178/97‑1 (MQ=255)
                                                                                                 |                                                                                              
AGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGT  >  NC_000913/2644072‑2644263

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