Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,966,637 A→G 39.0% S600P (TCA→CCA)  ptsP ← PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,966,6370AG39.0% 108.9 / 187.6 172S600P (TCA→CCA) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base A (56/49);  new base G (40/27);  total (96/76)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.30e-01

Only 100 of 173 total aligned reads displayed.

GTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTT  >  NC_000913/2966539‑2966735
                                                                                                  |                                                                                                  
 tttcGATCAACGGCCAGAATGTATTGATTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACg                                                                                                  <  1:3863072/97‑1 (MQ=255)
 ttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGATATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCg                                                                                                  <  2:2972721/100‑1 (MQ=255)
 ttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACg                                                                                                  >  2:280954/1‑100 (MQ=255)
 ttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACg                                                                                                  <  2:3510873/100‑1 (MQ=255)
 ttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACg                                                                                                  <  2:3371035/100‑1 (MQ=255)
 ttgCGATAAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACg                                                                                                  <  1:130986/100‑1 (MQ=255)
  tgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGTAGCATAAATACCATTGACgg                                                                                                 <  2:979672/100‑1 (MQ=255)
  tgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCg                                                                                                  >  1:3708724/1‑99 (MQ=255)
   gCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGCCCCGCTTTGCCAGATGCGGCAGCATAAATACCATTggcggc                                                                                                <  2:1419120/100‑1 (MQ=255)
   gCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGc                                                                                                <  2:3061870/100‑1 (MQ=255)
     gATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCAc                                                                                              >  2:2256861/1‑100 (MQ=255)
     gATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCAc                                                                                              >  1:1354092/1‑100 (MQ=255)
     gATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCa                                                                                               >  1:137343/1‑99 (MQ=255)
     gATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCAc                                                                                              <  1:4242745/100‑1 (MQ=255)
       tCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACt                                                                                             >  1:1538282/1‑99 (MQ=255)
        cAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTc                                                                                           >  1:82399/1‑100 (MQ=255)
         aaCGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCCGATGCGGCAGCATAAATACCATTGGCGGCACTTcc                                                                                          >  2:787539/1‑100 (MQ=255)
         aaCGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTcc                                                                                          <  1:1864102/100‑1 (MQ=255)
           cGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAg                                                                                        >  1:4096515/1‑100 (MQ=255)
           cGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAg                                                                                        <  2:3783051/100‑1 (MQ=255)
            gaccAGAATCTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGc                                                                                       <  1:3397909/98‑1 (MQ=255)
             gCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCa                                                                                      >  2:866789/1‑100 (MQ=255)
              ccAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCat                                                                                     <  1:1133599/100‑1 (MQ=255)
              ccAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCat                                                                                     >  1:2258027/1‑100 (MQ=255)
                 gAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCatgat                                                                                  >  2:18716/1‑100 (MQ=255)
                  aaTGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCatgat                                                                                  <  2:567985/99‑1 (MQ=255)
                    tGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGcc                                                                               <  2:4005973/100‑1 (MQ=255)
                    tGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGcc                                                                               >  2:869147/1‑100 (MQ=255)
                    tGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGcc                                                                               <  2:2260111/100‑1 (MQ=255)
                     gTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCg                                                                              <  2:3682206/100‑1 (MQ=255)
                     gTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCg                                                                              >  2:710145/1‑100 (MQ=255)
                      tATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGa                                                                             >  2:1735478/1‑100 (MQ=255)
                       gttGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGAt                                                                            >  1:2989761/2‑100 (MQ=255)
                          gAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACg                                                                         >  1:3201180/1‑100 (MQ=255)
                           aGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGt                                                                        >  1:2395291/1‑100 (MQ=255)
                            gTCAGATCGTTGGTGCCGACAGAGATGAAATCGAGCGGCTGTGCCAGATGCGCCAGCAGAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTg                                                                       <  1:3953825/100‑1 (MQ=255)
                                gATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGttt                                                                   >  1:3084197/1‑100 (MQ=255)
                                    gTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGa                                                               <  2:1916565/100‑1 (MQ=255)
                                    gTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGa                                                               >  1:4024710/1‑100 (MQ=255)
                                     ttGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGaa                                                              <  2:1316907/100‑1 (MQ=255)
                                         tGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTc                                                          >  2:2496491/1‑100 (MQ=255)
                                          gCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCg                                                         <  1:2214063/100‑1 (MQ=255)
                                             aaCAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTa                                                       <  1:2015415/99‑1 (MQ=255)
                                             aaCAGAGATGAAATCGACCCGCTTTGCCAGATGCGACAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTaa                                                      <  2:1131287/100‑1 (MQ=255)
                                              aCAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATGTCGTAAc                                                     <  2:2250168/100‑1 (MQ=255)
                                               cAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCATCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAAcc                                                    <  2:2996528/100‑1 (MQ=255)
                                               cAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAAcc                                                    >  2:3420930/1‑100 (MQ=255)
                                               cAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTCTGGGTATTTCGTAAcc                                                    >  2:3115921/1‑100 (MQ=255)
                                                agagATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCg                                                   <  1:1941933/100‑1 (MQ=255)
                                                  agaTGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGat                                                 >  2:4096338/1‑100 (MQ=255)
                                                    aTGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGatca                                               <  1:1941952/100‑1 (MQ=255)
                                                    aTGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGatca                                               >  2:2143801/1‑100 (MQ=255)
                                                    aTGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGatca                                               <  1:1501798/100‑1 (MQ=255)
                                                      gAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAACATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGatcatc                                             >  1:2737182/1‑100 (MQ=255)
                                                      gAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGatcatc                                             >  1:1477955/1‑100 (MQ=255)
                                                        aaTCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCAtctc                                           >  1:44450/1‑100 (MQ=255)
                                                         aTCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATctcc                                          <  2:1354096/100‑1 (MQ=255)
                                                         aTCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATctcc                                          >  1:3581704/1‑100 (MQ=255)
                                                         aTCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATctcc                                          >  1:1730875/1‑100 (MQ=255)
                                                          tCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATctcct                                         <  1:1467506/100‑1 (MQ=255)
                                                           cGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATctcctc                                        <  1:2087722/100‑1 (MQ=255)
                                                            gACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCg                                       <  2:1826738/100‑1 (MQ=255)
                                                            gACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCg                                       <  2:935836/100‑1 (MQ=255)
                                                            gACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCg                                       <  2:293251/100‑1 (MQ=255)
                                                               ccGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACt                                    >  2:652261/1‑100 (MQ=255)
                                                               ccGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACt                                    <  2:1480854/100‑1 (MQ=255)
                                                                cGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACtt                                   >  2:4271964/1‑100 (MQ=255)
                                                                     tGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGt                              >  2:2716950/1‑100 (MQ=255)
                                                                     tGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTGTGGGAATTTCGTAACCGATCATCTCCTCGGCTTCACg                               >  2:1943804/1‑99 (MQ=255)
                                                                      gCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTc                             <  2:2162815/100‑1 (MQ=255)
                                                                       ccAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTcc                            >  1:3908805/1‑100 (MQ=255)
                                                                       ccAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTcc                            >  1:884383/1‑100 (MQ=255)
                                                                        cAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCg                           >  2:2571760/1‑100 (MQ=255)
                                                                        cAGATGCGGCAGCATAAATACCATTGACGGCACTTAAAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCg                           >  1:2107895/1‑100 (MQ=255)
                                                                         cttccgatctgcaTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCgg                          <  1:502956/90‑1 (MQ=255)
                                                                          gATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGc                         >  1:2448126/1‑100 (MQ=255)
                                                                          gATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGc                         <  2:1656665/100‑1 (MQ=255)
                                                                          gATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGc                         >  1:533390/1‑100 (MQ=255)
                                                                           aTGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGGATTTCGTAACCGATCATCTCCTCTACTTCAcggccggc                         <  1:3993576/99‑1 (MQ=255)
                                                                            tGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCAc                       >  1:1735663/1‑100 (MQ=255)
                                                                                gcagcaTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTc                   >  1:779520/1‑100 (MQ=255)
                                                                                  agcaTAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGTAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCa                  <  1:972955/99‑1 (MQ=255)
                                                                                  agcaTAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCa                  <  2:2018601/99‑1 (MQ=255)
                                                                                  agcaTAAATACCATTGACGGCACTTCCAGCATGGTGCCGATACGTGGTTTCGGAACTTCGTAACCGATCATCTCCGCGACTTCACGTCCCGCa                        <  1:2029363/93‑1 (MQ=255)
                                                                                  agcaTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCaa                 <  1:28320/100‑1 (MQ=255)
                                                                                   gcaTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCggagatcggaag                >  2:502958/1‑90 (MQ=255)
                                                                                   gcaTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAAt                >  2:1273622/1‑100 (MQ=255)
                                                                                    caTAAATACCATTGACGGCACTTCCAGCATGATGTCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATc               >  2:1207922/1‑100 (MQ=255)
                                                                                     aTAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCa              <  2:723376/100‑1 (MQ=255)
                                                                                     aTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCa              <  2:3449437/100‑1 (MQ=255)
                                                                                     aTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCa              >  1:4095605/1‑100 (MQ=255)
                                                                                       aaaTACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATc               >  1:165167/1‑97 (MQ=255)
                                                                                        aaTACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAggc           >  1:940837/1‑100 (MQ=255)
                                                                                        aaTACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAggc           >  2:4249888/1‑100 (MQ=255)
                                                                                        aaTACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAggc           >  1:1332754/1‑100 (MQ=255)
                                                                                         aTACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAggcg          <  2:2455428/100‑1 (MQ=255)
                                                                                          tACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAggcgg         >  2:3288050/1‑100 (MQ=255)
                                                                                            ccATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGgcggcg       >  1:83135/1‑100 (MQ=255)
                                                                                                tGACGGCACTTCCAGCCTGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCTTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCt   <  2:3301765/100‑1 (MQ=255)
                                                                                                 gACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCtt  >  1:1006831/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
GTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTT  >  NC_000913/2966539‑2966735

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: ‑