Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,577,657 A→G 16.6% intergenic (‑65/+74) gntK ← / ← gntR gluconate kinase 2/d‑gluconate inducible gluconate regulon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,577,6570AG16.6% 463.9 / 77.6 204intergenic (‑65/+74)gntK/gntRgluconate kinase 2/d‑gluconate inducible gluconate regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (82/88);  new base G (15/19);  total (97/107)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 209 total aligned reads displayed.

ACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGG  >  NC_000913/3577558‑3577756
                                                                                                   |                                                                                                   
 cGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAg                                                                                                    >  2:224198/1‑100 (MQ=255)
    aaaTGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCg                                                                                                 <  2:2963527/100‑1 (MQ=255)
    aaaTGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCg                                                                                                 <  1:1751778/100‑1 (MQ=255)
     aaTGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCgg                                                                                                <  1:637245/100‑1 (MQ=255)
       tgtgGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGAc                                                                                              <  2:362829/100‑1 (MQ=255)
        gtggtgATCATGGTTAGTCGTGCTTAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGAc                                                                                              <  2:3298477/99‑1 (MQ=255)
          ggtgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAat                                                                                           >  2:764937/1‑100 (MQ=255)
           gtgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACAatg                                                                                          <  1:676394/100‑2 (MQ=255)
           gtgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACAatg                                                                                          <  2:3736713/100‑2 (MQ=255)
            tgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACAatgt                                                                                         <  1:1503910/100‑3 (MQ=255)
             gATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATc                                                                                        <  2:2258313/100‑1 (MQ=255)
             gATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATc                                                                                        <  2:394525/100‑1 (MQ=255)
             gATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATc                                                                                        >  1:1949738/1‑100 (MQ=255)
              aTCATGGTTAGTCGTGCCCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACAatgt                                                                                         <  1:2877897/98‑3 (MQ=255)
               tCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACAatgttac                                                                                      <  2:1470880/100‑6 (MQ=255)
               tCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACAatgttac                                                                                      <  1:998864/100‑6 (MQ=255)
               tCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAg                                                                                      >  2:3113298/1‑100 (MQ=255)
                 aTGGTTAGTCGTGCTCAGAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACAatgttaccg                                                                                    <  2:3722942/100‑8 (MQ=255)
                  tGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACAatgttaccga                                                                                   <  2:2122930/100‑9 (MQ=38)
                   ggTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGc                                                                                  <  1:3304891/100‑1 (MQ=255)
                      tAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAat                                                                               <  1:3579518/100‑1 (MQ=255)
                        gTCGTGCTCAACGTGCTGCCCCAACAATAGAAGACATAAAACTGTTACGGGTCACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAatat                                                                             <  1:1045898/100‑1 (MQ=255)
                        gTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAatat                                                                             >  1:2547786/1‑100 (MQ=255)
                        gTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAatat                                                                             >  1:837535/1‑100 (MQ=255)
                        gTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAatat                                                                             >  1:4048223/1‑100 (MQ=255)
                            tGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATatccat                                                                         <  1:3449872/100‑1 (MQ=255)
                             gCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATatccatc                                                                        <  2:1253252/100‑1 (MQ=255)
                             gCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATatccatc                                                                        >  1:1300511/1‑100 (MQ=255)
                              cTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATatccatcc                                                                       <  2:426238/100‑1 (MQ=255)
                               tCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATAtccatcca                                                                      <  1:4032772/100‑1 (MQ=255)
                                cAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATccatccat                                                                     >  1:562112/1‑100 (MQ=255)
                                        tgCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAg                                                             >  1:2143550/1‑100 (MQ=255)
                                        tgCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAg                                                             <  2:2437163/100‑1 (MQ=255)
                                        tgCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAg                                                             >  2:1695836/1‑100 (MQ=255)
                                          cccccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTTCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAgtg                                                           >  2:2642263/1‑100 (MQ=255)
                                          cccccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAgtg                                                           <  1:606836/100‑1 (MQ=255)
                                           ccccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAgtgt                                                          >  2:241808/1‑100 (MQ=255)
                                            cccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGta                                                         <  2:1103886/100‑1 (MQ=255)
                                            cccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGta                                                         <  2:2342158/100‑1 (MQ=255)
                                            cccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGta                                                         >  2:2568654/1‑100 (MQ=255)
                                            cccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGta                                                         <  2:3856616/100‑1 (MQ=255)
                                             ccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGtat                                                        <  1:2794829/100‑1 (MQ=255)
                                             ccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGtat                                                        >  1:2278386/1‑100 (MQ=255)
                                              cacaATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGtata                                                       >  2:165020/1‑100 (MQ=255)
                                              cacaATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGtata                                                       <  1:1345599/100‑1 (MQ=255)
                                                caATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGCCAATATCAGATGCAATATCCATCCATGCGTTAAGTGtatatg                                                     <  1:2827941/100‑3 (MQ=255)
                                                 aaTACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAg                                                     >  2:4035390/1‑99 (MQ=255)
                                                  aTACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAgtg                                                   <  1:4144969/100‑1 (MQ=255)
                                                    aCAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGtgtg                                                 <  1:1069773/100‑1 (MQ=255)
                                                     cAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGa                                                <  1:2661315/100‑1 (MQ=255)
                                                        gAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCt                                             >  1:1306769/1‑100 (MQ=255)
                                                        gAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCt                                             >  2:1714334/1‑100 (MQ=255)
                                                         aaTTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTa                                            >  1:2651950/1‑100 (MQ=255)
                                                          aTTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTAc                                           >  1:3008289/1‑100 (MQ=255)
                                                                 aTGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAAt                                    <  2:3451733/100‑1 (MQ=255)
                                                                 aTGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAAt                                    <  2:366955/100‑1 (MQ=255)
                                                                   gTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAAttt                                  >  1:3949010/1‑100 (MQ=255)
                                                                       cGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtgtg                               <  2:2890322/99‑1 (MQ=255)
                                                                       cGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATAAGAAGCAATATACATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTgg                              >  2:3597582/1‑100 (MQ=255)
                                                                        gggTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTggg                             <  2:204371/100‑1 (MQ=255)
                                                                         ggTAACTGTTATCGGTAACCTTGTCCACCATGATAATGTTAGTAGCCATATTCTTCCTTAGCTTTATTGTTTAGGTGTTACCTACTTCCAATTtgtg                               >  2:2937984/1‑97 (MQ=255)
                                                                         ggTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGc                            >  1:2433383/1‑100 (MQ=255)
                                                                          gTAACTGTTATCGGTAACATTGTCCAGGCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCt                           <  1:728558/100‑1 (MQ=255)
                                                                              cTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTaaa                       >  2:1305108/1‑100 (MQ=255)
                                                                              cTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTaaa                       <  1:2890790/100‑1 (MQ=255)
                                                                              cTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTa                         <  2:3359401/98‑1 (MQ=255)
                                                                                gTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATa                     >  1:971548/1‑100 (MQ=255)
                                                                                gTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATa                     >  1:2232994/1‑100 (MQ=255)
                                                                                gTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATa                     >  1:158783/1‑100 (MQ=255)
                                                                                gTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATa                     <  2:114227/100‑1 (MQ=255)
                                                                                gTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATa                     <  2:2294713/100‑1 (MQ=255)
                                                                                  tATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGa                   <  1:790434/100‑1 (MQ=255)
                                                                                   aTCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGAt                  <  1:2013073/100‑1 (MQ=255)
                                                                                   aTCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGAt                  <  1:3924372/100‑1 (MQ=255)
                                                                                    tCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGTCTTAAATAGATc                 <  1:2749115/100‑1 (MQ=255)
                                                                                    tCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATc                 >  2:1371359/1‑100 (MQ=255)
                                                                                    tCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATc                 <  2:4077090/100‑1 (MQ=255)
                                                                                    atcttgataTTGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATc                 <  2:608149/92‑1 (MQ=255)
                                                                                     ttctgataTTGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATcc                >  1:238070/8‑100 (MQ=255)
                                                                                     cGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATcc                >  2:920845/1‑100 (MQ=255)
                                                                                       gTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGc              >  1:2611114/1‑100 (MQ=255)
                                                                                       gTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGc              >  2:2908425/1‑100 (MQ=255)
                                                                                        tgataTTGTCCGGCTGGACAATGTCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGcc             >  2:3425160/5‑100 (MQ=255)
                                                                                        tgataTTGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGcc             >  2:1965742/5‑100 (MQ=255)
                                                                                        tgataTTGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGcc             >  2:830986/5‑100 (MQ=255)
                                                                                        tgataTTGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGcc             >  1:3035896/5‑100 (MQ=255)
                                                                                        tgataTTGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGcc             >  1:2184500/5‑100 (MQ=255)
                                                                                        tgataTTGTCCGGCTGGACAATATCAGAAGCAATATACATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGcc             >  2:3421692/5‑100 (MQ=255)
                                                                                        tgataTTGTCCGGCTGGACAATATCAGAAGAAAGATCCATCCGTGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCg               >  2:1553540/5‑98 (MQ=255)
                                                                                        tAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGATACTTCAAATTTGTGGGCTTAAATAGATCAGcc             >  2:1411218/1‑100 (MQ=255)
                                                                                         gataTTGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGccc            >  1:3061354/4‑100 (MQ=255)
                                                                                         gataTTGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGccc            >  2:2663692/4‑100 (MQ=255)
                                                                                           cATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCgg          <  2:2801387/100‑1 (MQ=255)
                                                                                           cATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCgg          >  2:3638753/1‑100 (MQ=255)
                                                                                           cATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCg           >  2:1997671/1‑99 (MQ=255)
                                                                                            aTTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCGTAAATAGATCCGCCCgg          >  1:3740330/1‑99 (MQ=255)
                                                                                             tcgTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTg        <  2:2550411/98‑1 (MQ=255)
                                                                                               gTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGAc      <  1:1913164/100‑1 (MQ=255)
                                                                                               gTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGAc      >  1:1993543/1‑100 (MQ=255)
                                                                                                tCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACa     >  2:654218/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGG  >  NC_000913/3577558‑3577756

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