Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,410,053 (T)9→8 100% intergenic (+47/‑80) rlmB → / → yjfI 23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,410,0450T.98.0% 685.0 / 3.1 153intergenic (+39/‑88)rlmB/yjfI23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base . (72/78);  total (74/79)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 179 total aligned reads displayed.

TGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATG  >  NC_000913/4409948‑4410140
                                                                                                 |                                                                                               
tGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAAATCTGGATGGctt                                                                                               <  2:53267/100‑1 (MQ=255)
    cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                           <  1:1004717/100‑1 (MQ=255)
    cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                           >  2:1567702/1‑100 (MQ=255)
    cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                           >  2:1833761/1‑100 (MQ=255)
     gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttttt                                                                                          >  2:2869135/1‑100 (MQ=255)
     gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAAATCTGGATGGcttttcat                                                                                          <  2:1680986/100‑4 (MQ=38)
     gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcat                                                                                          <  1:811227/100‑4 (MQ=38)
     gTTTCGGTTCCGACCTTAATTTGCTTATTTAATACGGTACGCCAGCGCAGCTAATTTCTCAGTTTTATGACAAGCCATCCAGATTTGTATCGtttttttt                                                                                          <  1:3062352/100‑1 (MQ=255)
      tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGTGCCAGCGCAGCTGATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                         <  1:2768746/100‑2 (MQ=255)
      tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                         <  1:2504849/100‑2 (MQ=255)
      tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                         <  2:475099/100‑2 (MQ=255)
      tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                         <  1:1253488/100‑2 (MQ=255)
      tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttttt                                                                                          >  1:907541/1‑99 (MQ=255)
      ttgcGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcata                                                                                         <  1:401803/97‑5 (MQ=37)
       ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgt                                                                                        <  1:1658133/100‑3 (MQ=255)
       ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgt                                                                                        >  2:1228456/1‑98 (MQ=255)
       ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                         <  1:482540/99‑2 (MQ=255)
       ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                         <  1:571323/99‑2 (MQ=255)
        tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCTGAATTATGAATAACCATCCAAATCTGGAAGGcttttcatagt                                                                                       <  2:2065612/100‑7 (MQ=25)
         cGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtct                                                                                      <  1:503853/100‑5 (MQ=255)
         cGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                       <  1:640597/99‑4 (MQ=255)
          ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAAATCTGGATGGcttttcataattc                                                                                     <  2:334244/100‑9 (MQ=37)
                 aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCACAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTTTAGCt                                                                               >  1:2463776/1‑98 (MQ=255)
                  ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGt                                                                            >  2:232695/1‑100 (MQ=255)
                   cGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTa                                                                           <  1:616011/100‑1 (MQ=255)
                   cGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTa                                                                           >  2:1173836/1‑100 (MQ=255)
                   cGGAATTTGCTTATTTGAAGCGGTGAGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTa                                                                           <  1:1513545/100‑1 (MQ=255)
                    ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTaa                                                                          >  1:706974/1‑100 (MQ=255)
                    ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTaa                                                                          <  2:146101/100‑1 (MQ=255)
                    ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTaa                                                                          <  1:3035975/100‑1 (MQ=255)
                    ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTaa                                                                          >  1:1081736/1‑100 (MQ=255)
                     gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAg                                                                         <  1:1800259/100‑1 (MQ=255)
                     gAATTTGCATATTAGAAGAGGTGAGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCAGCCAGATTTGGATGGTTTTTTTAGTCTATAGCTGGTaa                                                                          <  1:294860/99‑1 (MQ=255)
                        tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATa                                                                      <  2:72166/100‑1 (MQ=255)
                        tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATa                                                                      <  1:2497531/100‑1 (MQ=255)
                          tGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAAt                                                                    >  1:1318445/1‑100 (MQ=255)
                          tGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCACAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAAt                                                                    >  1:1633163/1‑100 (MQ=255)
                            cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                  <  2:1748675/100‑1 (MQ=255)
                            cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                  >  2:2499717/1‑100 (MQ=255)
                            cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                  >  2:2769988/1‑100 (MQ=255)
                              tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                                >  2:1409445/1‑100 (MQ=255)
                               aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                               >  2:1275764/1‑100 (MQ=255)
                                tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTa                                                              <  2:406538/100‑1 (MQ=255)
                                tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTa                                                              >  2:2248227/1‑100 (MQ=255)
                                 ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                             >  1:406798/1‑100 (MQ=255)
                                 ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                             >  2:704611/1‑100 (MQ=255)
                                   gAAGCGGTTCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTg                                                           <  1:2710252/100‑1 (MQ=255)
                                    aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                          <  1:368009/100‑1 (MQ=255)
                                    aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                          >  2:2966751/1‑100 (MQ=255)
                                    aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                          <  2:2117128/100‑1 (MQ=255)
                                     aGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCa                                                         <  1:2615752/100‑1 (MQ=255)
                                       cGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaaa                                                       <  1:1126372/100‑1 (MQ=255)
                                        ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAat                                                      >  1:18566/1‑100 (MQ=255)
                                              ccAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                                <  2:1448205/100‑1 (MQ=255)
                                               cAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                               >  2:926643/1‑100 (MQ=255)
                                                aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                              >  2:3001324/1‑100 (MQ=255)
                                                aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                              <  2:406802/100‑1 (MQ=255)
                                                 gcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTAGAGCTGGTAAGATAATTACGTATTGCAAATATTCAtttt                                             <  1:200814/100‑1 (MQ=255)
                                                  cgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTtta                                            >  2:434273/1‑100 (MQ=255)
                                                  cgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTtta                                            >  2:1639659/1‑100 (MQ=255)
                                                   gcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTATATAGCTGGTAACATAATTACGTATTGTAAATATTGATTttat                                           >  1:3048495/1‑100 (MQ=255)
                                                    cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                          <  1:2643132/100‑1 (MQ=255)
                                                    cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                          <  1:2821053/100‑1 (MQ=255)
                                                        tAATTTCTTATAATTATGAAAATCCATACAGATTTGGATGGTTTTTTCTGTCTATAGCTGGTAACATAATTACGTATTGCAAATATTCATTTTATTaa                                        <  1:2354652/98‑1 (MQ=255)
                                                          aTTTCTCATAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAtt                                    <  1:1726013/100‑1 (MQ=255)
                                                          aTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAtt                                    >  1:886569/1‑100 (MQ=255)
                                                              ctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTTATTTTATTAATATTTaaa                                >  1:522184/1‑100 (MQ=255)
                                                              ctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                <  1:2380190/100‑1 (MQ=255)
                                                              ctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                <  2:2924513/100‑1 (MQ=255)
                                                              ctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                >  2:2836786/1‑100 (MQ=255)
                                                               tcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                               <  2:886574/100‑1 (MQ=255)
                                                                  gAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTa                            >  2:1935718/1‑100 (MQ=255)
                                                                    gacaaaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGTTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAATTAAt                          >  2:101686/7‑100 (MQ=37)
                                                                     ttATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                         <  2:1081742/100‑1 (MQ=255)
                                                                      tATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                        >  1:2031701/1‑100 (MQ=255)
                                                                      caaaaaaaACCATCCAAATCTGGATGGCTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                          >  1:2473900/5‑98 (MQ=25)
                                                                      caaaaaaaACCATCCAAATCTGGATGGCTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                        >  2:2321886/5‑100 (MQ=25)
                                                                       caaaaaaACCATCCAAATCTGGATGGCTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                       >  2:274219/4‑100 (MQ=37)
                                                                       aTGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                       <  1:1549778/100‑1 (MQ=255)
                                                                         gAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                     <  2:17147/100‑1 (MQ=255)
                                                                          aaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat                    >  2:20215/1‑100 (MQ=255)
                                                                           aaaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                   <  1:1155550/100‑1 (MQ=255)
                                                                             aaccATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGa                 >  2:288711/3‑100 (MQ=255)
                                                                              gCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATggn                 <  2:792748/99‑2 (MQ=255)
                                                                              accATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaa                >  1:2631269/2‑100 (MQ=255)
                                                                               ccATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat               <  1:1261315/100‑1 (MQ=255)
                                                                                 aTCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataa             <  1:926940/100‑1 (MQ=255)
                                                                                   ccAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata           <  1:2345956/100‑1 (MQ=255)
                                                                                    cAGATTTGGATGGTTTTTTTTTTCTATAGCTGGTAAGATAATTACGTATTGTAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat          >  1:1769160/1‑100 (MQ=255)
                                                                                    cAGATTTGGATGGTTTTTTTT‑‑TCTATAGCTGGTAAGATAATTACGTGTTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatata         >  1:1695296/1‑100 (MQ=255)
                                                                                    cAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat          <  2:2638424/100‑1 (MQ=255)
                                                                                    cAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat          >  1:2882655/1‑100 (MQ=255)
                                                                                     aGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatata         <  2:1638295/100‑1 (MQ=255)
                                                                                       aTTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATg       <  2:2718465/100‑1 (MQ=255)
                                                                                       aTTTGGATGGTTTTTTTTGTCTATAGCTGCTAAGATAATTACGTATTGCAAATATTCATTTTATTATTATTTAAACTTAATCTATATGGAATAATATATg       <  1:1833750/100‑1 (MQ=255)
                                                                                        tttGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGa      <  1:2292413/100‑1 (MQ=255)
                                                                                         ttGGATTGTTTTTTTTGTTTTTAGTTGTTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata           <  1:2078363/94‑1 (MQ=255)
                                                                                          tGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACa    >  1:1802106/1‑100 (MQ=255)
                                                                                          tGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACa    <  1:2038580/100‑1 (MQ=255)
                                                                                            gATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATg  >  2:1685019/1‑100 (MQ=255)
                                                                                                 |                                                                                               
TGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATG  >  NC_000913/4409948‑4410140

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