Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,888,213 Δ1 bp 15.2% 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.15.2% 456.7 / 80.3 158intergenic (+20/‑223)mnmE/tnaCtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase/tryptophanase leader peptide
Reads supporting (aligned to +/- strand):  ref base T (80/53);  new base . (9/15);  total (90/68)
Fisher's exact test for biased strand distribution p-value = 4.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.55e-01

Only 100 of 165 total aligned reads displayed.

ACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTA  >  NC_000913/3888113‑3888311
                                                                                                   |                                                                                                    
aCTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCAt                                                                                                      >  2:1295786/1‑100 (MQ=255)
aCTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCAt                                                                                                      >  2:1317706/1‑100 (MQ=255)
aCTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCGCCAGCTTATGTATTGGTAAGTAACCGCGCTTACAAAGCCGCAt                                                                                                      >  2:3120837/1‑100 (MQ=255)
 cTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCAtt                                                                                                     >  1:1583559/1‑100 (MQ=255)
 cTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCAt                                                                                                      <  1:1184724/99‑1 (MQ=255)
  ttAAGCGAAATAACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGAATTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCACATTc                                                                                                    >  2:1653634/1‑100 (MQ=255)
      gCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGAc                                                                                                >  1:806482/1‑100 (MQ=255)
      gCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGAc                                                                                                >  1:1322563/1‑100 (MQ=255)
       cGAAATCACCGGGGAATTTACTTCATACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGTTAAGTAACCGCGCTTACTAAGCCGCATTCTGAc                                                                                                <  2:1471875/99‑1 (MQ=255)
       cGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACt                                                                                               >  1:1608778/1‑100 (MQ=255)
       cGAAATCACCGGGGAATTTACTTCAGACGACCGGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGCAGCCGCATTCTGAc                                                                                                <  2:1216553/99‑1 (MQ=255)
         aaaTCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGt                                                                                             >  2:2497969/1‑100 (MQ=255)
          aaTCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTc                                                                                            <  2:3598617/100‑1 (MQ=255)
          aaTCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTc                                                                                            <  1:1012483/100‑1 (MQ=255)
             cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAg                                                                                          <  1:3320977/99‑1 (MQ=255)
             cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑a                                                                                        <  1:2380386/100‑1 (MQ=255)
             cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAG‑aa                                                                                       <  1:2601210/100‑2 (MQ=255)
                 ggggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCg                                                                                    <  2:1407318/100‑1 (MQ=255)
                 ggggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCg                                                                                    <  1:1961106/100‑1 (MQ=255)
                  gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCgg                                                                                   >  1:1019534/1‑100 (MQ=255)
                  gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCgg                                                                                   <  2:2361492/100‑1 (MQ=38)
                  gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCgg                                                                                   <  1:523320/100‑1 (MQ=38)
                  gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCgg                                                                                   <  1:424760/100‑1 (MQ=38)
                   ggAATTTACTTCAGAGGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgc                                                                                  <  1:1618524/100‑1 (MQ=38)
                   ggAATTTACTTCAGACGACCTGCTGGTGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgc                                                                                  <  2:2314521/100‑1 (MQ=38)
                   ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgc                                                                                  <  1:2323887/100‑1 (MQ=255)
                   ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgc                                                                                  <  1:279968/100‑1 (MQ=255)
                   ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgc                                                                                  <  1:1369740/100‑1 (MQ=38)
                   ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgc                                                                                  <  1:3136156/100‑1 (MQ=38)
                   ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgc                                                                                  <  1:1151031/100‑1 (MQ=38)
                    gAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgct                                                                                 <  2:2765590/100‑1 (MQ=38)
                     aaTTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgctt                                                                                <  2:278420/100‑1 (MQ=255)
                       tttACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcg                                                                              <  2:31053/100‑1 (MQ=255)
                        ttACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACTAAACCGCATTCTGACTGTCAG‑ATGCGgcttcgc                                                                             >  2:2536936/1‑100 (MQ=255)
                         tACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgct                                                                            <  2:3574735/100‑1 (MQ=255)
                         tACTTCAGACGACCTGCTGGCGCGGATTTTCCCCAGCTTCTGTATTGGTAAGTAACCGCTCTTACGAAGCCGCA‑TATGACAGTCAGAATGCGgcttcg                                                                              <  2:2886251/98‑1 (MQ=25)
                          aCTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgctt                                                                           <  2:2158747/100‑1 (MQ=255)
                           cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgcttc                                                                          >  1:2602625/1‑100 (MQ=255)
                           cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgcttc                                                                          >  1:1989269/1‑100 (MQ=255)
                           cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgcttc                                                                          >  2:534732/1‑100 (MQ=255)
                           cTTCAGACGACCTGCTCGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgcttc                                                                          <  2:2839434/100‑1 (MQ=255)
                            ttCAGACGACATGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGTCGCATTCTGGCTGTCAG‑ATGCGGCTTCGCTTCa                                                                         <  2:2399033/100‑1 (MQ=255)
                              cAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCAtt                                                                       >  1:2406602/1‑100 (MQ=255)
                               aGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTg                                                                      <  1:751881/100‑1 (MQ=255)
                               aGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTg                                                                      >  1:1423972/1‑100 (MQ=255)
                                gacgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGt                                                                     <  2:505231/100‑1 (MQ=255)
                                 acgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGtt                                                                    <  1:257975/100‑1 (MQ=255)
                                  cgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTa                                                                   >  2:1837698/1‑100 (MQ=255)
                                    acCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTAcc                                                                 >  1:2873871/1‑100 (MQ=255)
                                    acCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTAcc                                                                 <  1:2554057/100‑1 (MQ=255)
                                     ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCa                                                                >  1:934760/1‑100 (MQ=255)
                                     ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCa                                                                >  2:1120565/1‑100 (MQ=255)
                                       tgctgGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACt                                                              >  2:2381790/1‑100 (MQ=255)
                                        gctgGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTc                                                             >  1:901412/1‑100 (MQ=255)
                                            gggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGt                                                         <  1:1965899/100‑1 (MQ=255)
                                            gggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGt                                                         <  2:1925086/100‑1 (MQ=255)
                                              gCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTa                                                       >  2:1540276/1‑100 (MQ=255)
                                                ggATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAtt                                                     >  1:634837/1‑100 (MQ=255)
                                                   ttttCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCt                                                  >  1:2232312/1‑100 (MQ=255)
                                                   ttttCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCt                                                  <  1:336742/100‑1 (MQ=255)
                                                       ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                              >  2:459422/1‑100 (MQ=255)
                                                       ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                              >  1:2912918/1‑100 (MQ=255)
                                                        tcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAcc                                             >  1:2554375/1‑100 (MQ=255)
                                                          cAGCTTCAGTAGTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGGTTCATTGTTACCACTCCTGTTATTCCTCAACCCt                                           >  2:1910756/1‑100 (MQ=255)
                                                                 tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                    >  1:1305761/1‑100 (MQ=255)
                                                                   tATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTaaa                                  >  1:1086834/1‑100 (MQ=255)
                                                                   tATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTaaa                                  >  1:3264874/1‑100 (MQ=255)
                                                                    aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTGCCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                 >  1:1626679/1‑100 (MQ=255)
                                                                    aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                 <  2:1719116/100‑1 (MQ=255)
                                                                       ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACAtt                              >  2:3526886/1‑100 (MQ=255)
                                                                       ggTAAGTAACCGCGCTTACAAAGACGAATTCAGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGGTATTCCTCAACCCTTTTTTTaaaa                                 >  2:2202393/1‑96 (MQ=255)
                                                                         taagtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTaa                            >  1:754878/1‑100 (MQ=255)
                                                                           agtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTaaaa                          <  1:2708055/100‑1 (MQ=255)
                                                                           agtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTaaaa                          <  2:1941306/100‑1 (MQ=255)
                                                                             taaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAAtt                        >  2:1493952/1‑100 (MQ=255)
                                                                              aaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTc                       <  1:2171503/100‑1 (MQ=255)
                                                                                ccGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTACAATTCtt                     <  1:1916878/100‑1 (MQ=255)
                                                                                ccGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCtt                     >  1:3417894/1‑100 (MQ=255)
                                                                                  gcgcTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTAc                   >  1:1890404/1‑100 (MQ=255)
                                                                                  gcgcTTACGAAGCCGCATTCTGACTGATAG‑ATGTGGCTTCCCTTCATTATTACCACTCCTATTATTCCTCAACCCTTTTTTTAAACATTAAAATTCt                      <  2:3448971/97‑1 (MQ=255)
                                                                                   cgcTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACg                  >  1:1854030/1‑100 (MQ=255)
                                                                                     cTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTa                >  1:2613936/1‑100 (MQ=255)
                                                                                       tACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAACATTCTTACGTAAt              >  2:1743951/1‑100 (MQ=255)
                                                                                       tACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt              >  1:1411623/1‑100 (MQ=255)
                                                                                       agcGAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt              >  1:3556698/3‑100 (MQ=37)
                                                                                        gcGAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAtt             >  1:2876379/2‑100 (MQ=37)
                                                                                        aCGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAtt             >  2:3295409/1‑100 (MQ=255)
                                                                                         cGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAttt            >  2:2895086/1‑100 (MQ=255)
                                                                                         cGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACCTAAt              >  2:773576/1‑98 (MQ=255)
                                                                                         cGAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTAc                   >  2:288864/1‑93 (MQ=37)
                                                                                          gAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta           >  1:550884/1‑100 (MQ=255)
                                                                                          gAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta           >  2:1863111/1‑100 (MQ=38)
                                                                                           aaGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTtat          <  1:1787286/100‑1 (MQ=255)
                                                                                           aaGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTtat          >  2:1423520/1‑100 (MQ=255)
                                                                                            aGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTtata         >  1:1621993/1‑100 (MQ=255)
                                                                                               cGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATc      >  1:3316935/1‑100 (MQ=255)
                                                                                                 cATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCtt    >  1:1530686/1‑100 (MQ=255)
                                                                                                  aTTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTACTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCttt   <  2:2090477/100‑1 (MQ=255)
                                                                                                   ttCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTa  <  1:2758021/100‑1 (MQ=255)
                                                                                                   ttCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTa  <  2:3605912/100‑1 (MQ=255)
                                                                                                   |                                                                                                    
ACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTA  >  NC_000913/3888113‑3888311

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: ‑