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)
                   ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCta                                                                               >  2:1796737/1‑100 (MQ=255)
                   ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCta                                                                               >  2:1437996/1‑100 (MQ=255)
                   ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCta                                                                               >  2:2310102/1‑100 (MQ=255)
                   ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAAATCTGGATGGctttagatcggaa                                                                                >  1:415909/1‑91 (MQ=255)
                    gTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCTGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCtat                                                                              <  1:2047916/100‑1 (MQ=255)
                          aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTgg                                                                        <  2:128067/100‑1 (MQ=255)
                          aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTgg                                                                        <  2:1582758/100‑1 (MQ=255)
                           ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCAAATTTCTCAGAATTATGAAAAACCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGt                                                                       >  1:1801476/1‑100 (MQ=255)
                            cGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTa                                                                      >  2:900118/1‑100 (MQ=255)
                            cGGAATTTGCTTATTTGAAGCGGTGAGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGt                                                                       >  1:708466/1‑99 (MQ=255)
                             ggAATTTGCTTAATTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTaa                                                                     >  2:1919499/1‑100 (MQ=255)
                              gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAg                                                                    <  2:2133691/100‑1 (MQ=255)
                                  ttGCTTATTTGAAGCGGTGCTCCAGCGCAGCTAATTTCTCAGAATTAGGAAAAGCCATCCATATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATaa                                                                <  2:2432645/100‑1 (MQ=255)
                                   tGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAAt                                                               >  1:206049/1‑100 (MQ=255)
                                    gCTTATTTTGAGCGGTGGGCGAGCGCAGTTAATTTCTCAGAGTTAGGAAAAGGCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATa                                                                 <  1:2517599/97‑1 (MQ=255)
                                      ttATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTAc                                                            <  2:1925697/100‑1 (MQ=255)
                                       tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTTTAGCTGGTAAGATAATTACg                                                           >  1:2172063/1‑100 (MQ=255)
                                        aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                          >  2:1883916/1‑100 (MQ=255)
                                        aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                          >  1:1434223/1‑100 (MQ=255)
                                        aTTTGAAGCGGTGCGCCAGCGCAGCTAATCTCTCAGAATTATGAAAAGCCATCCAGATTAGGATGG‑TTTTTTTTGTCTTTAGTTGCTAAGATAATTACGt                                                          >  2:238215/1‑100 (MQ=255)
                                         tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTa                                                         >  2:272601/1‑100 (MQ=255)
                                         gttGAAACGGTGAGCCAGCGGAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATGTGCATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACgga                                                         <  1:2755125/99‑3 (MQ=255)
                                           tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                       <  2:2537585/100‑1 (MQ=255)
                                           tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                       >  2:142512/1‑100 (MQ=255)
                                           tGAAGCGGTGCGCCAGCACAGCTAATTTCTCAGAATTATGAAAAGCCATACAGCTTTGGCTGGTTTTTTTTT‑TTTCTATCTGGTAATATAATTa                                                             >  1:676784/1‑94 (MQ=255)
                                            gAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTg                                                      <  2:129877/100‑1 (MQ=255)
                                             aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                     >  2:1289910/1‑100 (MQ=255)
                                              aGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCa                                                    <  2:1177585/100‑1 (MQ=255)
                                                cGGTGCGCCAGCGCAGCTAAGTTCTCAGAATTATGAAAAGCCGTCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAATATAATTACGTATTGCaaa                                                  <  2:330489/100‑1 (MQ=255)
                                                 ggTGCGCCAGCGCGGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTT‑TATATAGCTGGTAAGATAATTACGTATTGCAAat                                                 >  1:2426916/1‑100 (MQ=255)
                                                 ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTT‑TTTATTGTTGGTAAGATAATTACGTATTTCAAat                                                 >  1:2485228/1‑100 (MQ=255)
                                                 ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAat                                                 >  2:1708887/1‑100 (MQ=255)
                                                 ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAat                                                 >  2:2305258/1‑100 (MQ=255)
                                                    gcgcCAGCGCAGCTAACTAATCAGAATAATGAACAGCCAGCGAGATTTGCATGG‑TCTATTTTGTCGATAGCTGTTAAGATAATTACGTATTACaaacact                                              >  1:2309218/1‑96 (MQ=255)
                                                      gcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTc                                             >  1:917332/1‑99 (MQ=255)
                                                       ccAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                           <  1:1596780/100‑1 (MQ=255)
                                                       ccAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                           <  2:667782/100‑1 (MQ=255)
                                                        cAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                          >  2:2852045/1‑100 (MQ=255)
                                                         aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                         <  1:2665957/100‑1 (MQ=255)
                                                         aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                         >  1:341676/1‑100 (MQ=255)
                                                         aGCGCAGCTAATGTGTCAGATTTATTAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                         <  2:537265/100‑1 (MQ=255)
                                                         aGCGCAGCGAATTTCTCAGTATTATGAAAAGCCATCCAGATGTGGCTGG‑TTGTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCa                                            <  1:2721946/97‑1 (MQ=255)
                                                          gcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtttt                                        >  2:2327624/1‑100 (MQ=255)
                                                          gcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtttt                                        <  2:1380493/100‑1 (MQ=255)
                                                            gcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                      <  2:1434228/100‑1 (MQ=255)
                                                              aGCTAATTTCTCTGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                    <  1:1228060/100‑1 (MQ=255)
                                                              aGCTAATTTCTCAGAATTATGAAGAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                     <  1:2363732/99‑1 (MQ=255)
                                                               gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTaa                                   <  1:1335479/100‑1 (MQ=255)
                                                                 tAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                 <  2:558731/100‑1 (MQ=255)
                                                                   aTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAttt                              >  2:1948318/1‑100 (MQ=255)
                                                                      tctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTT‑TTTTTTTCTGCTAAGATAATTCCGTATTGCAAATATTCATTTTATTAATATTTaa                            >  1:465221/1‑100 (MQ=255)
                                                                       ctcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                           >  2:2522640/1‑100 (MQ=255)
                                                                        tcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                          >  2:1645886/1‑100 (MQ=255)
                                                                        tcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTAATTTTATTAATATTTAAAc                          >  1:2666564/1‑100 (MQ=255)
                                                                          aGACTTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                          <  2:881109/98‑1 (MQ=255)
                                                                             gacaaaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                     >  1:368658/7‑100 (MQ=37)
                                                                                attaaaaACCATCCAAATTTGCAT‑ACTTTTTTTTGTTTATAACTGATAAGATAATTACATATTACAAATATTCATTTTATTAATATTTAAACTTAATc                    <  1:1688919/95‑1 (MQ=255)
                                                                                aTGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTTCAACTATTCACTTTATTAATATTTAAACTTAATCt                   >  1:1227952/1‑99 (MQ=255)
                                                                                 aaaaaaACCATCCAAATCTGGATGGCTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                  >  1:142635/3‑100 (MQ=38)
                                                                                 aaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                 >  2:887726/3‑100 (MQ=38)
                                                                                 aaaaaaACCATCCAAATCTGGATGG‑TTTCTTTGGTCTATAACTGGTAAGATAATTACGTATTGCAAATTATCATTATATTAATGTTTACACTTAGTCta                  >  1:2736624/3‑99 (MQ=19)
                                                                                  gAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                >  2:2658862/1‑100 (MQ=255)
                                                                                  gAAAAGCCATACAGATTTGGATGG‑TTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                <  1:196886/100‑1 (MQ=255)
                                                                                    aaaGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg              <  2:1955346/100‑1 (MQ=255)
                                                                                    aaaGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg              <  2:1955338/100‑1 (MQ=255)
                                                                                     aaaccATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg             >  2:1327224/4‑100 (MQ=255)
                                                                                     aaGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                 >  1:2579536/1‑96 (MQ=255)
                                                                                      aGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTATTCTATATgg             >  1:15203/1‑99 (MQ=255)
                                                                                       gCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGTaa           >  1:1704254/1‑100 (MQ=255)
                                                                                       accATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaa           >  2:277448/2‑100 (MQ=255)
                                                                                        ccATCCAAATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat          <  1:520134/100‑1 (MQ=255)
                                                                                          aTCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCTTTTTATTAATATTTAAACTTAATTTATATGAtataa        >  1:2561960/1‑100 (MQ=255)
                                                                                           tcaaGATTTGATTTGTTTTTTTTT‑TTTATTGCTGTTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAGTCTATATg              <  1:1211787/90‑1 (MQ=255)
                                                                                            ccAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata      <  2:1457086/100‑1 (MQ=255)
                                                                                             cAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat     >  2:2753286/1‑100 (MQ=255)
                                                                                             cAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat     <  1:1720171/100‑1 (MQ=255)
                                                                                               gATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatat   >  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: ‑