Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,888,225:1 +A 17.0% 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.A17.0% 410.7 / 95.1 135intergenic (+32/‑211)mnmE/tnaCtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase/tryptophanase leader peptide
Reads supporting (aligned to +/- strand):  ref base . (62/50);  new base A (18/5);  total (80/55)
Fisher's exact test for biased strand distribution p-value = 6.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-03

Only 100 of 141 total aligned reads displayed.

CACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTT  >  NC_000913/3888126‑3888323
                                                                                                   |                                                                                                   
cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGa                                                                                                    <  1:897707/100‑1 (MQ=255)
cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGa                                                                                                    >  2:422135/1‑100 (MQ=255)
cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGa                                                                                                    >  1:2043243/1‑100 (MQ=255)
  ccGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATg                                                                                                  >  2:850919/1‑100 (MQ=255)
     gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCgg                                                                                               <  2:831930/100‑1 (MQ=38)
      ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgc                                                                                              <  2:1659702/100‑1 (MQ=255)
      ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgc                                                                                              <  2:1941844/100‑1 (MQ=38)
      ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgc                                                                                              <  1:2285812/100‑1 (MQ=38)
      ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgc                                                                                              <  2:1301400/100‑1 (MQ=38)
           ttACTTCAGACGACCTGCTGGGGCGTATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgc                                                                                         <  1:1513368/100‑1 (MQ=255)
              cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgcttc                                                                                      >  1:1981295/1‑100 (MQ=255)
               ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCa                                                                                     <  1:1300088/100‑1 (MQ=255)
               tgaaGACGACCTGCTGGGGTGGATTTTCTCCAGCTTCTGTTTTGGTACTTAACCGCGCTTACGAAGCCGCATTCTGTCTGTCAGATGCGGCTTCGCTTCa                                                                                     <  1:610319/97‑1 (MQ=255)
                  aGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTg                                                                                  >  1:363188/1‑100 (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)
                      gacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTAc                                                                              <  2:1713400/100‑1 (MQ=255)
                        ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGGTTCGCTTCATTGTTACCa                                                                            >  2:1182172/1‑100 (MQ=255)
                        ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCa                                                                            >  1:884999/1‑100 (MQ=255)
                           gctgGGGCGGATTTTCTCCAGCTTCTGTACTAGGAACTAACCGAGCTTACGACGCCGCATTCTGACTGTCAGATGCAGCTTCGCTTCATTGTTACCag                                                                           <  1:2162712/98‑2 (MQ=255)
                             tgGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCt                                                                       >  1:2414203/1‑100 (MQ=255)
                                ggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGtt                                                                    <  2:1264218/100‑1 (MQ=255)
                                ggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGtt                                                                    >  2:1441723/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)
                                    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)
                                         tctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCaa                                                           >  1:2129422/1‑100 (MQ=255)
                                          ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                          <  2:1130543/100‑1 (MQ=255)
                                          ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                          >  1:1873644/1‑100 (MQ=255)
                                           tcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGCAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                          >  1:1125454/1‑99 (MQ=255)
                                           tcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAcc                                                         >  2:219008/1‑100 (MQ=255)
                                             cAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCt                                                       >  1:2144414/1‑100 (MQ=255)
                                                   cTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCttttttt                                                 <  1:1239889/100‑1 (MQ=255)
                                                    tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                                >  2:2005277/1‑100 (MQ=255)
                                                    tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                                >  2:1113417/1‑100 (MQ=255)
                                                     gTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTaa                                               >  2:2135174/1‑100 (MQ=255)
                                                       aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                                >  1:2374931/1‑97 (MQ=255)
                                                       aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                             <  2:2339196/100‑1 (MQ=255)
                                                       aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                             <  2:631626/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:163721/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)
                                                                    cgcgCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCttn                                <  2:2018096/100‑2 (MQ=255)
                                                                    cgcgCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTa                                <  2:205691/100‑1 (MQ=255)
                                                                     gcgcTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTAc                               >  2:2584953/1‑100 (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)
                                                                          tccGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt                          <  2:1720013/98‑1 (MQ=255)
                                                                          tACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt                          >  1:689705/1‑100 (MQ=255)
                                                                          tACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt                          >  2:1555325/1‑100 (MQ=255)
                                                                          agcGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt                          >  2:2531900/3‑100 (MQ=37)
                                                                          agcGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt                          >  2:1633523/3‑100 (MQ=37)
                                                                          agcGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt                          >  1:1298864/3‑100 (MQ=37)
                                                                          agcGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt                          >  2:574080/3‑100 (MQ=37)
                                                                            cGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAttt                        >  2:1373248/1‑100 (MQ=255)
                                                                            cGAAGCCGCATTCTGACCGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAttt                        >  1:787537/1‑100 (MQ=255)
                                                                            cGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAttt                        >  1:2465590/1‑100 (MQ=37)
                                                                             gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGTTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta                       >  1:1687518/1‑100 (MQ=37)
                                                                             gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta                       >  2:1739574/1‑100 (MQ=38)
                                                                             gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta                       >  2:1966955/1‑100 (MQ=38)
                                                                             gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTAAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta                       >  2:1775091/1‑100 (MQ=37)
                                                                              aaGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTtat                      >  2:1210952/1‑100 (MQ=38)
                                                                               aGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTtata                     >  1:1641787/1‑100 (MQ=255)
                                                                               aGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTTTTTTTTTTAATTCTTACGTAATTtata                     >  1:1211507/1‑100 (MQ=21)
                                                                                gCCGCATTCTGACTGTCAGATGCGGCTTCGGTTGATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTtata                     >  1:1702234/1‑99 (MQ=255)
                                                                                gCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATaa                    >  2:1647551/1‑100 (MQ=255)
                                                                                gCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATaa                    >  2:748115/1‑100 (MQ=38)
                                                                                gCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATaa                    >  2:2074123/1‑100 (MQ=38)
                                                                                 ccGCATTCTGACTGTCAGATGCGGCTTCTCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATaa                    <  1:1059771/99‑1 (MQ=255)
                                                                                 ccGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAAt                   <  2:2346500/100‑1 (MQ=255)
                                                                                   cgcatCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATc                  >  2:1385946/5‑100 (MQ=38)
                                                                                     aTTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCttt               <  2:85028/100‑1 (MQ=255)
                                                                                        cTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaa            >  2:2372050/1‑100 (MQ=255)
                                                                                        cTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaa            >  2:72159/1‑100 (MQ=255)
                                                                                           aCTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaaaaa         >  2:848412/1‑100 (MQ=255)
                                                                                            cTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATTATCTTTaaaaaaa        >  1:1226603/1‑100 (MQ=255)
                                                                                             tGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAg       <  2:1984901/100‑1 (MQ=255)
                                                                                             tGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTATTTTTAAACATTAAAATGCTTACGTAATTTATAATCTTTaaaaaat        <  1:448204/99‑2 (MQ=255)
                                                                                              gTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGc      >  1:195587/1‑100 (MQ=255)
                                                                                               tCAGATGGCGATTCGCTTCATTGTTACGTCTCTTGTTATTTCTCAATTCTTTTTGTGAATTTTAGAATTGTTACTTGATTTATAATTTTTAAAAAAAGCa     <  1:210269/100‑1 (MQ=255)
                                                                                               tCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCa     <  1:2243851/100‑1 (MQ=255)
                                                                                                cAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCAt    >  1:272074/1‑100 (MQ=255)
                                                                                                  gATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCAttt  >  1:974110/1‑100 (MQ=255)
                                                                                                  gATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCAttt  >  1:2521458/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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: