Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,888,225:1 +A 15.3% intergenic (+33/‑210) mnmE → / → tnaC tRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase/tryptophanase leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,888,2241.A15.3% 644.8 / 125.3 196intergenic (+32/‑211)mnmE/tnaCtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase/tryptophanase leader peptide
Reads supporting (aligned to +/- strand):  ref base . (77/89);  new base A (14/16);  total (91/105)
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.20e-01

Only 100 of 201 total aligned reads displayed.

TCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTT  >  NC_000913/3888125‑3888323
                                                                                                    |                                                                                                    
 cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGa                                                                                                     >  2:3847167/1‑100 (MQ=255)
   ccGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATg                                                                                                   >  2:101026/1‑100 (MQ=255)
    cGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGc                                                                                                  >  1:1769963/1‑100 (MQ=255)
      gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCgg                                                                                                <  2:3616880/100‑1 (MQ=255)
      gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCgg                                                                                                <  1:3068317/100‑1 (MQ=38)
      gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCgg                                                                                                <  1:2233479/100‑1 (MQ=38)
       ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgc                                                                                               <  2:2186982/100‑1 (MQ=38)
        gAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCg                                                                                                 <  2:3187545/97‑1 (MQ=255)
         aaTTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgctt                                                                                             <  1:664733/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)
             tACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgct                                                                                         <  1:2931536/100‑1 (MQ=255)
               cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgcttc                                                                                       >  2:1803406/1‑100 (MQ=255)
               cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgcttc                                                                                       >  1:3576141/1‑100 (MQ=255)
               cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATACGgcttcgcttc                                                                                       >  1:1940006/1‑100 (MQ=255)
                ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCa                                                                                      <  1:2115837/100‑1 (MQ=255)
                ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCa                                                                                      <  2:619356/100‑1 (MQ=255)
                ttCAGACGACCTGCTGGGGCGGAATTTCTCCAGCTTCTGTATTGGTAAGTAACCGAGCTTACGAAGCCGCATGCTGACTGTCAGATGCGGCTTCGCTTCa                                                                                      >  2:3141351/1‑100 (MQ=255)
                  cAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCAtt                                                                                    >  1:3241626/1‑100 (MQ=255)
                    gacgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGt                                                                                  <  1:3353644/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)
                         ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTAcc                                                                              <  1:2901863/99‑1 (MQ=255)
                         ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCa                                                                             >  1:1898834/1‑100 (MQ=255)
                         ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCa                                                                             >  2:1786050/1‑100 (MQ=255)
                         ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCa                                                                             >  1:2283822/1‑100 (MQ=255)
                         ccTGCTGGGGCGGATTTTCTCCAGCTGCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCa                                                                             >  1:2387595/1‑100 (MQ=255)
                            gctgGGGCGGATTTTCTCCAGCTTCTGTATAGGAAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATAGTTACCACTc                                                                          >  2:2564974/1‑100 (MQ=255)
                               ggggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTg                                                                       >  2:1942290/1‑100 (MQ=255)
                                gggCGGATTTTCTCCAGTTTTTGTATTGGTAAGTAACCGCGTTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGt                                                                      <  1:2618380/100‑1 (MQ=255)
                                  gCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTa                                                                    >  2:336056/1‑100 (MQ=255)
                                   cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                                   >  2:3812627/1‑100 (MQ=255)
                                   cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                                   >  2:482638/1‑100 (MQ=255)
                                   cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                                   >  1:623541/1‑100 (MQ=255)
                                   cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                                   >  2:2773764/1‑100 (MQ=255)
                                     gATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTc                                                                 <  2:825071/100‑1 (MQ=255)
                                         ttCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCa                                                             <  2:2155098/100‑1 (MQ=255)
                                          tctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCaa                                                            <  2:3376153/100‑1 (MQ=255)
                                           ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCaag                                                           >  1:1843645/1‑99 (MQ=255)
                                           ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                           <  1:2809632/100‑1 (MQ=255)
                                            tcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAcc                                                          >  1:449801/1‑100 (MQ=255)
                                            tcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAcc                                                          <  2:835839/100‑1 (MQ=255)
                                              catcTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCTCTCCTGTTATTCCTCAAc                                                           <  2:1332129/94‑1 (MQ=255)
                                                 cTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCttt                                                      <  1:3441037/99‑1 (MQ=255)
                                                   tCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCtttttt                                                  <  1:1268367/100‑1 (MQ=255)
                                                   tCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCtttttt                                                  <  2:3741883/100‑1 (MQ=255)
                                                     tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTAAGATGCGGATTCGCTTCATAGTTAACAATCCTGTTATTCCTCAACCCTTTTTTTa                                                >  1:3679861/1‑100 (MQ=255)
                                                        aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCTGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                             <  2:2742269/100‑1 (MQ=255)
                                                        aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                             >  1:1577788/1‑100 (MQ=255)
                                                          tGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACAt                                           <  2:3550521/100‑1 (MQ=255)
                                                           ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACAtt                                          >  2:1738988/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)
                                                                      gcgcTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTAc                               <  2:2943804/100‑1 (MQ=255)
                                                                       cgcTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACg                              <  1:545610/100‑1 (MQ=255)
                                                                         cTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTa                            >  1:754826/1‑100 (MQ=255)
                                                                          aagcGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTaa                           >  1:2444461/4‑100 (MQ=37)
                                                                           tACGAAGCCGCATTCTGACTGTCAGATGCGTCTTTGCTTCATTGTTTCCATTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt                          <  1:1219044/100‑1 (MQ=255)
                                                                           tACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt                          >  1:3755482/1‑100 (MQ=255)
                                                                           tACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAt                          >  2:609427/1‑100 (MQ=255)
                                                                           tACGAAGCCGCATTCTGACTGTAAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTATTTTAAACATTAAAATTCTTACGTa                            >  2:3008303/1‑99 (MQ=255)
                                                                            aCGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTTCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAtt                         <  1:52329/100‑1 (MQ=255)
                                                                             cGAAGCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAttt                        >  1:569808/1‑100 (MQ=37)
                                                                              gAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta                       >  1:800505/1‑100 (MQ=255)
                                                                              gAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTta                       <  1:2884521/100‑1 (MQ=255)
                                                                              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                    >  2:3636169/1‑100 (MQ=255)
                                                                                 gCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATaa                    <  1:3509999/100‑1 (MQ=255)
                                                                                 gCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATaa                    >  1:3333425/1‑100 (MQ=38)
                                                                                  ccGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTTCGTAATTTATAAt                   <  2:2422977/100‑1 (MQ=255)
                                                                                  ccGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATaa                    >  1:3877057/1‑99 (MQ=255)
                                                                                  ccGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAAt                   >  2:1607871/1‑100 (MQ=255)
                                                                                   cGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATc                  >  1:2831057/1‑100 (MQ=255)
                                                                                      aTTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCttt               <  2:3755443/100‑1 (MQ=255)
                                                                                      aTTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCttt               <  1:1977625/100‑1 (MQ=255)
                                                                                      aTTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCttt               <  2:1118515/100‑1 (MQ=255)
                                                                                       ttCTGACTGTCAGATGCGGCTTTGCTTCATTGTTACTTCTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTa              <  1:2505299/100‑1 (MQ=255)
                                                                                         cTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaa            >  2:1555363/1‑100 (MQ=255)
                                                                                         cTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaa            >  2:181877/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       <  2:2913161/100‑1 (MQ=255)
                                                                                              tGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAg       <  1:1325330/100‑1 (MQ=255)
                                                                                               gTCAGATGCGGCTTCGCTTCATTGTTACGAATCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAg       <  2:2313643/99‑1 (MQ=255)
                                                                                                tCAGATGCGGCTTCGATTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATATTTAAAAAAAGCa     <  2:383232/100‑1 (MQ=255)
                                                                                                   gATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCAttt  >  1:3912599/1‑100 (MQ=255)
                                                                                                    |                                                                                                    
TCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTT  >  NC_000913/3888125‑3888323

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: