Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,707,320 A→T 20.0% intergenic (‑196/+2) lepA ← / ← rseC back‑translocating elongation factor EF4, GTPase/SoxR iron‑sulfur cluster reduction factor component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,707,3200AT20.0% 445.7 / 95.7 193intergenic (‑196/+2)lepA/rseCback‑translocating elongation factor EF4, GTPase/SoxR iron‑sulfur cluster reduction factor component
Reads supporting (aligned to +/- strand):  ref base A (89/65);  new base T (18/21);  total (107/86)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 198 total aligned reads displayed.

ACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAA  >  NC_000913/2707221‑2707419
                                                                                                   |                                                                                                   
acactacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGaa                                                                                                     >  2:607998/1‑100 (MQ=255)
   ctacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAAt                                                                                                   >  1:3468756/1‑99 (MQ=255)
   ctacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATc                                                                                                  >  1:4348873/1‑100 (MQ=255)
   ctacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATc                                                                                                  >  2:154866/1‑100 (MQ=255)
   caacaaccGTGAGGAAACGAAATGTTCATACTGTATGTACTCTCCTATCTGGAAACGTCCTCGCATTTGATGTGCAGAATGCAACAAACACAGTGAAAt                                                                                                   >  1:1131434/3‑99 (MQ=255)
    tacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCa                                                                                                 <  2:2652967/100‑1 (MQ=255)
     acaacaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCAc                                                                                                >  1:513668/1‑100 (MQ=255)
      caacaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACt                                                                                               <  2:1783497/100‑1 (MQ=255)
      caacaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACt                                                                                               <  1:2743873/100‑1 (MQ=255)
          aTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCt                                                                                           <  2:3646556/100‑1 (MQ=255)
           ttGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtt                                                                                          <  1:1469753/100‑2 (MQ=255)
            tGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCttt                                                                                         <  2:2540742/100‑3 (MQ=255)
             gAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttg                                                                                        <  2:4146684/100‑4 (MQ=255)
             gAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTcgc                                                                                        >  2:2354454/1‑100 (MQ=255)
               ggCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttgtt                                                                                      <  2:4275849/100‑6 (MQ=255)
               ggCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttgtt                                                                                      <  2:1231553/100‑6 (MQ=255)
               ggCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGt                                                                                      >  1:832161/1‑100 (MQ=255)
                gccaGGAACTGTTCATCCCACATGTAATGTGCTATGGGGAAACGTCCCCGCATTTGTTATGCAAAATGGAACAAAGCCAGTGATTTCACTggagtt                                                                                         <  1:2062882/93‑5 (MQ=255)
                gCGAGGAAATGTTCATACCGTATGGATTGTTGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttgtt                                                                                      <  2:862007/99‑6 (MQ=38)
                gCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGTAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCttt                                                                                         <  1:1293798/96‑3 (MQ=255)
                gCGAGGAAACTTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCCCGCATTTGTTATGCTAATGGCAACAAAGCCAGTGATTTCACTGGCtttgttg                                                                                     <  1:2530666/100‑7 (MQ=25)
                 cGAGGAAATGTTCATACCGTATGGATTGTTGTATCTGGACACTGCCTCGCATTTGTTATGCAAAATGCAACACAGCCAGTGATTTCACTGGCtttgt                                                                                       <  1:2776146/97‑5 (MQ=255)
                 cGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCt                                                                                    <  2:1101303/100‑1 (MQ=255)
                   aGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTc                                                                                  >  2:1977207/1‑100 (MQ=255)
                         tGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGa                                                                            <  1:3251144/100‑1 (MQ=255)
                         tGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGa                                                                            >  2:443227/1‑100 (MQ=255)
                           ttCATACCGTATGTATTTTTATATCTACAAACGTCCTTGCTTTCGTTATCCAATATGCAACAAAGCCAGTGATTTCACTGTCtttgttgc                                                                                    <  1:3851497/90‑8 (MQ=12)
                           ttCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATg                                                                          <  1:2488496/100‑1 (MQ=255)
                           ttCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATg                                                                          <  2:144860/100‑1 (MQ=255)
                            tCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATg                                                                          <  1:1231536/99‑1 (MQ=255)
                            tCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGt                                                                         >  2:2799497/1‑100 (MQ=255)
                             cATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGtt                                                                        >  2:3141901/1‑100 (MQ=255)
                             cATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGtt                                                                        >  1:2826632/1‑100 (MQ=255)
                              aTACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGttt                                                                       >  2:1279277/1‑100 (MQ=255)
                               tACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTc                                                                      >  1:1537940/1‑100 (MQ=255)
                               naCCGTATGGATTGTGGTATCTGGAACCGTCCTCGCATTTGTTATGCAAAATGCAACCAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTc                                                                      >  2:4285010/2‑100 (MQ=255)
                                aCCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCa                                                                     <  2:1498169/100‑1 (MQ=255)
                                 ccGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCaa                                                                    >  1:4263311/1‑100 (MQ=255)
                                 ccGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCaa                                                                    >  2:2695841/1‑100 (MQ=255)
                                   gTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAAt                                                                  <  1:224166/100‑1 (MQ=255)
                                    tATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATc                                                                 >  1:3402974/1‑100 (MQ=255)
                                      tGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGc                                                               <  2:3197191/100‑1 (MQ=255)
                                      tGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGc                                                               <  2:2303025/100‑1 (MQ=255)
                                       ggATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGca                                                              <  2:693299/100‑1 (MQ=255)
                                       ggATTGTGGTATCTGGAAACGTCCTCGCATGTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGc                                                               <  2:4395837/99‑1 (MQ=255)
                                        gATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGcac                                                             >  2:2418627/1‑100 (MQ=255)
                                          ttGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACaa                                                           >  1:2941232/1‑100 (MQ=255)
                                           tgtgGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAg                                                          >  1:3607839/1‑100 (MQ=255)
                                              ggTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGc                                                        >  1:597139/1‑99 (MQ=255)
                                               gTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCa                                                      >  2:3554669/1‑100 (MQ=255)
                                                    tGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAggcgg                                                 >  1:2452501/1‑100 (MQ=255)
                                                     ggAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAggcggc                                                <  1:295181/100‑1 (MQ=255)
                                                     ggAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGCGAAATCACTGGCTCGCGTCTTACGAAGATGTTTCAAATCGCACAAGGCAAggcggc                                                >  1:3714467/1‑100 (MQ=255)
                                                      gAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCa                                               <  2:1889171/100‑1 (MQ=255)
                                                         aCGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATg                                            <  2:2463321/100‑1 (MQ=255)
                                                           gTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGcc                                          <  2:3823027/100‑1 (MQ=255)
                                                             ccTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCAc                                        <  2:330595/100‑1 (MQ=255)
                                                              cTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCATCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCAc                                        >  1:2237674/1‑99 (MQ=255)
                                                               tCGCATTTGTTATGCAAAATGCAACAAAGCCCGTGAAATCACTGGCTCGCGTGTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGc                                      >  1:1893994/1‑100 (MQ=255)
                                                               tCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGc                                      >  1:1648168/1‑100 (MQ=255)
                                                                cGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCt                                     >  2:3297756/1‑100 (MQ=255)
                                                                   aTTTGTTATGCAAAATGCAACAAAGCCAGAGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAATTCGCACAAGGCCAGGCGGCAATGCCACGCTTaa                                  <  1:2635297/100‑1 (MQ=255)
                                                                     ttGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAga                                <  2:832165/100‑1 (MQ=255)
                                                                      tGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAgat                               >  2:2893299/1‑100 (MQ=255)
                                                                      tGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAgat                               <  1:4179811/100‑1 (MQ=255)
                                                                         tATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAgatgat                            >  2:1429741/1‑100 (MQ=255)
                                                                            gCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCgg                         <  2:3601867/100‑1 (MQ=255)
                                                                             cAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCgg                         >  1:1482479/1‑99 (MQ=255)
                                                                             cAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGt                        >  2:790447/1‑100 (MQ=255)
                                                                               aaaTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTg                      <  1:1023796/100‑1 (MQ=255)
                                                                               aaaTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTg                      >  1:1108512/1‑100 (MQ=255)
                                                                                aaTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGc                     <  2:3784058/100‑1 (MQ=255)
                                                                                 aTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTg                      <  1:1412960/98‑1 (MQ=255)
                                                                                 aTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGcc                    >  1:4397243/1‑100 (MQ=255)
                                                                                  tGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCCCGCGGCAATGCCACGCTTAAGCTGATCGGTTGcc                    >  1:1240292/1‑99 (MQ=255)
                                                                                    caacaaAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCAtt                 >  1:3384821/1‑100 (MQ=255)
                                                                                    caacaaAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCAtt                 >  2:1716318/1‑100 (MQ=255)
                                                                                      acaaAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCt               <  2:2175585/100‑1 (MQ=255)
                                                                                       gcgaGCCAGTGATTTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTg              >  2:3042916/4‑100 (MQ=255)
                                                                                       caaAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTg              <  1:62471/100‑1 (MQ=255)
                                                                                       caaAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTg              <  1:3113981/100‑1 (MQ=255)
                                                                                        cgaGCCAGTGATTTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGc             >  2:3488154/3‑100 (MQ=255)
                                                                                        cgaGCCAGTGATTTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGc             >  1:260195/3‑100 (MQ=255)
                                                                                        cgaGCCAGTGATTTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGc             >  2:1357536/3‑100 (MQ=255)
                                                                                        cgaGCCAGTGATTTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCAAGCTTAAGATGATCGGTTGCCATTCTGc             >  1:2345374/3‑100 (MQ=255)
                                                                                         aaGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAAGCCAGGCGGCAATGCCACGCTTAAGATAATCGGTTGCCATTCTGcc            >  1:2276631/1‑100 (MQ=255)
                                                                                         aaGCCAGTGAAATCACTGGCTCGCGTCTCCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGcc            <  1:380208/100‑1 (MQ=255)
                                                                                          aGCCAGTGATTTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGACGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGccc           >  1:396478/1‑100 (MQ=255)
                                                                                          aGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTCGCCATTCTGccc           >  1:1446679/1‑100 (MQ=255)
                                                                                           gCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGcc            >  2:2497311/1‑98 (MQ=255)
                                                                                           gCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCg          <  2:2770121/100‑1 (MQ=255)
                                                                                            ccAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGccc           >  1:714850/1‑98 (MQ=255)
                                                                                             cAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCggg        >  2:3052435/1‑100 (MQ=255)
                                                                                              aGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGgc       >  1:2018515/1‑100 (MQ=255)
                                                                                              aGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGATTAAGATGATCGGTTGCCATTCTGCCCGGgc       >  2:2065751/1‑100 (MQ=255)
                                                                                               gTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGgca      <  2:1452730/100‑1 (MQ=255)
                                                                                                 gAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGgcagc    >  1:1942970/1‑100 (MQ=255)
                                                                                                 gAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGgcagc    >  1:1803328/1‑100 (MQ=255)
                                                                                                   aaTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCaa  >  2:1694690/1‑100 (MQ=255)
                                                                                                   aaTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCaa  >  2:3422093/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAA  >  NC_000913/2707221‑2707419

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