Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,966,637 A→G 84.2% 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,6370AG84.2% 425.7 / 59.5 177S600P (TCA→CCA) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base A (14/14);  new base G (67/82);  total (81/96)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

Only 100 of 181 total aligned reads displayed.

TTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGC  >  NC_000913/2966540‑2966733
                                                                                                 |                                                                                                
ttgCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCg                                                                                                >  1:2192986/1‑100 (MQ=255)
  gCGATCAACGGGCAGAATGTATTGTGTCAGATCGTTGGTGCCAACAGAGATGAAATCGCCCCGCTTTGCCAGGTACGGCAGCAGAGATACCATTGGCgg                                                                                               >  2:1726454/1‑99 (MQ=255)
  gCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTggcggc                                                                                              <  2:1504585/100‑1 (MQ=255)
  gCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTggcggc                                                                                              <  1:241830/100‑1 (MQ=255)
  gCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTggcggc                                                                                              >  2:806509/1‑100 (MQ=255)
   cGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCa                                                                                             >  2:3449051/1‑100 (MQ=255)
   cGAGCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTggcggc                                                                                              <  2:3726139/99‑1 (MQ=255)
    gATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCAc                                                                                            >  2:3193123/1‑100 (MQ=255)
     aTCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACt                                                                                           <  1:1474763/100‑1 (MQ=255)
          cGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAg                                                                                      <  1:1557819/100‑1 (MQ=255)
           ggCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGAAGCGGCAGCATAAATACCATTGGCGGCACTTCCAGc                                                                                     >  2:2108923/1‑100 (MQ=255)
           ggCCAGAATGTATAGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGc                                                                                     >  2:2582833/1‑100 (MQ=255)
              cAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCatg                                                                                  <  2:2449211/100‑1 (MQ=255)
              cAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCatg                                                                                  <  1:2178769/100‑1 (MQ=255)
               aGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCatga                                                                                 >  1:199517/1‑100 (MQ=255)
                 aaTGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAAGCGGCCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCatg                                                                                  <  2:357438/97‑1 (MQ=255)
                  aTGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGc                                                                              <  2:3730258/100‑1 (MQ=255)
                   tGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGcc                                                                             <  1:2984858/100‑1 (MQ=255)
                   tGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGcc                                                                             <  1:994514/100‑1 (MQ=255)
                    gTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCg                                                                            <  2:1879056/100‑1 (MQ=255)
                     tATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCg                                                                            >  1:1286521/1‑99 (MQ=255)
                      aTTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGAt                                                                          <  2:652936/100‑1 (MQ=255)
                       ttGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATa                                                                         <  2:2808074/100‑1 (MQ=255)
                       ttGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACAATTGGCGGCACTTCCAGCATGATGCCGATa                                                                         >  2:1858565/1‑100 (MQ=255)
                        tGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATAc                                                                        <  2:20518/100‑1 (MQ=255)
                         gAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACg                                                                       >  1:829480/1‑100 (MQ=255)
                          aGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGt                                                                      >  1:4028141/1‑100 (MQ=255)
                           gTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTg                                                                     >  1:1784765/1‑100 (MQ=255)
                           gTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTg                                                                     >  2:1865178/1‑100 (MQ=255)
                           gTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTg                                                                     >  1:4154274/1‑100 (MQ=255)
                            tCAGATCGTTCGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTgg                                                                    <  2:2913689/100‑1 (MQ=255)
                              aGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGtt                                                                  >  2:2155233/1‑100 (MQ=255)
                               gATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGttt                                                                 <  1:1520872/100‑1 (MQ=255)
                                aTCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTg                                                                <  2:841499/100‑1 (MQ=255)
                                aTCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTg                                                                >  2:2726955/1‑100 (MQ=255)
                                aTCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTg                                                                <  1:2505675/100‑1 (MQ=255)
                                 tCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTgg                                                               <  1:1144975/100‑1 (MQ=255)
                                 tCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTgg                                                               <  2:4077327/100‑1 (MQ=255)
                                  cGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTggg                                                              <  2:199514/100‑1 (MQ=255)
                                  cGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTggg                                                              >  1:1647486/1‑100 (MQ=255)
                                   gTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGa                                                             <  2:855714/100‑1 (MQ=255)
                                   gTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGa                                                             <  2:2421277/100‑1 (MQ=255)
                                      ggTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAtt                                                          >  2:88989/1‑100 (MQ=255)
                                      ggTGCCAACAAAGATGAAGTCCACCCGCTTTGCCAGATGCGGCAGCATAAATACCAATGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGaa                                                            <  2:252460/98‑1 (MQ=39)
                                         gCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCg                                                       >  2:2975906/1‑100 (MQ=255)
                                           cAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTa                                                     <  1:3457270/100‑1 (MQ=255)
                                            aaCAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATATATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTaa                                                    <  2:2665813/100‑1 (MQ=255)
                                            aaCAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTaa                                                    >  1:95219/1‑100 (MQ=255)
                                              cAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTa                                                     >  2:3319153/1‑97 (MQ=255)
                                                gagaTGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGa                                                >  1:2084905/1‑100 (MQ=255)
                                                  gaTGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGatc                                              <  2:3803645/100‑1 (MQ=255)
                                                   aTGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGatca                                             <  2:2687461/100‑1 (MQ=255)
                                                   aTGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGatca                                             <  1:692036/100‑1 (MQ=255)
                                                   aTGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGatca                                             >  1:2283863/1‑100 (MQ=255)
                                                    tGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGatcat                                            >  1:3633116/1‑100 (MQ=255)
                                                     gAAATCGACCCGCTTTGCCAGATGCGGGAGCATAAATACCATTGGCGGCACTGCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGatca                                             <  2:1133099/98‑1 (MQ=255)
                                                     gAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGatcatc                                           <  2:63843/100‑1 (MQ=255)
                                                     gAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGatcatc                                           >  1:1805707/1‑100 (MQ=255)
                                                          cGACCCGCTTTGCCAGATGCGGCAGCATCAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATctcctc                                      >  2:1935264/1‑100 (MQ=255)
                                                           gACCCGCTTTGCCAGATGCGGCAGCATAAATCCCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCg                                     >  1:2814160/1‑100 (MQ=255)
                                                           gACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCg                                     <  2:2731278/100‑1 (MQ=255)
                                                           gACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCg                                     <  1:392862/100‑1 (MQ=255)
                                                           gACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCg                                     >  1:3477132/1‑100 (MQ=255)
                                                             cccGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGAc                                   <  1:3647156/100‑1 (MQ=255)
                                                              ccGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACt                                  <  2:3488466/100‑1 (MQ=255)
                                                              ccGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACt                                  >  1:3176841/1‑100 (MQ=255)
                                                               cGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACtt                                 >  2:1518920/1‑100 (MQ=255)
                                                               cGCTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACtt                                 <  2:3165684/100‑1 (MQ=255)
                                                                 cTTTGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCa                               <  2:1074494/100‑1 (MQ=255)
                                                                  tttGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCa                               <  1:2419924/99‑1 (MQ=255)
                                                                  tttGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCAc                              <  2:1676305/100‑1 (MQ=255)
                                                                   ttGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACg                             <  2:4028111/100‑1 (MQ=255)
                                                                   ttGCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACg                             <  2:4076403/100‑1 (MQ=255)
                                                                   ttGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACg                             <  1:256751/100‑1 (MQ=255)
                                                                     gCCAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACTTGGTTGGGGAATTTCGTAACCGATCATCTCCTCGACTTCACg                             >  2:1963933/1‑98 (MQ=255)
                                                                      ccAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTcc                          >  1:1236775/1‑100 (MQ=255)
                                                                       cAGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCg                         >  1:3706789/1‑100 (MQ=255)
                                                                        aGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCgg                        <  1:3375931/100‑1 (MQ=255)
                                                                        aGATGCGGCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCgg                        <  2:4130760/100‑1 (MQ=255)
                                                                        aGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCgg                        >  1:3671414/1‑100 (MQ=255)
                                                                            gCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACg                    <  1:3886638/100‑1 (MQ=255)
                                                                            gCGGCAGCATAAATACCATGGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACg                    <  1:3909877/100‑1 (MQ=255)
                                                                             cGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGt                   <  1:3374653/100‑1 (MQ=255)
                                                                              ggCAGCATAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGtt                  <  1:2855742/100‑1 (MQ=255)
                                                                               gcagcaTAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTc                 >  1:385878/1‑100 (MQ=255)
                                                                               gcagcaTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGGCTTCACGTCCGGCACATTc                 >  2:1251191/1‑100 (MQ=255)
                                                                                  ttccgatctccATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAAt              <  2:1115161/91‑1 (MQ=255)
                                                                                  gcaTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAAt              <  1:1725712/100‑1 (MQ=255)
                                                                                   caTAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATc             <  2:4104020/100‑1 (MQ=255)
                                                                                   caTAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATc             >  1:845219/1‑100 (MQ=255)
                                                                                   caTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATc             <  1:3357544/100‑1 (MQ=255)
                                                                                    aTAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCa            >  2:1407404/1‑100 (MQ=255)
                                                                                    aTAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCa            <  1:1518929/100‑1 (MQ=255)
                                                                                    aTAAATACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCa            >  2:1119734/1‑100 (MQ=255)
                                                                                      aaaTACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAgg          >  2:2365188/1‑100 (MQ=255)
                                                                                        aTACCATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAggc         >  2:1382128/1‑99 (MQ=255)
                                                                                           ccATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAAtagatcggaa     >  1:1115166/1‑91 (MQ=255)
                                                                                           ccATTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGgcggcg     <  2:428845/100‑1 (MQ=255)
                                                                                             aTTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTg   <  2:3919037/100‑1 (MQ=255)
                                                                                             aTTGGCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTg   >  2:1315309/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGC  >  NC_000913/2966540‑2966733

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