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

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

Only 100 of 198 total aligned reads displayed.

ACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTA  >  NC_000913/3888113‑3888311
                                                                                                   |                                                                                                     
aCTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCAt                                                                                                       >  1:361368/1‑100 (MQ=255)
   tAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCt                                                                                                    >  1:3550560/1‑100 (MQ=255)
    aaGCGAAATCACCGGGGAATTTACTAAAGAAGACCTGCTGGGGCGGATTTTATCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTg                                                                                                   >  2:92430/1‑100 (MQ=255)
      gagAAATCACCGGGGTATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGAc                                                                                                 <  1:3867823/98‑1 (MQ=255)
      gCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGATAAGTAACCGAGCTTACGAAGCCGCATTCTGAc                                                                                                 >  1:248974/1‑100 (MQ=255)
        gAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTg                                                                                               <  2:1639656/100‑1 (MQ=255)
        gAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTg                                                                                               <  2:1526750/100‑1 (MQ=255)
         aaaTCACCGGTGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGt                                                                                              <  1:3736740/100‑1 (MQ=255)
         aaaTCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGt                                                                                              >  2:2931510/1‑100 (MQ=255)
             cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑a                                                                                         >  2:3847167/1‑100 (MQ=255)
             cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑a                                                                                         >  2:2316204/1‑100 (MQ=255)
              aCCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTAAG‑At                                                                                        >  1:1728180/1‑100 (MQ=255)
               ccGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑At                                                                                        >  2:782241/1‑99 (MQ=255)
               ccGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATg                                                                                       >  2:101026/1‑100 (MQ=255)
                 ggggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCg                                                                                     <  1:1951896/100‑1 (MQ=255)
                  gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCgg                                                                                    <  1:3068317/100‑1 (MQ=38)
                  gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCgg                                                                                    <  1:2233479/100‑1 (MQ=38)
                   ggAATTTACTTCAGACGACCTGCTGGGTCGTATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgc                                                                                   <  1:396282/100‑1 (MQ=37)
                   ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgc                                                                                   <  2:2186982/100‑1 (MQ=38)
                   ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgc                                                                                   <  1:3506462/100‑1 (MQ=38)
                    gAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgct                                                                                  >  1:1840048/1‑100 (MQ=255)
                    gAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGCCAGAATGCGgct                                                                                  <  1:2582427/100‑1 (MQ=38)
                    gAAGTTACTTCAGACCACCTGCTGGGGCGGATTTTCTCCGGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgc                                                                                   <  2:1115530/99‑1 (MQ=37)
                     aaTTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgctt                                                                                 <  1:664733/100‑1 (MQ=38)
                     aaTTCACTTCAGACGATCTGCTGGGGCGAATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAATCCGCA‑TCTGACAGTCAGAATGCGgc                                                                                   <  2:539371/98‑1 (MQ=37)
                       tttACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcg                                                                               <  2:1902161/100‑1 (MQ=255)
                        ttACTTCAGATGACCTGCTGGTGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgcttcgt                                                                              <  1:1867085/100‑2 (MQ=25)
                         tACTTCAGACGACCTGCTGGGGCGGATTTTTTCCATCTTTTTTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgct                                                                             <  1:458201/100‑1 (MQ=255)
                         tACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgct                                                                             <  1:2931536/100‑1 (MQ=255)
                          aCTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTTTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgctt                                                                            <  1:2152256/100‑1 (MQ=255)
                           cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCGG‑ATGCGgcttcgcttc                                                                           >  1:3363077/1‑100 (MQ=255)
                           cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgcttc                                                                           >  1:1936560/1‑100 (MQ=255)
                           cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgcttc                                                                           >  2:1803406/1‑100 (MQ=255)
                           cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTATCAG‑ATGCGgcttcgcttc                                                                           >  2:2452407/1‑100 (MQ=255)
                            ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCa                                                                          <  1:2115817/100‑1 (MQ=255)
                            ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCa                                                                          <  2:619356/100‑1 (MQ=255)
                            ttCAGACGACCTGCTGGGGCGGAATTTCTCCAGCTTCTGTATTGGTAAGTAACCGAGCTTACGAAGCCGCATGCTGACTGTCAG‑ATGCGGCTTCGCTTCa                                                                          >  2:3141351/1‑100 (MQ=255)
                              cAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCAtt                                                                        >  1:3241626/1‑100 (MQ=255)
                               aGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTg                                                                       <  1:2843723/100‑1 (MQ=255)
                                  cgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTa                                                                    <  2:1114733/100‑1 (MQ=255)
                                    acCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTAcc                                                                  <  2:1936544/100‑1 (MQ=255)
                                     ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTAcc                                                                  <  1:2901863/99‑1 (MQ=255)
                                     ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCa                                                                 >  2:1786050/1‑100 (MQ=255)
                                     ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCa                                                                 >  1:2283822/1‑100 (MQ=255)
                                     ccTGCTGGGGCGGATTTTCTCCAGCTGCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCa                                                                 >  1:2387595/1‑100 (MQ=255)
                                          tgGGGCGGATTTTCTCCAGCTTCCGAATTGGTACGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑GTGCGGCTCCGCTTCATTGTTACCACTCCt                                                            >  2:2356186/1‑100 (MQ=255)
                                           ggggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTg                                                           >  2:1942290/1‑100 (MQ=255)
                                             ggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGtt                                                         >  2:1599202/1‑100 (MQ=255)
                                             ggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGtt                                                         <  2:1167362/100‑1 (MQ=255)
                                              gCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTa                                                        >  2:2430929/1‑100 (MQ=255)
                                               cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                       >  2:482638/1‑100 (MQ=255)
                                               cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                       >  1:623541/1‑100 (MQ=255)
                                               cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                       >  2:2773764/1‑100 (MQ=255)
                                               cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                       <  2:3265186/100‑1 (MQ=255)
                                                  aTTTTCTCCAGCTTCTGTATTGGTAAGTAATCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTcc                                                    <  1:1758217/100‑1 (MQ=255)
                                                  aTTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTcc                                                    >  1:3911156/1‑100 (MQ=255)
                                                  aTTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTcc                                                    <  2:361362/100‑1 (MQ=255)
                                                     ttCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCa                                                 <  2:2155098/100‑1 (MQ=255)
                                                     ctctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCa                                                 <  2:1843630/99‑1 (MQ=255)
                                                       ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                               <  1:2809632/100‑1 (MQ=255)
                                                       ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                               >  2:3830370/1‑100 (MQ=255)
                                                          catcTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCTCTCCTGTTATTCCTCAAc                                               <  2:1332129/94‑1 (MQ=255)
                                                                 tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTa                                    >  2:853907/1‑100 (MQ=255)
                                                                    aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCTG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAAc                                 <  2:2742269/100‑1 (MQ=255)
                                                                    aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAAc                                 <  2:3496974/100‑1 (MQ=255)
                                                                    aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAAc                                 >  1:1577788/1‑100 (MQ=255)
                                                                     ttGGTAAGTAACCGCGGTTACGAAGCCTCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCATCCCTTT‑TTTTAAAc                                 <  2:3361580/99‑1 (MQ=255)
                                                                     ttGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACa                                >  2:363415/1‑100 (MQ=255)
                                                                       ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACAtt                              >  1:282009/1‑100 (MQ=255)
                                                                       ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACAtt                              >  2:1738988/1‑100 (MQ=255)
                                                                        gtaagtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTa                             >  1:130816/1‑100 (MQ=255)
                                                                            gtaaCCGCGCTTACGAAGCCGGATTCTGACTGTCAG‑ATGCGGCTTCGCTTTATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAAt                         <  2:2339913/100‑1 (MQ=255)
                                                                             taaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTAAACCCTTT‑TTTTAAACATTAAAAtt                        >  1:1058170/1‑100 (MQ=255)
                                                                             taaCCGCGCTTACGAAGCCACATTCTGACAGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAAtt                        >  2:907382/1‑100 (MQ=255)
                                                                              aaCCGTGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTc                       <  1:856117/100‑1 (MQ=255)
                                                                                ccGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTTAAATTCt                      >  2:1232056/1‑99 (MQ=255)
                                                                                 cgcgCTTACGAAGCCGCATTCTGACTGTCAG‑CTGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTa                    >  1:1199902/1‑100 (MQ=255)
                                                                                  gcgcTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTa                    <  2:827291/99‑1 (MQ=255)
                                                                                  gcgcTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTAc                   <  2:2943804/100‑1 (MQ=255)
                                                                                   cgcTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACg                  <  1:545610/100‑1 (MQ=255)
                                                                                     cTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTa                >  1:754826/1‑100 (MQ=255)
                                                                                      ttATGATGCCGCATTCTGACTGTCAA‑ATGCGGCTTCGCTTGATTGTTACCTCTCCTGTTATTAATTAACCCTTT‑TTTTAAAGATGAATATACTTACg                  <  2:3736567/97‑1 (MQ=255)
                                                                                      aagcGAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTaa               >  1:2444461/4‑100 (MQ=37)
                                                                                       tACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAAt              >  1:3755482/1‑100 (MQ=255)
                                                                                       tACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAAt              >  2:609427/1‑100 (MQ=255)
                                                                                       tACGAAGCCGCATTCTGACTGTAAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTATTTTAAACATTAAAATTCTTACGTa                >  2:3008303/1‑99 (MQ=255)
                                                                                       agcGAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAAt              >  1:1223715/3‑100 (MQ=37)
                                                                                         cGAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAAttt            >  1:569808/1‑100 (MQ=37)
                                                                                          gAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTta           <  1:2730811/100‑1 (MQ=255)
                                                                                          gAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTta           >  1:800505/1‑100 (MQ=255)
                                                                                          gAAGCCGCATCCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCATTC‑TTTTAAACATTAAAATTCTTACGTAATTta           <  2:2569236/100‑1 (MQ=255)
                                                                                          gAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTTCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTta           >  1:1058335/1‑100 (MQ=37)
                                                                                          gAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTta           >  1:1467840/1‑100 (MQ=38)
                                                                                          gAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTta           >  1:1262071/1‑100 (MQ=38)
                                                                                          gAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTta           >  2:2203782/1‑100 (MQ=38)
                                                                                             gCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATaa        >  1:3333425/1‑100 (MQ=38)
                                                                                             gCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTAATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACTTTAAAATTCTTACGTAATTTATaa        >  2:3744938/1‑100 (MQ=37)
                                                                                               cacATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAAt       >  2:1399935/3‑99 (MQ=255)
                                                                                                gtaTTCTGACTGTCAG‑ATGCTGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAATCCATT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCt     <  2:3810075/98‑1 (MQ=255)
                                                                                                   ttCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTa  <  1:1098807/100‑1 (MQ=255)
                                                                                                   |                                                                                                     
ACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTA  >  NC_000913/3888113‑3888311

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