Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,888,225:1 +A 15.5% intergenic (+33/‑210) 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,2241.A15.5% 693.1 / 141.3 213intergenic (+32/‑211)mnmE/tnaCtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase/tryptophanase leader peptide
Reads supporting (aligned to +/- strand):  ref base . (93/87);  new base A (17/16);  total (110/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 216 total aligned reads displayed.

CACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTT  >  NC_000913/3888126‑3888323
                                                                                                   |                                                                                                   
    ggggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCg                                                                                                >  1:2217171/1‑100 (MQ=255)
    ggggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCg                                                                                                >  1:967879/1‑100 (MQ=255)
     gggTATTTACTTCAGACGACCTGCTGGGGCGGATGTTCTCCCGCTTCTGTATTGGTAAGTAACCGTGCTTACGAAGCCGCATCTGACAGTCAGAATGCgg                                                                                               <  2:1504395/100‑1 (MQ=37)
      ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgc                                                                                              <  2:1049282/100‑1 (MQ=38)
       gAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgct                                                                                             >  2:3374326/1‑100 (MQ=255)
       gAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCgg                                                                                               <  2:2878180/98‑1 (MQ=38)
       gAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCgg                                                                                               <  2:2983428/98‑1 (MQ=38)
          tttACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgcttcg                                                                                          <  2:3228635/100‑1 (MQ=37)
             aCTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgcttcgtaa                                                                                       <  1:3071635/100‑4 (MQ=37)
             aCTTCAGACGACCTGCTGGGGCGGATTTTCCCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgctt                                                                                       <  1:3549971/100‑1 (MQ=255)
              cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGACGCCGCATTCTGACTGTCAGATGCGgcttcgcttc                                                                                      >  1:3614482/1‑100 (MQ=255)
              cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgcttc                                                                                      <  2:1657332/100‑1 (MQ=255)
               ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCa                                                                                     <  2:3265853/100‑1 (MQ=255)
               ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCa                                                                                     <  2:552321/100‑1 (MQ=255)
                tCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCAt                                                                                    <  2:3323051/100‑1 (MQ=255)
                tCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCAt                                                                                    >  1:3410064/1‑100 (MQ=255)
                tCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCAt                                                                                    >  1:1210126/1‑100 (MQ=255)
                tCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTGCTGTATTGGTAAGAAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTGCGCTTCa                                                                                     <  2:791615/99‑1 (MQ=255)
                 cAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCAtt                                                                                   <  1:261208/100‑1 (MQ=255)
                 cAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCAtt                                                                                   <  2:1743852/100‑1 (MQ=255)
                  aGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTg                                                                                  >  1:3393555/1‑100 (MQ=255)
                  aGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTg                                                                                  <  2:819335/100‑1 (MQ=255)
                  aGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTg                                                                                  <  1:1442336/100‑1 (MQ=255)
                      gacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTAc                                                                              <  2:2838069/100‑1 (MQ=255)
                      gacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTAc                                                                              <  2:3276547/100‑1 (MQ=255)
                       acCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTAcc                                                                             <  1:1545612/100‑1 (MQ=255)
                       acCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTAcc                                                                             <  2:1629606/100‑1 (MQ=255)
                       acCTGCTGGGGCGGATTCTCTCCAGAATCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATGGTTatc                                                                             >  2:1592260/1‑98 (MQ=255)
                         ctgctgGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCAc                                                                           >  1:812742/1‑100 (MQ=255)
                           gctgGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTc                                                                         <  1:232134/100‑1 (MQ=255)
                             tgGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCt                                                                       <  1:748006/100‑1 (MQ=255)
                              ggggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTg                                                                      <  2:3050532/100‑1 (MQ=255)
                               gggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGt                                                                     <  1:722498/100‑1 (MQ=255)
                               gggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGt                                                                     <  2:183054/100‑1 (MQ=255)
                                ggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGtt                                                                    >  1:2741132/1‑100 (MQ=255)
                                ggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGtt                                                                    <  1:2846776/100‑1 (MQ=255)
                                  cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                                  <  1:1894985/100‑1 (MQ=255)
                                  cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                                  >  1:1604876/1‑100 (MQ=255)
                                  cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                                  >  2:1339969/1‑100 (MQ=255)
                                    gATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTc                                                                >  1:676371/1‑100 (MQ=255)
                                     aTTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTcc                                                               <  2:2399150/100‑1 (MQ=255)
                                     aTTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTcc                                                               >  2:2047198/1‑100 (MQ=255)
                                      ttttCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCt                                                              <  1:2180757/100‑1 (MQ=255)
                                         tctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCaa                                                           <  1:3680964/100‑1 (MQ=255)
                                          ctgcAGCTTCTGTATTGGTAAGGAACCGCGCTTACGAAGCCGCATTCCGACTGTCAGATGCGGCTTCGCTTCAGTGTTACCACCCCTGTTATTCCTCaa                                                           <  2:1227168/96‑1 (MQ=255)
                                          ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                          <  1:1106128/100‑1 (MQ=255)
                                           tcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAcc                                                         >  1:2554188/1‑100 (MQ=255)
                                            ccAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAccc                                                        >  1:1743573/1‑100 (MQ=255)
                                            ccAGCCTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAccc                                                        >  1:3776026/1‑100 (MQ=255)
                                             cAGCTTCTGTATTGGTAAGTAACCGCGCTTACGCAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCAGTGTTACCACTCCTGTTATTCCTCAACCCt                                                       >  2:2297462/1‑100 (MQ=255)
                                             cAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCt                                                       >  1:1871824/1‑100 (MQ=255)
                                             cAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCt                                                       >  1:2822490/1‑100 (MQ=255)
                                                cTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCtttt                                                    <  2:3302729/100‑1 (MQ=255)
                                                  tCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCtttttt                                                  <  1:3313352/100‑1 (MQ=255)
                                                    tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                                >  2:1204753/1‑100 (MQ=255)
                                                    tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                                <  1:1959617/100‑1 (MQ=255)
                                                    tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                                <  1:2388167/100‑1 (MQ=255)
                                                     gTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTaa                                               >  2:1401841/1‑100 (MQ=255)
                                                     gTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTaa                                               <  1:3184456/100‑1 (MQ=255)
                                                       aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                             <  1:1064573/100‑1 (MQ=255)
                                                        ttGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACa                                            <  2:2182071/100‑1 (MQ=255)
                                                         aggTAAGTACCCGCGCTTACAAAGCCGCATTCTGACTGTCAGATGAGTCTTCGCTTCATTGTTACCACTCCTGTTATTTCTCAACTCTTTTTTTAAAAAt                                           >  2:2489795/2‑100 (MQ=255)
                                                          ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTAGCACTCCTGTTATTCCTCAACCCTTTTTTTAAACAtt                                          >  1:2683277/1‑100 (MQ=255)
                                                          ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACAtt                                          >  1:2460936/1‑100 (MQ=255)
                                                          ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACAtt                                          >  1:2030617/1‑100 (MQ=255)
                                                          ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACAtt                                          >  1:3229889/1‑100 (MQ=255)
                                                          ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACAtt                                          >  1:187892/1‑100 (MQ=255)
                                                          ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGGTATTCCTCAACCCTTTTTTTAAATAtt                                          >  2:356149/1‑100 (MQ=255)
                                                           gtaagtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTa                                         >  1:1803157/1‑100 (MQ=255)
                                                               gtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAAt                                     >  2:695744/1‑100 (MQ=255)
                                                               gtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCCGCTTCGCTTTATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAAt                                     <  2:2710088/100‑1 (MQ=255)
                                                                 aaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTc                                   <  1:1675689/100‑1 (MQ=255)
                                                                  aCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCt                                  >  1:200398/1‑100 (MQ=255)
                                                                  aCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCt                                  >  1:1380530/1‑100 (MQ=255)
                                                                   ccGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCtt                                 >  1:2127349/1‑100 (MQ=255)
                                                                    cgcgCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTa                                <  1:2759719/100‑1 (MQ=255)
                                                                     gcgcTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTTAACCCTTTTTTTAAACATTAAAATTCTTa                                >  2:1762583/1‑99 (MQ=255)
                                                                      cgcTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACg                              <  2:599272/100‑1 (MQ=255)
                                                                       gcTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGt                             >  1:1475380/1‑100 (MQ=255)
                                                                           aCGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAtt                         <  2:2253897/100‑1 (MQ=255)
                                                                            cGATGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAttt                        <  2:3481068/100‑1 (MQ=255)
                                                                            cGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAttt                        >  2:1689512/1‑100 (MQ=37)
                                                                             gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCGCTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTa                            >  2:2119634/1‑95 (MQ=37)
                                                                             gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta                       >  2:876498/1‑100 (MQ=38)
                                                                             gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta                       >  1:3039171/1‑100 (MQ=38)
                                                                             gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta                       >  1:2888419/1‑100 (MQ=38)
                                                                              aaGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTtat                      >  2:2640103/1‑100 (MQ=38)
                                                                               aGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTtata                     >  1:392561/1‑100 (MQ=38)
                                                                               aGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTtata                     >  1:2914411/1‑100 (MQ=38)
                                                                                gCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATaa                    >  1:2283806/1‑100 (MQ=255)
                                                                                       tCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaa             <  2:373332/100‑1 (MQ=255)
                                                                                         tGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaaa           >  2:3784680/1‑100 (MQ=255)
                                                                                            cTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAATATTAAAATTCTTACGTAATTTATAATCTTTaaaaaaa        >  2:62331/1‑100 (MQ=255)
                                                                                            cTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaaaaaa        >  2:3501873/1‑100 (MQ=255)
                                                                                              gTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGc      <  1:1823158/100‑1 (MQ=255)
                                                                                              gTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGc      >  1:1410827/1‑100 (MQ=255)
                                                                                               tCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCa     <  1:3602369/100‑1 (MQ=255)
                                                                                                cAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCAt    >  1:2663918/1‑100 (MQ=255)
                                                                                                  gATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCAttt  >  1:323866/1‑100 (MQ=255)
                                                                                                  gATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCAttt  >  1:79223/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTT  >  NC_000913/3888126‑3888323

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