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.

TCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTT  >  NC_000913/3888125‑3888323
                                                                                                    |                                                                                                    
 cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑a                                                                                                     >  2:3847167/1‑100 (MQ=255)
   ccGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATg                                                                                                   >  2:101026/1‑100 (MQ=255)
    cGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGc                                                                                                  >  1:1769963/1‑100 (MQ=255)
      gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCgg                                                                                                <  2:3616880/100‑1 (MQ=255)
      gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCgg                                                                                                <  1:3068317/100‑1 (MQ=38)
      gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCgg                                                                                                <  1:2233479/100‑1 (MQ=38)
       ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgc                                                                                               <  2:2186982/100‑1 (MQ=38)
        gAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCg                                                                                                 <  2:3187545/97‑1 (MQ=255)
         aaTTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgctt                                                                                             <  1:664733/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)
             tACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgct                                                                                         <  1:2931536/100‑1 (MQ=255)
               cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgcttc                                                                                       >  2:1803406/1‑100 (MQ=255)
               cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgcttc                                                                                       >  1:3576141/1‑100 (MQ=255)
               cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATACGgcttcgcttc                                                                                       >  1:1940006/1‑100 (MQ=255)
                ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCa                                                                                      <  1:2115837/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)
                    gacgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGt                                                                                  <  1:3353644/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‑ATGCGGCTTCGCTTCATTGTTAcc                                                                              <  1:2901863/99‑1 (MQ=255)
                         ccTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCa                                                                             >  1:1898834/1‑100 (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)
                            gctgGGGCGGATTTTCTCCAGCTTCTGTATAGGAAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATAGTTACCACTc                                                                          >  2:2564974/1‑100 (MQ=255)
                               ggggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTg                                                                       >  2:1942290/1‑100 (MQ=255)
                                gggCGGATTTTCTCCAGTTTTTGTATTGGTAAGTAACCGCGTTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGt                                                                      <  1:2618380/100‑1 (MQ=255)
                                  gCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTa                                                                    >  2:336056/1‑100 (MQ=255)
                                   cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                                   >  2:3812627/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)
                                     gATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTc                                                                 <  2:825071/100‑1 (MQ=255)
                                         ttCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCa                                                             <  2:2155098/100‑1 (MQ=255)
                                          tctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCaa                                                            <  2:3376153/100‑1 (MQ=255)
                                           ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCaag                                                           >  1:1843645/1‑99 (MQ=255)
                                           ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                           <  1:2809632/100‑1 (MQ=255)
                                            tcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAcc                                                          >  1:449801/1‑100 (MQ=255)
                                            tcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAcc                                                          <  2:835839/100‑1 (MQ=255)
                                              catcTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCTCTCCTGTTATTCCTCAAc                                                           <  2:1332129/94‑1 (MQ=255)
                                                 cTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCttt                                                      <  1:3441037/99‑1 (MQ=255)
                                                   tCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCttt‑ttt                                                  <  1:1268367/100‑1 (MQ=255)
                                                   tCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCttt‑ttt                                                  <  2:3741883/100‑1 (MQ=255)
                                                     tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTAAG‑ATGCGGATTCGCTTCATAGTTAACAATCCTGTTATTCCTCAACCCTTT‑TTTTa                                                >  1:3679861/1‑100 (MQ=255)
                                                        aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCTG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAAc                                             <  2:2742269/100‑1 (MQ=255)
                                                        aTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAAc                                             >  1:1577788/1‑100 (MQ=255)
                                                          tGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACAt                                           <  2:3550521/100‑1 (MQ=255)
                                                           ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACAtt                                          >  2:1738988/1‑100 (MQ=255)
                                                             taagtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTaa                                        >  2:704327/1‑100 (MQ=255)
                                                               agtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTaaaa                                      >  2:2442872/1‑100 (MQ=255)
                                                               agtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACAACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTaaaa                                      >  1:1165266/1‑100 (MQ=255)
                                                                gtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAAt                                     <  2:2338883/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)
                                                                      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)
                                                                          aagcGAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTaa                           >  1:2444461/4‑100 (MQ=37)
                                                                           tACGAAGCCGCATTCTGACTGTCAG‑ATGCGTCTTTGCTTCATTGTTTCCATTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAAt                          <  1:1219044/100‑1 (MQ=255)
                                                                           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)
                                                                            aCGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTTCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAAtt                         <  1:52329/100‑1 (MQ=255)
                                                                             cGAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAAttt                        >  1:569808/1‑100 (MQ=37)
                                                                              gAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTta                       >  1:800505/1‑100 (MQ=255)
                                                                              gAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTta                       <  1:2884521/100‑1 (MQ=255)
                                                                              gAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTta                       >  2:2203782/1‑100 (MQ=38)
                                                                              gAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTta                       >  1:1467840/1‑100 (MQ=38)
                                                                              gAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCATTT‑TTTTAAACATTAAAATTCTTACGTAATTta                       >  1:1250690/1‑100 (MQ=37)
                                                                              gAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTACTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTta                       >  1:1089561/1‑100 (MQ=37)
                                                                                 gCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATaa                    >  2:3636169/1‑100 (MQ=255)
                                                                                 gCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATaa                    <  1:3509999/100‑1 (MQ=255)
                                                                                 gCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATaa                    >  1:3333425/1‑100 (MQ=38)
                                                                                  ccGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTTCGTAATTTATAAt                   <  2:2422977/100‑1 (MQ=255)
                                                                                  ccGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATaa                    >  1:3877057/1‑99 (MQ=255)
                                                                                  ccGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAAt                   >  2:1607871/1‑100 (MQ=255)
                                                                                   cGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATc                  >  1:2831057/1‑100 (MQ=255)
                                                                                      aTTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCttt               <  2:3755443/100‑1 (MQ=255)
                                                                                      aTTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCttt               <  1:1977625/100‑1 (MQ=255)
                                                                                      aTTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCttt               <  2:1118515/100‑1 (MQ=255)
                                                                                       ttCTGACTGTCAG‑ATGCGGCTTTGCTTCATTGTTACTTCTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTa              <  1:2505299/100‑1 (MQ=255)
                                                                                         cTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaa            >  2:1555363/1‑100 (MQ=255)
                                                                                         cTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaa            >  2:181877/1‑100 (MQ=255)
                                                                                             cTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaaaaaa        <  2:1866467/100‑1 (MQ=255)
                                                                                             cTGTCAG‑ATGCGACTTCGCTTCATTGTTACCACACCTATTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTAGGTAATTTATAATCTTTaaaaaaa        >  2:2588184/1‑100 (MQ=255)
                                                                                              tGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCATCCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAg       <  1:2153585/100‑1 (MQ=255)
                                                                                              tGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAg       <  2:2913161/100‑1 (MQ=255)
                                                                                              tGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAg       <  1:1325330/100‑1 (MQ=255)
                                                                                               gTCAG‑ATGCGGCTTCGCTTCATTGTTACGAATCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAg       <  2:2313643/99‑1 (MQ=255)
                                                                                                tCAG‑ATGCGGCTTCGATTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATATTTAAAAAAAGCa     <  2:383232/100‑1 (MQ=255)
                                                                                                   g‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCAttt  >  1:3912599/1‑100 (MQ=255)
                                                                                                    |                                                                                                    
TCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTT‑TTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTT  >  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: ‑