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.90.7% 390.9 / 30.5 106intergenic (+39/‑88)rlmB/yjfI23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base T (7/1);  new base . (60/36);  total (68/38)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

CCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAT  >  NC_000913/4409946‑4410139
                                                                                                   |                                                                                              
ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGt                                                                                                >  2:574324/1‑100 (MQ=255)
 cTGAACGCTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGt                                                                                                >  1:488301/1‑99 (MQ=255)
   gAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttt                                                                                             <  2:67245/100‑1 (MQ=255)
     aCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttc                                                                                           <  2:208440/100‑2 (MQ=38)
      cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                          >  2:1618739/1‑100 (MQ=255)
      cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttca                                                                                          <  2:2192/100‑3 (MQ=38)
       gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcat                                                                                         <  1:190953/100‑4 (MQ=37)
       gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcat                                                                                         <  2:531872/100‑4 (MQ=38)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                        <  2:48488/100‑2 (MQ=255)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                      <  2:447613/100‑4 (MQ=255)
           cGGTTGCGACCGGAATTTGCTTATTTTAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtct                                                                                     <  1:2700614/100‑5 (MQ=255)
            ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtct                                                                                     >  1:963378/1‑95 (MQ=255)
            ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCta                                                                                   >  2:1437996/1‑100 (MQ=255)
            ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCta                                                                                   >  2:2310102/1‑100 (MQ=255)
            ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCta                                                                                   >  2:1796737/1‑100 (MQ=255)
            ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAAATCTGGATGGctttagatcggaa                                                                                    >  1:415909/1‑91 (MQ=255)
             gTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCTGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCtat                                                                                  <  1:2047916/100‑1 (MQ=255)
                gCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGc                                                                               <  2:2259306/100‑1 (MQ=255)
                  gACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTg                                                                             >  2:618725/1‑100 (MQ=255)
                  gACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTg                                                                             >  2:1437875/1‑100 (MQ=255)
                   aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTgg                                                                            <  2:128067/100‑1 (MQ=255)
                    ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGt                                                                           >  1:1321614/1‑100 (MQ=255)
                     cGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTa                                                                          >  2:900118/1‑100 (MQ=255)
                     cGGAATTTGCTTATTTGAAGCGGTGAGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGt                                                                           >  1:708466/1‑99 (MQ=255)
                      ggAATTTGCTTAATTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTaa                                                                         >  2:1919499/1‑100 (MQ=255)
                       gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAg                                                                        <  2:2133691/100‑1 (MQ=255)
                           ttGCTTATTTGAAGCGGTGCTCCAGCGCAGCTAATTTCTCAGAATTAGGAAAAGCCATCCATATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATaa                                                                    <  2:2432645/100‑1 (MQ=255)
                             gCTTATTTTGAGCGGTGGGCGAGCGCAGTTAATTTCTCAGAGTTAGGAAAAGGCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATa                                                                     <  1:2517599/97‑1 (MQ=255)
                                tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTTTAGCTGGTAAGATAATTACg                                                               >  1:2172063/1‑100 (MQ=255)
                                tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                               >  2:1217635/1‑100 (MQ=255)
                                 aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                              >  2:1883916/1‑100 (MQ=255)
                                 aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                              >  1:1434223/1‑100 (MQ=255)
                                  tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTa                                                             >  2:272601/1‑100 (MQ=255)
                                  gttGAAACGGTGAGCCAGCGGAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATGTGCATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACgga                                                             <  1:2755125/99‑3 (MQ=255)
                                    tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                           >  2:142512/1‑100 (MQ=255)
                                    tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                           <  2:2537585/100‑1 (MQ=255)
                                    tGAAGCGGTGCGCCAGCACAGCTAATTTCTCAGAATTATGAAAAGCCATACAGCTTTGGCTGGTTTTTTTTTTTTCTATCTGGTAATATAATTa                                                                 >  1:676784/1‑94 (MQ=255)
                                      aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                         >  2:1289910/1‑100 (MQ=255)
                                       aGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCa                                                        <  2:1177585/100‑1 (MQ=255)
                                         cGGTGCGCCAGCGCAGCTAAGTTCTCAGAATTATGAAAAGCCGTCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAATATAATTACGTATTGCaaa                                                      <  2:330489/100‑1 (MQ=255)
                                          ggTGCGCCAGCGCGGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTATATAGCTGGTAAGATAATTACGTATTGCAAat                                                     >  1:2426916/1‑100 (MQ=255)
                                          ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTTTATTGTTGGTAAGATAATTACGTATTTCAAat                                                     >  1:2485228/1‑100 (MQ=255)
                                          ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAat                                                     >  2:1708887/1‑100 (MQ=255)
                                          ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAat                                                     >  2:2305258/1‑100 (MQ=255)
                                             gcgcCAGCGCAGCTAACTAATCAGAATAATGAACAGCCAGCGAGATTTGCATGGTCTATTTTGTCGATAGCTGTTAAGATAATTACGTATTACaaacact                                                  >  1:2309218/1‑96 (MQ=255)
                                               gcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTc                                                 >  1:917332/1‑99 (MQ=255)
                                               gcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCa                                                >  2:2759802/1‑100 (MQ=255)
                                                ccAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                               <  2:667782/100‑1 (MQ=255)
                                                ccAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                               <  1:1596780/100‑1 (MQ=255)
                                                 cAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                              >  2:2852045/1‑100 (MQ=255)
                                                  aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                             >  1:341676/1‑100 (MQ=255)
                                                  aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                             <  1:2665957/100‑1 (MQ=255)
                                                  aGCGCAGCTAATGTGTCAGATTTATTAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                             <  2:537265/100‑1 (MQ=255)
                                                   gcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtttt                                            >  2:2327624/1‑100 (MQ=255)
                                                   gcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtttt                                            <  2:1380493/100‑1 (MQ=255)
                                                     gcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                          >  2:1825741/1‑100 (MQ=255)
                                                       aGCTAATTTCTCTGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                        <  1:1228060/100‑1 (MQ=255)
                                                       aGCTAATTTCTCAGAATTATGAAGAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                         <  1:2363732/99‑1 (MQ=255)
                                                        gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                        >  2:2419596/1‑99 (MQ=255)
                                                            aTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAttt                                  >  2:1948318/1‑100 (MQ=255)
                                                               tctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTTTTTTTCTGCTAAGATAATTCCGTATTGCAAATATTCATTTTATTAATATTTaa                                >  1:465221/1‑100 (MQ=255)
                                                               tctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaa                                <  1:1791133/100‑1 (MQ=255)
                                                                ctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                               <  2:1321619/100‑1 (MQ=255)
                                                                ctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                               >  1:2128472/1‑100 (MQ=255)
                                                                ctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                               >  2:2522640/1‑100 (MQ=255)
                                                                 tcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                              >  2:1645886/1‑100 (MQ=255)
                                                                 tcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTAATTTTATTAATATTTAAAc                              >  1:2666564/1‑100 (MQ=255)
                                                                  cAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                             >  1:976616/1‑100 (MQ=255)
                                                                   aGACTTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                              <  2:881109/98‑1 (MQ=255)
                                                                    gcatgaTGAATAGCCCTTCAGATAGGGATGGTTTTTTTTGTGTATGTCTGGAAAGATAATTACGTATTGCAAATATTCATTTTATTAATAGTTAAACtt                            <  1:770268/94‑1 (MQ=255)
                                                                      gacaaaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                         >  2:1938437/7‑100 (MQ=37)
                                                                        caaaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                       >  2:1568730/5‑100 (MQ=38)
                                                                        caaaaaaaACCATCCAAATCTGGATGGCTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                       >  2:873139/5‑100 (MQ=25)
                                                                         attaaaaACCATCCAAATTTGCATACTTTTTTTTGTTTATAACTGATAAGATAATTACATATTACAAATATTCATTTTATTAATATTTAAACTTAATc                        <  1:1688919/95‑1 (MQ=255)
                                                                         aTGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTTCAACTATTCACTTTATTAATATTTAAACTTAATCt                       >  1:1227952/1‑99 (MQ=255)
                                                                         aTGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                      >  2:2148130/1‑100 (MQ=255)
                                                                          aaaaaaACCATCCAAATCTGGATGGCTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                      >  1:142635/3‑100 (MQ=38)
                                                                          aaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                     >  1:305873/3‑100 (MQ=38)
                                                                          aaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                     >  2:887726/3‑100 (MQ=38)
                                                                          aaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                     >  1:195883/3‑100 (MQ=38)
                                                                          aaaaaaACCATCCAAATCTGGATGGTTTCTTTGGTCTATAACTGGTAAGATAATTACGTATTGCAAATTATCATTATATTAATGTTTACACTTAGTCta                      >  1:2736624/3‑99 (MQ=19)
                                                                           gAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                    >  2:2658862/1‑100 (MQ=255)
                                                                           gAAAAGCCATACAGATTTGGATGGTTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                    <  1:196886/100‑1 (MQ=255)
                                                                           aaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                    >  2:1650038/2‑100 (MQ=38)
                                                                             aaaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                  <  2:1955346/100‑1 (MQ=255)
                                                                             aaaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                  <  2:1955338/100‑1 (MQ=255)
                                                                              aaaccATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                 >  2:1327224/4‑100 (MQ=255)
                                                                              aaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                     >  1:2579536/1‑96 (MQ=255)
                                                                               aGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTATTCTATATgg                 >  1:15203/1‑99 (MQ=255)
                                                                                gCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGTaa               >  1:1704254/1‑100 (MQ=255)
                                                                                 ccATCCAAATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat              <  1:520134/100‑1 (MQ=255)
                                                                                    tcaaGATTTGATTTGTTTTTTTTTTTTATTGCTGTTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAGTCTATATg                  <  1:1211787/90‑1 (MQ=255)
                                                                                     ccAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata          <  2:1457086/100‑1 (MQ=255)
                                                                                      cAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat         <  1:1720171/100‑1 (MQ=255)
                                                                                      cAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat         >  2:2753286/1‑100 (MQ=255)
                                                                                        gATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatat       >  2:1039826/1‑100 (MQ=255)
                                                                                         aTTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATg      <  1:2524083/100‑1 (MQ=255)
                                                                                            tGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACa   <  2:308426/100‑1 (MQ=255)
                                                                                             ggATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAt  >  2:2792224/1‑100 (MQ=255)
                                                                                              gATGGTTTTTTTTTTTTTTTTCTGGTAAAATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATTACAt  >  1:946018/1‑99 (MQ=255)
                                                                                                   |                                                                                              
CCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAT  >  NC_000913/4409946‑4410139

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