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.

GAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAG  >  NC_000913/3888121‑3888318
                                                                                                  |                                                                                                    
gAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTg                                                                                                     <  2:1639656/100‑1 (MQ=255)
    tCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAg                                                                                                 >  1:585112/1‑100 (MQ=255)
     cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑a                                                                                               >  1:3743364/1‑100 (MQ=255)
     cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑a                                                                                               >  2:3847167/1‑100 (MQ=255)
       ccGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATg                                                                                             >  2:101026/1‑100 (MQ=255)
       ccGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATg                                                                                             >  1:663155/1‑100 (MQ=255)
        cGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGc                                                                                            >  1:1769963/1‑100 (MQ=255)
          gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCgg                                                                                          <  1:2233479/100‑1 (MQ=38)
           ggAATTTACTTCAGACGACCTGCTGGGTCGTATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgc                                                                                         <  1:396282/100‑1 (MQ=37)
           ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgc                                                                                         <  1:3506462/100‑1 (MQ=38)
           ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgc                                                                                         <  2:2186982/100‑1 (MQ=38)
             aaTTTACTTCAGACGACCTGCTGGGGCGGATTTGCTCCAGCTTCTGTATGGGTAAGTAACCGCGCTTCCGAAGCCGCA‑TCTGACAGTCAGAATGCGgctt                                                                                       <  2:234624/100‑1 (MQ=25)
             aaTTCACTTCAGACGATCTGCTGGGGCGAATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAATCCGCA‑TCTGACAGTCAGAATGCGgc                                                                                         <  2:539371/98‑1 (MQ=37)
              aTTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttc                                                                                      <  2:1504206/100‑1 (MQ=255)
               tttACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGTTTCTGTATTGGTAAGTAACCGCGCTGACGAAGCCGCATTCTGACTGTCAG‑AGGCGgcttc                                                                                      <  1:1695846/99‑1 (MQ=255)
               tttACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcg                                                                                     <  2:1902161/100‑1 (MQ=255)
               tttACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgcttcg                                                                                     <  1:1733407/100‑1 (MQ=38)
                 tACTTCAGACGACCTGCTGGGGCGGATTTTTTCCATCTTTTTTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgct                                                                                   <  1:458201/100‑1 (MQ=255)
                   cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCGG‑ATGCGgcttcgcttc                                                                                 >  1:3363077/1‑100 (MQ=255)
                   cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgcttc                                                                                 >  2:1803406/1‑100 (MQ=255)
                   cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATACGgcttcgcttc                                                                                 >  1:1940006/1‑100 (MQ=255)
                    ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCa                                                                                <  2:619356/100‑1 (MQ=255)
                    ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCa                                                                                <  1:2115837/100‑1 (MQ=255)
                       aGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCAtt                                                                              <  1:759493/99‑1 (MQ=255)
                       aGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTg                                                                             <  1:2843723/100‑1 (MQ=255)
                          cgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTa                                                                          <  2:1114733/100‑1 (MQ=255)
                           gacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTAc                                                                         <  2:1219294/100‑1 (MQ=255)
                            acCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTAcc                                                                        <  2:1936544/100‑1 (MQ=255)
                             ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCa                                                                       >  1:2283822/1‑100 (MQ=255)
                             ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCa                                                                       >  1:1898834/1‑100 (MQ=255)
                             ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCa                                                                       >  2:1786050/1‑100 (MQ=255)
                             ccTGCTGGGGCGGATTTTCTCCAGCTGCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCa                                                                       >  1:2387595/1‑100 (MQ=255)
                              ctgctgGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTCCCAc                                                                      <  2:2955098/100‑1 (MQ=255)
                                gctgGGGCGGATTTTCTCCAGCTTCTGTATAGGAAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATAGTTACCACTc                                                                    >  2:2564974/1‑100 (MQ=255)
                                   ggggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTg                                                                 >  2:1942290/1‑100 (MQ=255)
                                    gggCGGATTTTCTCCAGTTTTTGTATTGGTAAGTAACCGCGTTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGt                                                                <  1:2618380/100‑1 (MQ=255)
                                    gggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGt                                                                <  1:2922435/100‑1 (MQ=255)
                                     ggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGtt                                                               >  2:1599202/1‑100 (MQ=255)
                                       cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                             >  2:482638/1‑100 (MQ=255)
                                       cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                             <  2:3265186/100‑1 (MQ=255)
                                       cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                             >  2:2773764/1‑100 (MQ=255)
                                         gATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTc                                                           <  2:825071/100‑1 (MQ=255)
                                          aTTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTcc                                                          >  1:3911156/1‑100 (MQ=255)
                                          aTTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTcc                                                          >  1:71244/1‑100 (MQ=255)
                                          aTTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTcc                                                          <  2:361362/100‑1 (MQ=255)
                                             ttCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCa                                                       <  2:2155098/100‑1 (MQ=255)
                                               ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCaag                                                     >  1:1843645/1‑99 (MQ=255)
                                               ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                     >  2:3830370/1‑100 (MQ=255)
                                               ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                     <  1:2809632/100‑1 (MQ=255)
                                                tcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAcc                                                    <  2:835839/100‑1 (MQ=255)
                                                  catcTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCTCTCCTGTTATTCCTCAAc                                                     <  2:1332129/94‑1 (MQ=255)
                                                  cAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCt                                                  >  2:3068183/1‑100 (MQ=255)
                                                       tCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCtttttt                                             <  1:1268367/100‑1 (MQ=255)
                                                         tGTATTGGTAAGTAACCGCGCTTACGAAGCTGCATTTTGACTGTCAG‑ATGCGGGTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                           <  1:563963/100‑1 (MQ=255)
                                                         tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                           <  2:1531628/100‑1 (MQ=255)
                                                         tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTAAG‑ATGCGGATTCGCTTCATAGTTAACAATCCTGTTATTCCTCAACCCTTTTTTTa                                           >  1:3679861/1‑100 (MQ=255)
                                                            aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCTG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                        <  2:2742269/100‑1 (MQ=255)
                                                            aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                        <  2:3496974/100‑1 (MQ=255)
                                                            aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                        <  2:2935386/100‑1 (MQ=255)
                                                            aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                        >  1:1577788/1‑100 (MQ=255)
                                                            aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                        >  1:3846792/1‑100 (MQ=255)
                                                             ttGGTAAGTAACCGCGGTTACGAAGCCTCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCATCCCTTTTTTTAAAc                                        <  2:3361580/99‑1 (MQ=255)
                                                             ttGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACa                                       >  1:3317040/1‑100 (MQ=255)
                                                             ttGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACa                                       >  2:363415/1‑100 (MQ=255)
                                                               ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACAtt                                     >  2:1738988/1‑100 (MQ=255)
                                                               ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACAtt                                     >  1:282009/1‑100 (MQ=255)
                                                                 taagtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTaa                                   >  2:704327/1‑100 (MQ=255)
                                                                   agtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTaaaa                                 >  2:2442872/1‑100 (MQ=255)
                                                                   agtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACAACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTaaaa                                 >  1:1165266/1‑100 (MQ=255)
                                                                    gtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAAt                                <  2:2338883/100‑1 (MQ=255)
                                                                     taaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTAAACCCTTTTTTTAAACATTAAAAtt                               >  1:1058170/1‑100 (MQ=255)
                                                                     taaCCGCGCTTACGAAGCCACATTCTGACAGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAAtt                               >  2:907382/1‑100 (MQ=255)
                                                                      aaCCGTGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTc                              <  1:856117/100‑1 (MQ=255)
                                                                              ttATGATGCCGCATTCTGACTGTCAA‑ATGCGGCTTCGCTTGATTGTTACCTCTCCTGTTATTAATTAACCCTTTTTTTAAAGATGAATATACTTACg                         <  2:3736567/97‑1 (MQ=255)
                                                                              aagcGAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTaa                      >  1:2444461/4‑100 (MQ=37)
                                                                               tACGAAGCCGCATTCTGACTGTCAG‑ATGCGTCTTTGCTTCATTGTTTCCATTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt                     <  1:1219044/100‑1 (MQ=255)
                                                                               tACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt                     >  2:609427/1‑100 (MQ=255)
                                                                                 cGAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAttt                   >  1:3089563/1‑100 (MQ=37)
                                                                                  gAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta                  >  2:2203782/1‑100 (MQ=38)
                                                                                  gAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta                  >  1:1467840/1‑100 (MQ=38)
                                                                                  gAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCATTTTTTTAAACATTAAAATTCTTACGTAATTta                  >  1:1250690/1‑100 (MQ=37)
                                                                                  gAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTACTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta                  >  1:1089561/1‑100 (MQ=37)
                                                                                     gCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATaa               >  1:296646/1‑100 (MQ=255)
                                                                                     gCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTAATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACTTTAAAATTCTTACGTAATTTATaa               >  2:3744938/1‑100 (MQ=37)
                                                                                      ccGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTTCGTAATTTATAAt              <  2:2422977/100‑1 (MQ=255)
                                                                                        gtaTTCTGACTGTCAG‑ATGCTGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAATCCATTTTTTAAACATTAAAATTCTTACGTAATTTATAATCt            <  2:3810075/98‑1 (MQ=255)
                                                                                          aTTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCttt          <  2:1118515/100‑1 (MQ=255)
                                                                                          aTTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCttt          <  1:1977625/100‑1 (MQ=255)
                                                                                           ttCTGACTGTCAG‑ATGCGGCTTTGCTTCATTGTTACTTCTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTa         <  1:2505299/100‑1 (MQ=255)
                                                                                             cTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaa       <  1:534262/100‑1 (MQ=255)
                                                                                             cTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaa       >  1:2887818/1‑100 (MQ=255)
                                                                                             cTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaa       >  2:181877/1‑100 (MQ=255)
                                                                                             cTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaa       <  1:3337544/100‑1 (MQ=255)
                                                                                               gACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaaaa     >  2:3150822/1‑100 (MQ=255)
                                                                                                 cTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaaaaaa   <  2:1866467/100‑1 (MQ=255)
                                                                                                 cTGTCAG‑ATGCGACTTCGCTTCATTGTTACCACACCTATTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTAGGTAATTTATAATCTTTaaaaaaa   >  2:2588184/1‑100 (MQ=255)
                                                                                                  tGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCATCCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAg  <  1:2153585/100‑1 (MQ=255)
                                                                                                  tGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAg  <  1:3520050/100‑1 (MQ=255)
                                                                                                  tGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAg  <  2:2913161/100‑1 (MQ=255)
                                                                                                  tGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAg  <  1:1325330/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
GAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAG  >  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: ‑