Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,888,219 T→A 16.2% intergenic (+27/‑216) 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,2190TA16.2% 492.4 / 81.8 196intergenic (+27/‑216)mnmE/tnaCtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase/tryptophanase leader peptide
Reads supporting (aligned to +/- strand):  ref base T (80/84);  new base A (15/17);  total (95/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-01

Only 100 of 201 total aligned reads displayed.

GAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAG  >  NC_000913/3888121‑3888318
                                                                                                  |                                                                                                    
gAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTg                                                                                                     <  2:1639656/100‑1 (MQ=255)
    tCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAg                                                                                                 >  1:585112/1‑100 (MQ=255)
     cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGa                                                                                               >  1:3743364/1‑100 (MQ=255)
     cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGa                                                                                               >  2:3847167/1‑100 (MQ=255)
       ccGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATg                                                                                             >  2:101026/1‑100 (MQ=255)
       ccGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATg                                                                                             >  1:663155/1‑100 (MQ=255)
        cGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGc                                                                                            >  1:1769963/1‑100 (MQ=255)
          gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCgg                                                                                          <  1:2233479/100‑1 (MQ=38)
           ggAATTTACTTCAGACGACCTGCTGGGTCGTATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgc                                                                                         <  1:396282/100‑1 (MQ=37)
           ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgc                                                                                         <  1:3506462/100‑1 (MQ=38)
           ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgc                                                                                         <  2:2186982/100‑1 (MQ=38)
             aaTTTACTTCAGACGACCTGCTGGGGCGGATTTGCTCCAGCTTCTGTATGGGTAAGTAACCGCGCTTCCGAAGCCGCATCTGACAGTCAGAATGCGgctt                                                                                       <  2:234624/100‑1 (MQ=25)
             aaTTCACTTCAGACGATCTGCTGGGGCGAATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAATCCGCATCTGACAGTCAGAATGCGgc                                                                                         <  2:539371/98‑1 (MQ=37)
              aTTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttc                                                                                      <  2:1504206/100‑1 (MQ=255)
               tttACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGTTTCTGTATTGGTAAGTAACCGCGCTGACGAAGCCGCATTCTGACTGTCAGAGGCGgcttc                                                                                      <  1:1695846/99‑1 (MQ=255)
               tttACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcg                                                                                     <  2:1902161/100‑1 (MQ=255)
               tttACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgcttcg                                                                                     <  1:1733407/100‑1 (MQ=38)
                 tACTTCAGACGACCTGCTGGGGCGGATTTTTTCCATCTTTTTTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgct                                                                                   <  1:458201/100‑1 (MQ=255)
                   cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCGGATGCGgcttcgcttc                                                                                 >  1:3363077/1‑100 (MQ=255)
                   cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgcttc                                                                                 >  2:1803406/1‑100 (MQ=255)
                   cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATACGgcttcgcttc                                                                                 >  1:1940006/1‑100 (MQ=255)
                    ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCa                                                                                <  2:619356/100‑1 (MQ=255)
                    ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCa                                                                                <  1:2115837/100‑1 (MQ=255)
                       aGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCAtt                                                                              <  1:759493/99‑1 (MQ=255)
                       aGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTg                                                                             <  1:2843723/100‑1 (MQ=255)
                          cgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTa                                                                          <  2:1114733/100‑1 (MQ=255)
                           gacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTAc                                                                         <  2:1219294/100‑1 (MQ=255)
                            acCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTAcc                                                                        <  2:1936544/100‑1 (MQ=255)
                             ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCa                                                                       >  1:2283822/1‑100 (MQ=255)
                             ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCa                                                                       >  1:1898834/1‑100 (MQ=255)
                             ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCa                                                                       >  2:1786050/1‑100 (MQ=255)
                             ccTGCTGGGGCGGATTTTCTCCAGCTGCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCa                                                                       >  1:2387595/1‑100 (MQ=255)
                              ctgctgGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTCCCAc                                                                      <  2:2955098/100‑1 (MQ=255)
                                gctgGGGCGGATTTTCTCCAGCTTCTGTATAGGAAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATAGTTACCACTc                                                                    >  2:2564974/1‑100 (MQ=255)
                                   ggggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTg                                                                 >  2:1942290/1‑100 (MQ=255)
                                    gggCGGATTTTCTCCAGTTTTTGTATTGGTAAGTAACCGCGTTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGt                                                                <  1:2618380/100‑1 (MQ=255)
                                    gggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGt                                                                <  1:2922435/100‑1 (MQ=255)
                                     ggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGtt                                                               >  2:1599202/1‑100 (MQ=255)
                                       cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                             >  2:482638/1‑100 (MQ=255)
                                       cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                             <  2:3265186/100‑1 (MQ=255)
                                       cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                             >  2:2773764/1‑100 (MQ=255)
                                         gATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTc                                                           <  2:825071/100‑1 (MQ=255)
                                          aTTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTcc                                                          >  1:3911156/1‑100 (MQ=255)
                                          aTTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTcc                                                          >  1:71244/1‑100 (MQ=255)
                                          aTTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTcc                                                          <  2:361362/100‑1 (MQ=255)
                                             ttCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCa                                                       <  2:2155098/100‑1 (MQ=255)
                                               ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCaag                                                     >  1:1843645/1‑99 (MQ=255)
                                               ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                     >  2:3830370/1‑100 (MQ=255)
                                               ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                     <  1:2809632/100‑1 (MQ=255)
                                                tcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAcc                                                    <  2:835839/100‑1 (MQ=255)
                                                  catcTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCTCTCCTGTTATTCCTCAAc                                                     <  2:1332129/94‑1 (MQ=255)
                                                  cAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCt                                                  >  2:3068183/1‑100 (MQ=255)
                                                       tCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCtttttt                                             <  1:1268367/100‑1 (MQ=255)
                                                         tGTATTGGTAAGTAACCGCGCTTACGAAGCTGCATTTTGACTGTCAGATGCGGGTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                           <  1:563963/100‑1 (MQ=255)
                                                         tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                           <  2:1531628/100‑1 (MQ=255)
                                                         tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTAAGATGCGGATTCGCTTCATAGTTAACAATCCTGTTATTCCTCAACCCTTTTTTTa                                           >  1:3679861/1‑100 (MQ=255)
                                                            aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCTGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                        <  2:2742269/100‑1 (MQ=255)
                                                            aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                        <  2:3496974/100‑1 (MQ=255)
                                                            aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                        <  2:2935386/100‑1 (MQ=255)
                                                            aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                        >  1:1577788/1‑100 (MQ=255)
                                                            aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                        >  1:3846792/1‑100 (MQ=255)
                                                             ttGGTAAGTAACCGCGGTTACGAAGCCTCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCATCCCTTTTTTTAAAc                                        <  2:3361580/99‑1 (MQ=255)
                                                             ttGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACa                                       >  1:3317040/1‑100 (MQ=255)
                                                             ttGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACa                                       >  2:363415/1‑100 (MQ=255)
                                                               ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACAtt                                     >  2:1738988/1‑100 (MQ=255)
                                                               ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACAtt                                     >  1:282009/1‑100 (MQ=255)
                                                                 taagtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTaa                                   >  2:704327/1‑100 (MQ=255)
                                                                   agtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTaaaa                                 >  2:2442872/1‑100 (MQ=255)
                                                                   agtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACAACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTaaaa                                 >  1:1165266/1‑100 (MQ=255)
                                                                    gtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAAt                                <  2:2338883/100‑1 (MQ=255)
                                                                     taaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTAAACCCTTTTTTTAAACATTAAAAtt                               >  1:1058170/1‑100 (MQ=255)
                                                                     taaCCGCGCTTACGAAGCCACATTCTGACAGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAAtt                               >  2:907382/1‑100 (MQ=255)
                                                                      aaCCGTGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTc                              <  1:856117/100‑1 (MQ=255)
                                                                              ttATGATGCCGCATTCTGACTGTCAAATGCGGCTTCGCTTGATTGTTACCTCTCCTGTTATTAATTAACCCTTTTTTTAAAGATGAATATACTTACg                         <  2:3736567/97‑1 (MQ=255)
                                                                              aagcGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTaa                      >  1:2444461/4‑100 (MQ=37)
                                                                               tACGAAGCCGCATTCTGACTGTCAGATGCGTCTTTGCTTCATTGTTTCCATTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt                     <  1:1219044/100‑1 (MQ=255)
                                                                               tACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt                     >  2:609427/1‑100 (MQ=255)
                                                                                 cGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAttt                   >  1:3089563/1‑100 (MQ=37)
                                                                                  gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta                  >  2:2203782/1‑100 (MQ=38)
                                                                                  gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta                  >  1:1467840/1‑100 (MQ=38)
                                                                                  gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCATTTTTTTAAACATTAAAATTCTTACGTAATTta                  >  1:1250690/1‑100 (MQ=37)
                                                                                  gAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTACTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta                  >  1:1089561/1‑100 (MQ=37)
                                                                                     gCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATaa               >  1:296646/1‑100 (MQ=255)
                                                                                     gCCGCATCTGACAGTCAGAATGCGGCTTCGCTTAATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACTTTAAAATTCTTACGTAATTTATaa               >  2:3744938/1‑100 (MQ=37)
                                                                                      ccGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTTCGTAATTTATAAt              <  2:2422977/100‑1 (MQ=255)
                                                                                        gtaTTCTGACTGTCAGATGCTGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAATCCATTTTTTAAACATTAAAATTCTTACGTAATTTATAATCt            <  2:3810075/98‑1 (MQ=255)
                                                                                          aTTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCttt          <  2:1118515/100‑1 (MQ=255)
                                                                                          aTTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCttt          <  1:1977625/100‑1 (MQ=255)
                                                                                           ttCTGACTGTCAGATGCGGCTTTGCTTCATTGTTACTTCTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTa         <  1:2505299/100‑1 (MQ=255)
                                                                                             cTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaa       <  1:534262/100‑1 (MQ=255)
                                                                                             cTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaa       >  1:2887818/1‑100 (MQ=255)
                                                                                             cTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaa       >  2:181877/1‑100 (MQ=255)
                                                                                             cTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaa       <  1:3337544/100‑1 (MQ=255)
                                                                                               gACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaaaa     >  2:3150822/1‑100 (MQ=255)
                                                                                                 cTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaaaaaa   <  2:1866467/100‑1 (MQ=255)
                                                                                                 cTGTCAGATGCGACTTCGCTTCATTGTTACCACACCTATTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTAGGTAATTTATAATCTTTaaaaaaa   >  2:2588184/1‑100 (MQ=255)
                                                                                                  tGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCATCCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAg  <  1:2153585/100‑1 (MQ=255)
                                                                                                  tGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAg  <  1:3520050/100‑1 (MQ=255)
                                                                                                  tGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAg  <  2:2913161/100‑1 (MQ=255)
                                                                                                  tGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAg  <  1:1325330/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
GAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAG  >  NC_000913/3888121‑3888318

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: