Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,410,037 T→C 15.6% intergenic (+31/‑96) 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,0370TC15.6% 269.9 / 45.8 131intergenic (+31/‑96)rlmB/yjfI23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base T (62/46);  new base C (13/7);  total (75/56)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-02

Only 100 of 132 total aligned reads displayed.

TTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATG  >  NC_000913/4409939‑4410135
                                                                                                  |                                                                                                  
tttCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGAttt                                                                                                   <  1:1212014/100‑1 (MQ=255)
 ttcttcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTg                                                                                                  <  2:1059702/100‑1 (MQ=255)
    ttcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGAt                                                                                               <  2:1560647/100‑1 (MQ=255)
    ttcCCTGAACGTTGCGGTTGCGACCGGAAGTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGAt                                                                                               <  1:1871236/100‑1 (MQ=255)
     tcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATg                                                                                              <  1:624930/100‑1 (MQ=255)
     tcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATg                                                                                              >  1:14691/1‑100 (MQ=255)
     tcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATg                                                                                              >  1:950512/1‑100 (MQ=255)
      cccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATgg                                                                                             <  2:511004/100‑1 (MQ=255)
       ccTGAACGTTTCGGTTGCGACCGGTATTTGCTTATTTGAAGCGGTTCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATgg                                                                                             >  2:2239673/1‑99 (MQ=255)
       ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATg                                                                                              >  2:2053622/1‑98 (MQ=255)
       ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGt                                                                                            >  2:574324/1‑100 (MQ=255)
        cTGAACGCTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGt                                                                                            >  1:488301/1‑99 (MQ=255)
          gAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGTAAAGCCATCCAAATCTGGATGGcttt                                                                                         <  2:1771939/100‑1 (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)
              gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttttt                                                                                     >  2:70394/1‑100 (MQ=255)
              gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcat                                                                                     <  1:190953/100‑4 (MQ=37)
               tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                    <  2:48488/100‑2 (MQ=255)
               tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                    <  2:2244938/100‑2 (MQ=255)
                 tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                  <  2:447613/100‑4 (MQ=255)
                   ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtct                                                                                 >  1:963378/1‑95 (MQ=255)
                   ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCta                                                                               >  2:1796737/1‑100 (MQ=255)
                   ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCta                                                                               >  2:1437996/1‑100 (MQ=255)
                   ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCta                                                                               >  2:2310102/1‑100 (MQ=255)
                   ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAAATCTGGATGGctttagatcggaa                                                                                >  1:415909/1‑91 (MQ=255)
                    gTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCTGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCtat                                                                              <  1:2047916/100‑1 (MQ=255)
                          aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTgg                                                                        <  2:128067/100‑1 (MQ=255)
                          aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTgg                                                                        <  2:1582758/100‑1 (MQ=255)
                           ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCAAATTTCTCAGAATTATGAAAAACCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGt                                                                       >  1:1801476/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)
                                   tGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAAt                                                               >  1:206049/1‑100 (MQ=255)
                                    gCTTATTTTGAGCGGTGGGCGAGCGCAGTTAATTTCTCAGAGTTAGGAAAAGGCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATa                                                                 <  1:2517599/97‑1 (MQ=255)
                                      ttATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTAc                                                            <  2:1925697/100‑1 (MQ=255)
                                       tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTTTAGCTGGTAAGATAATTACg                                                           >  1:2172063/1‑100 (MQ=255)
                                        aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                          >  2:1883916/1‑100 (MQ=255)
                                        aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                          >  1:1434223/1‑100 (MQ=255)
                                        aTTTGAAGCGGTGCGCCAGCGCAGCTAATCTCTCAGAATTATGAAAAGCCATCCAGATTAGGATGGTTTTTTTTGTCTTTAGTTGCTAAGATAATTACGt                                                          >  2:238215/1‑100 (MQ=255)
                                         tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTa                                                         >  2:272601/1‑100 (MQ=255)
                                         gttGAAACGGTGAGCCAGCGGAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATGTGCATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACgga                                                         <  1:2755125/99‑3 (MQ=255)
                                           tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                       <  2:2537585/100‑1 (MQ=255)
                                           tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                       >  2:142512/1‑100 (MQ=255)
                                           tGAAGCGGTGCGCCAGCACAGCTAATTTCTCAGAATTATGAAAAGCCATACAGCTTTGGCTGGTTTTTTTTTTTTCTATCTGGTAATATAATTa                                                             >  1:676784/1‑94 (MQ=255)
                                            gAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTg                                                      <  2:129877/100‑1 (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)
                                                       ccAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                           <  1:1596780/100‑1 (MQ=255)
                                                       ccAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                           <  2:667782/100‑1 (MQ=255)
                                                        cAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                          >  2:2852045/1‑100 (MQ=255)
                                                         aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                         <  1:2665957/100‑1 (MQ=255)
                                                         aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                         >  1:341676/1‑100 (MQ=255)
                                                         aGCGCAGCTAATGTGTCAGATTTATTAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                         <  2:537265/100‑1 (MQ=255)
                                                         aGCGCAGCGAATTTCTCAGTATTATGAAAAGCCATCCAGATGTGGCTGGTTGTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCa                                            <  1:2721946/97‑1 (MQ=255)
                                                          gcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtttt                                        >  2:2327624/1‑100 (MQ=255)
                                                          gcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtttt                                        <  2:1380493/100‑1 (MQ=255)
                                                            gcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                      <  2:1434228/100‑1 (MQ=255)
                                                              aGCTAATTTCTCTGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                    <  1:1228060/100‑1 (MQ=255)
                                                              aGCTAATTTCTCAGAATTATGAAGAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                     <  1:2363732/99‑1 (MQ=255)
                                                               gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTaa                                   <  1:1335479/100‑1 (MQ=255)
                                                                 tAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                 <  2:558731/100‑1 (MQ=255)
                                                                   aTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAttt                              >  2:1948318/1‑100 (MQ=255)
                                                                      tctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTTTTTTTCTGCTAAGATAATTCCGTATTGCAAATATTCATTTTATTAATATTTaa                            >  1:465221/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)
                                                                          aGACTTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                          <  2:881109/98‑1 (MQ=255)
                                                                             gacaaaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                     >  1:368658/7‑100 (MQ=37)
                                                                                attaaaaACCATCCAAATTTGCATACTTTTTTTTGTTTATAACTGATAAGATAATTACATATTACAAATATTCATTTTATTAATATTTAAACTTAATc                    <  1:1688919/95‑1 (MQ=255)
                                                                                aTGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTTCAACTATTCACTTTATTAATATTTAAACTTAATCt                   >  1:1227952/1‑99 (MQ=255)
                                                                                 aaaaaaACCATCCAAATCTGGATGGCTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                  >  1:142635/3‑100 (MQ=38)
                                                                                 aaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                 >  2:887726/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)
                                                                                    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)
                                                                                       accATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaa           >  2:277448/2‑100 (MQ=255)
                                                                                        ccATCCAAATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat          <  1:520134/100‑1 (MQ=255)
                                                                                          aTCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCTTTTTATTAATATTTAAACTTAATTTATATGAtataa        >  1:2561960/1‑100 (MQ=255)
                                                                                           tcaaGATTTGATTTGTTTTTTTTTTTTATTGCTGTTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAGTCTATATg              <  1:1211787/90‑1 (MQ=255)
                                                                                            ccAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata      <  2:1457086/100‑1 (MQ=255)
                                                                                             cAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat     >  2:2753286/1‑100 (MQ=255)
                                                                                             cAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat     <  1:1720171/100‑1 (MQ=255)
                                                                                               gATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatat   >  2:1039826/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
TTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATG  >  NC_000913/4409939‑4410135

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: