Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,888,213 Δ1 bp 16.5% 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.5% 387.6 / 74.4 139intergenic (+20/‑223)mnmE/tnaCtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase/tryptophanase leader peptide
Reads supporting (aligned to +/- strand):  ref base T (61/55);  new base . (17/6);  total (78/61)
Fisher's exact test for biased strand distribution p-value = 6.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

Only 100 of 143 total aligned reads displayed.

CTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTT  >  NC_000913/3888114‑3888310
                                                                                                  |                                                                                                   
cTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCAtt                                                                                                    >  1:403046/1‑100 (MQ=255)
cTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCAtc                                                                                                    >  1:573484/1‑99 (MQ=255)
  tAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCt                                                                                                  >  2:446908/1‑100 (MQ=255)
    aGCGAAATCACCGGGGAATTTCCTTCAGACGACCTGCTGGGGCGGATGTCCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGa                                                                                                >  1:1269229/1‑100 (MQ=255)
    aGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGa                                                                                                >  1:2471193/1‑100 (MQ=255)
      cGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACt                                                                                              >  2:2072463/1‑100 (MQ=255)
      cGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACt                                                                                              <  2:2570405/100‑1 (MQ=255)
        aaaTCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGt                                                                                            <  1:887253/100‑1 (MQ=255)
        aaaTCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTc                                                                                           <  2:573486/100‑1 (MQ=255)
         aaTCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTc                                                                                           <  1:1991551/100‑1 (MQ=255)
         aaTCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTTGACAGTCa                                                                                          <  1:171218/100‑1 (MQ=255)
          aTCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCa                                                                                          <  2:1929209/100‑1 (MQ=255)
          aTCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCa                                                                                          <  2:342685/100‑1 (MQ=255)
            cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGa                                                                                       >  1:2043243/1‑100 (MQ=255)
            cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGa                                                                                       <  1:774469/100‑1 (MQ=255)
            cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGa                                                                                       >  2:422135/1‑100 (MQ=255)
              ccGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATg                                                                                     >  2:850919/1‑100 (MQ=255)
                ggggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGc                                                                                    <  1:1391883/99‑1 (MQ=255)
                 gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCgg                                                                                  <  2:831930/100‑1 (MQ=38)
                  ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgc                                                                                 <  2:1659702/100‑1 (MQ=255)
                  ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgc                                                                                 <  2:1301400/100‑1 (MQ=38)
                   gAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgct                                                                                <  2:1942075/100‑1 (MQ=38)
                             cAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCAt                                                                       >  1:2575238/1‑99 (MQ=255)
                              aGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTg                                                                     <  1:2318017/100‑1 (MQ=255)
                               gacgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGt                                                                    <  1:1672503/100‑1 (MQ=255)
                                acgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGtt                                                                   <  1:1753258/100‑1 (MQ=255)
                                 cgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAAGCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTa                                                                  >  1:1439999/1‑100 (MQ=255)
                                 cgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTa                                                                  >  1:1121338/1‑100 (MQ=255)
                                  gacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTan                                                                 <  2:1201844/100‑2 (MQ=255)
                                    ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGGTTCGCTTCATTGTTACCa                                                               >  2:1182172/1‑100 (MQ=255)
                                    ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCa                                                               >  1:884999/1‑100 (MQ=255)
                                       gctgGGGCGGATTTTCTCCAGCTTCTGTACTAGGAACTAACCGAGCTTACGACGCCGCATTCTGACTGTCAGATGCAGCTTCGCTTCATTGTTACCag                                                              <  1:2162712/98‑2 (MQ=255)
                                        ctgGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTcc                                                           <  2:1244490/100‑1 (MQ=255)
                                         tgGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCt                                                          >  1:2414203/1‑100 (MQ=255)
                                           gggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGt                                                        <  1:1274914/100‑1 (MQ=255)
                                            ggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGtt                                                       <  2:1264218/100‑1 (MQ=255)
                                            ggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGtt                                                       >  2:2243866/1‑100 (MQ=255)
                                             gCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGACGCATTCTGACTGTCAGATGAGGCTTCGCTTCATTGTTACCACTCCTGTTa                                                      >  1:1041745/1‑100 (MQ=255)
                                             gCGGATATTCTCCAGCTTCTGTACTCGTCCGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACCCCAGTTa                                                      >  1:2360929/1‑100 (MQ=255)
                                              cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                     >  1:599352/1‑100 (MQ=255)
                                               ggATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAtt                                                    >  1:2216298/1‑100 (MQ=255)
                                                gATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTTTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTc                                                   >  1:2098120/1‑100 (MQ=255)
                                                 aTTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTcc                                                  >  1:1520345/1‑100 (MQ=255)
                                                  ttttCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCt                                                 <  2:1028782/100‑1 (MQ=255)
                                                     tctcCGGCTTCTGTATGGGTAAGTAAGCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATGCCTCaa                                              >  1:1808358/1‑100 (MQ=255)
                                                      ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                             >  1:1873644/1‑100 (MQ=255)
                                                      ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                             <  2:1130543/100‑1 (MQ=255)
                                                       tcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGCAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                             >  1:1125454/1‑99 (MQ=255)
                                                         cAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCt                                          >  1:2144414/1‑100 (MQ=255)
                                                               cTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCttttttt                                    <  1:1239889/100‑1 (MQ=255)
                                                                tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                   >  2:1113417/1‑100 (MQ=255)
                                                                tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                   >  2:2005277/1‑100 (MQ=255)
                                                                 gTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTaa                                  >  2:2135174/1‑100 (MQ=255)
                                                                  tATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTaaa                                 <  2:237195/100‑1 (MQ=255)
                                                                   aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                <  2:2339196/100‑1 (MQ=255)
                                                                   aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                <  2:631626/100‑1 (MQ=255)
                                                                   aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                <  2:884090/100‑1 (MQ=255)
                                                                    ttGGTAAGTATCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACa                               <  2:2072898/100‑1 (MQ=255)
                                                                    ttGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACa                               >  2:2371392/1‑100 (MQ=255)
                                                                      ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACAtt                             >  1:1873623/1‑100 (MQ=255)
                                                                      ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACAtt                             >  1:2109580/1‑100 (MQ=255)
                                                                         aagtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTaaa                          >  2:1893482/1‑100 (MQ=255)
                                                                          agtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTaaaa                         >  2:110633/1‑100 (MQ=255)
                                                                           gtaaCCGCGCTTACGAAGCCGCCTTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAAt                        >  2:2389174/1‑100 (MQ=255)
                                                                           gtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTTAACATTAAAAt                        >  1:1025286/1‑100 (MQ=255)
                                                                              aCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCt                     >  2:1774277/1‑100 (MQ=255)
                                                                               ccGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCtt                    >  2:302509/1‑100 (MQ=255)
                                                                               ccGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCtt                    <  2:1440007/100‑1 (MQ=255)
                                                                                cgcgCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCttn                   <  2:2018096/100‑2 (MQ=255)
                                                                                 aatgaagcGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGTTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTAc                  >  2:1054827/8‑100 (MQ=21)
                                                                                  cgcTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACg                 <  2:1137775/100‑1 (MQ=255)
                                                                                  cgcTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACg                 <  2:1863165/100‑1 (MQ=255)
                                                                                  atgaagcGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACg                 >  1:653570/7‑100 (MQ=25)
                                                                                     ttACGAAGCCGCATTCTTACTGTCATGTGCGGCTTCGCTTTATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTaa              <  1:2448133/100‑1 (MQ=255)
                                                                                      tccGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt             <  2:1720013/98‑1 (MQ=255)
                                                                                      tACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt             >  1:689705/1‑100 (MQ=255)
                                                                                      tACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt             >  2:1555325/1‑100 (MQ=255)
                                                                                      agcGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt             >  1:1298864/3‑100 (MQ=37)
                                                                                      agcGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt             >  2:1633523/3‑100 (MQ=37)
                                                                                      agcGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt             >  2:2531900/3‑100 (MQ=37)
                                                                                        cGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAttt           >  2:1373248/1‑100 (MQ=255)
                                                                                        cGAAGCCGCATTCTGACTGTCACACGCAGCTTCGCTTCATTGTTACCACGCCTGTTATTCCTCAACCCTTTTTTTAAACCTTAAAATTCTTACGTaa              >  1:921044/1‑97 (MQ=255)
                                                                                        cGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAttt           >  2:206254/1‑100 (MQ=37)
                                                                                        cGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAttt           >  1:2465590/1‑100 (MQ=37)
                                                                                         gAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta          >  2:2115982/1‑100 (MQ=255)
                                                                                         gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGTTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta          >  1:1687518/1‑100 (MQ=37)
                                                                                         gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAATTATTAAAATTCTTACGTAATTta          >  1:1949987/1‑100 (MQ=37)
                                                                                         gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta          >  2:1966955/1‑100 (MQ=38)
                                                                                         gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta          >  2:1739574/1‑100 (MQ=38)
                                                                                         gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTAAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta          >  2:1775091/1‑100 (MQ=37)
                                                                                           aGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTtata        >  1:1641787/1‑100 (MQ=255)
                                                                                           aGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTTTTTTTTTTAATTCTTACGTAATTtata        >  1:1211507/1‑100 (MQ=21)
                                                                                            gCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATaa       >  2:1647551/1‑100 (MQ=255)
                                                                                            gCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATaa       >  2:748115/1‑100 (MQ=38)
                                                                                             ccGCATTCTGACTGTCAGATGCGGCTTCTCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATaa       <  1:1059771/99‑1 (MQ=255)
                                                                                             ccGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAAt      >  1:162409/1‑100 (MQ=255)
                                                                                               gaaTTCTGACTGTCAGAAGCGGCTTCAATTTATTGTTACCAATCCTGTTATTCCTAAACCTTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCt    <  1:478188/98‑1 (MQ=255)
                                                                                               gCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCt    <  2:1693502/100‑1 (MQ=255)
                                                                                                 aTTCTGACTGTCATTTGCGGCTTCCCTTCATTTTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCt    <  1:1731256/98‑1 (MQ=255)
                                                                                                 aTTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCttt  <  2:84906/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTT  >  NC_000913/3888114‑3888310

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