Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,888,213 Δ1 bp 16.0% intergenic (+21/‑222) mnmE → / → tnaC tRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase/tryptophanase leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,888,2120T.16.0% 547.8 / 103.6 194intergenic (+20/‑223)mnmE/tnaCtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase/tryptophanase leader peptide
Reads supporting (aligned to +/- strand):  ref base T (81/82);  new base . (14/17);  total (95/99)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.88e-01

Only 100 of 198 total aligned reads displayed.

ACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTA  >  NC_000913/3888113‑3888311
                                                                                                   |                                                                                                     
aCTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCAt                                                                                                       >  1:361368/1‑100 (MQ=255)
   tAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCt                                                                                                    >  1:3550560/1‑100 (MQ=255)
    aaGCGAAATCACCGGGGAATTTACTAAAGAAGACCTGCTGGGGCGGATTTTATCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTg                                                                                                   >  2:92430/1‑100 (MQ=255)
      gagAAATCACCGGGGTATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGAc                                                                                                 <  1:3867823/98‑1 (MQ=255)
      gCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGATAAGTAACCGAGCTTACGAAGCCGCATTCTGAc                                                                                                 >  1:248974/1‑100 (MQ=255)
        gAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTg                                                                                               <  2:1639656/100‑1 (MQ=255)
        gAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTg                                                                                               <  2:1526750/100‑1 (MQ=255)
         aaaTCACCGGTGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGt                                                                                              <  1:3736740/100‑1 (MQ=255)
         aaaTCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGt                                                                                              >  2:2931510/1‑100 (MQ=255)
             cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGa                                                                                         >  2:3847167/1‑100 (MQ=255)
             cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGa                                                                                         >  2:2316204/1‑100 (MQ=255)
              aCCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTAAGAt                                                                                        >  1:1728180/1‑100 (MQ=255)
               ccGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGAt                                                                                        >  2:782241/1‑99 (MQ=255)
               ccGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATg                                                                                       >  2:101026/1‑100 (MQ=255)
                 ggggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCg                                                                                     <  1:1951896/100‑1 (MQ=255)
                  gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCgg                                                                                    <  1:3068317/100‑1 (MQ=38)
                  gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCgg                                                                                    <  1:2233479/100‑1 (MQ=38)
                   ggAATTTACTTCAGACGACCTGCTGGGTCGTATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgc                                                                                   <  1:396282/100‑1 (MQ=37)
                   ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgc                                                                                   <  2:2186982/100‑1 (MQ=38)
                   ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgc                                                                                   <  1:3506462/100‑1 (MQ=38)
                    gAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgct                                                                                  >  1:1840048/1‑100 (MQ=255)
                    gAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGCCAGAATGCGgct                                                                                  <  1:2582427/100‑1 (MQ=38)
                    gAAGTTACTTCAGACCACCTGCTGGGGCGGATTTTCTCCGGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgc                                                                                   <  2:1115530/99‑1 (MQ=37)
                     aaTTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgctt                                                                                 <  1:664733/100‑1 (MQ=38)
                     aaTTCACTTCAGACGATCTGCTGGGGCGAATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAATCCGCATCTGACAGTCAGAATGCGgc                                                                                   <  2:539371/98‑1 (MQ=37)
                       tttACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcg                                                                               <  2:1902161/100‑1 (MQ=255)
                        ttACTTCAGATGACCTGCTGGTGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgcttcgt                                                                              <  1:1867085/100‑2 (MQ=25)
                         tACTTCAGACGACCTGCTGGGGCGGATTTTTTCCATCTTTTTTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgct                                                                             <  1:458201/100‑1 (MQ=255)
                         tACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgct                                                                             <  1:2931536/100‑1 (MQ=255)
                          aCTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTTTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgctt                                                                            <  1:2152256/100‑1 (MQ=255)
                           cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCGGATGCGgcttcgcttc                                                                           >  1:3363077/1‑100 (MQ=255)
                           cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgcttc                                                                           >  1:1936560/1‑100 (MQ=255)
                           cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgcttc                                                                           >  2:1803406/1‑100 (MQ=255)
                           cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTATCAGATGCGgcttcgcttc                                                                           >  2:2452407/1‑100 (MQ=255)
                            ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCa                                                                          <  1:2115817/100‑1 (MQ=255)
                            ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCa                                                                          <  2:619356/100‑1 (MQ=255)
                            ttCAGACGACCTGCTGGGGCGGAATTTCTCCAGCTTCTGTATTGGTAAGTAACCGAGCTTACGAAGCCGCATGCTGACTGTCAGATGCGGCTTCGCTTCa                                                                          >  2:3141351/1‑100 (MQ=255)
                              cAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCAtt                                                                        >  1:3241626/1‑100 (MQ=255)
                               aGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTg                                                                       <  1:2843723/100‑1 (MQ=255)
                                  cgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTa                                                                    <  2:1114733/100‑1 (MQ=255)
                                    acCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTAcc                                                                  <  2:1936544/100‑1 (MQ=255)
                                     ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTAcc                                                                  <  1:2901863/99‑1 (MQ=255)
                                     ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCa                                                                 >  2:1786050/1‑100 (MQ=255)
                                     ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCa                                                                 >  1:2283822/1‑100 (MQ=255)
                                     ccTGCTGGGGCGGATTTTCTCCAGCTGCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCa                                                                 >  1:2387595/1‑100 (MQ=255)
                                          tgGGGCGGATTTTCTCCAGCTTCCGAATTGGTACGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGGTGCGGCTCCGCTTCATTGTTACCACTCCt                                                            >  2:2356186/1‑100 (MQ=255)
                                           ggggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTg                                                           >  2:1942290/1‑100 (MQ=255)
                                             ggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGtt                                                         >  2:1599202/1‑100 (MQ=255)
                                             ggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGtt                                                         <  2:1167362/100‑1 (MQ=255)
                                              gCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTa                                                        >  2:2430929/1‑100 (MQ=255)
                                               cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                       >  2:482638/1‑100 (MQ=255)
                                               cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                       >  1:623541/1‑100 (MQ=255)
                                               cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                       >  2:2773764/1‑100 (MQ=255)
                                               cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                       <  2:3265186/100‑1 (MQ=255)
                                                  aTTTTCTCCAGCTTCTGTATTGGTAAGTAATCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTcc                                                    <  1:1758217/100‑1 (MQ=255)
                                                  aTTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTcc                                                    >  1:3911156/1‑100 (MQ=255)
                                                  aTTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTcc                                                    <  2:361362/100‑1 (MQ=255)
                                                     ttCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCa                                                 <  2:2155098/100‑1 (MQ=255)
                                                     ctctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCa                                                 <  2:1843630/99‑1 (MQ=255)
                                                       ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                               <  1:2809632/100‑1 (MQ=255)
                                                       ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                               >  2:3830370/1‑100 (MQ=255)
                                                          catcTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCTCTCCTGTTATTCCTCAAc                                               <  2:1332129/94‑1 (MQ=255)
                                                                 tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                    >  2:853907/1‑100 (MQ=255)
                                                                    aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCTGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                 <  2:2742269/100‑1 (MQ=255)
                                                                    aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                 <  2:3496974/100‑1 (MQ=255)
                                                                    aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                 >  1:1577788/1‑100 (MQ=255)
                                                                     ttGGTAAGTAACCGCGGTTACGAAGCCTCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCATCCCTTTTTTTAAAc                                 <  2:3361580/99‑1 (MQ=255)
                                                                     ttGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACa                                >  2:363415/1‑100 (MQ=255)
                                                                       ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACAtt                              >  1:282009/1‑100 (MQ=255)
                                                                       ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACAtt                              >  2:1738988/1‑100 (MQ=255)
                                                                        gtaagtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTa                             >  1:130816/1‑100 (MQ=255)
                                                                            gtaaCCGCGCTTACGAAGCCGGATTCTGACTGTCAGATGCGGCTTCGCTTTATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAAt                         <  2:2339913/100‑1 (MQ=255)
                                                                             taaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTAAACCCTTTTTTTAAACATTAAAAtt                        >  1:1058170/1‑100 (MQ=255)
                                                                             taaCCGCGCTTACGAAGCCACATTCTGACAGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAAtt                        >  2:907382/1‑100 (MQ=255)
                                                                              aaCCGTGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTc                       <  1:856117/100‑1 (MQ=255)
                                                                                ccGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTTAAATTCt                      >  2:1232056/1‑99 (MQ=255)
                                                                                 cgcgCTTACGAAGCCGCATTCTGACTGTCAGCTGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTa                    >  1:1199902/1‑100 (MQ=255)
                                                                                  gcgcTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTa                    <  2:827291/99‑1 (MQ=255)
                                                                                  gcgcTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTAc                   <  2:2943804/100‑1 (MQ=255)
                                                                                   cgcTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACg                  <  1:545610/100‑1 (MQ=255)
                                                                                     cTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTa                >  1:754826/1‑100 (MQ=255)
                                                                                      ttATGATGCCGCATTCTGACTGTCAAATGCGGCTTCGCTTGATTGTTACCTCTCCTGTTATTAATTAACCCTTTTTTTAAAGATGAATATACTTACg                  <  2:3736567/97‑1 (MQ=255)
                                                                                      aagcGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTaa               >  1:2444461/4‑100 (MQ=37)
                                                                                       tACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt              >  1:3755482/1‑100 (MQ=255)
                                                                                       tACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt              >  2:609427/1‑100 (MQ=255)
                                                                                       tACGAAGCCGCATTCTGACTGTAAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTATTTTAAACATTAAAATTCTTACGTa                >  2:3008303/1‑99 (MQ=255)
                                                                                       agcGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt              >  1:1223715/3‑100 (MQ=37)
                                                                                         cGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAttt            >  1:569808/1‑100 (MQ=37)
                                                                                          gAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta           <  1:2730811/100‑1 (MQ=255)
                                                                                          gAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta           >  1:800505/1‑100 (MQ=255)
                                                                                          gAAGCCGCATCCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCATTCTTTTAAACATTAAAATTCTTACGTAATTta           <  2:2569236/100‑1 (MQ=255)
                                                                                          gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTTCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta           >  1:1058335/1‑100 (MQ=37)
                                                                                          gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta           >  1:1467840/1‑100 (MQ=38)
                                                                                          gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta           >  1:1262071/1‑100 (MQ=38)
                                                                                          gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta           >  2:2203782/1‑100 (MQ=38)
                                                                                             gCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATaa        >  1:3333425/1‑100 (MQ=38)
                                                                                             gCCGCATCTGACAGTCAGAATGCGGCTTCGCTTAATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACTTTAAAATTCTTACGTAATTTATaa        >  2:3744938/1‑100 (MQ=37)
                                                                                               cacATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAAt       >  2:1399935/3‑99 (MQ=255)
                                                                                                gtaTTCTGACTGTCAGATGCTGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAATCCATTTTTTAAACATTAAAATTCTTACGTAATTTATAATCt     <  2:3810075/98‑1 (MQ=255)
                                                                                                   ttCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTa  <  1:1098807/100‑1 (MQ=255)
                                                                                                   |                                                                                                     
ACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTA  >  NC_000913/3888113‑3888311

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: