Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,620,224 C→G 15.0% intergenic (+17/+14) marB → / ← eamA periplasmic mar operon regulator/cysteine and O‑acetyl‑L‑serine efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,620,2240CG15.0% 366.7 / 54.8 133intergenic (+17/+14)marB/eamAperiplasmic mar operon regulator/cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base C (58/55);  new base G (10/10);  total (68/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

Only 100 of 134 total aligned reads displayed.

CGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAA  >  NC_000913/1620125‑1620321
                                                                                                   |                                                                                                 
    tGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtct                                                                                               >  1:2622008/1‑100 (MQ=255)
     gATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtctc                                                                                              >  2:3750752/1‑100 (MQ=255)
      aTTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCg                                                                                             >  2:1295835/1‑100 (MQ=255)
       tttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCgg                                                                                            >  2:383704/1‑100 (MQ=255)
        ttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCggg                                                                                           <  1:2546641/100‑1 (MQ=255)
         tAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCgggg                                                                                          >  1:530823/1‑100 (MQ=255)
         tAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCgggg                                                                                          >  2:478056/1‑100 (MQ=255)
         tAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCgggg                                                                                          <  1:683658/100‑1 (MQ=255)
         tAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCgggg                                                                                          >  1:1249667/1‑100 (MQ=255)
          aaaTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGc                                                                                         <  1:282882/100‑1 (MQ=255)
          aaaTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGc                                                                                         >  2:562512/1‑100 (MQ=255)
           aaTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGcc                                                                                        <  1:1183379/100‑2 (MQ=255)
             tCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTTTTCTGGCCCCGAGATGTCGGGGccag                                                                                      <  1:2516149/100‑4 (MQ=255)
             tCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCtta                                                                                      <  1:1414866/100‑1 (MQ=255)
             tCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCtta                                                                                      >  2:3016464/1‑100 (MQ=255)
             tCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCtta                                                                                      >  2:3194279/1‑100 (MQ=255)
             tCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTAGCCCCGAGATGTCGGGGccag                                                                                      <  1:597257/100‑4 (MQ=255)
              cacaTGGGTACTGGCAGTGATAAGTCGGATGCTCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccaga                                                                                     <  2:3197262/100‑5 (MQ=255)
              cacaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGcca                                                                                       <  2:1529585/98‑3 (MQ=255)
              cacaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttat                                                                                     >  1:2665282/1‑100 (MQ=255)
                caTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttatta                                                                                   <  1:939085/100‑1 (MQ=255)
                caTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                    >  2:3593594/1‑99 (MQ=255)
                  tGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAAc                                                                                 >  2:308496/1‑100 (MQ=255)
                   gggTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCTGGGCTTATTaaat                                                                                <  1:2803140/100‑3 (MQ=255)
                   gggTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACt                                                                                >  2:1612124/1‑100 (MQ=255)
                   ggctACTGGCAGTGATAAGTCGGCTGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAAc                                                                                 <  2:1895144/96‑1 (MQ=255)
                     gTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGTCCCGACATCTCGGGGCTTATTAACTTc                                                                              <  2:774883/100‑1 (MQ=255)
                     gTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTc                                                                              <  1:2993282/100‑1 (MQ=255)
                     gTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTc                                                                              <  2:3737511/100‑1 (MQ=255)
                       aCTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTccc                                                                            <  1:4001778/100‑1 (MQ=255)
                       aCTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTccc                                                                            <  1:3999891/100‑1 (MQ=255)
                        cTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCACCACGGTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTccc                                                                            <  2:2268027/99‑1 (MQ=255)
                        cTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCa                                                                           >  1:1316615/1‑100 (MQ=255)
                         tGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCAc                                                                          <  2:2061223/100‑1 (MQ=255)
                          ggCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCAcc                                                                         >  1:1795841/1‑100 (MQ=255)
                             aGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCttt                                                                      >  1:2418073/1‑100 (MQ=255)
                              gTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTa                                                                     >  1:1036980/1‑100 (MQ=255)
                                gATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAcc                                                                   <  2:2745659/100‑1 (MQ=255)
                                gATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAcc                                                                   >  2:3344346/1‑100 (MQ=255)
                                gATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAcc                                                                   >  1:3414391/1‑100 (MQ=255)
                                    aGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCtt                                                               >  1:2439560/1‑100 (MQ=255)
                                     gTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTACCTTCCCACCTTTACCGCttt                                                              <  2:2531226/100‑1 (MQ=255)
                                        ggATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACg                                                           <  1:1643305/100‑1 (MQ=255)
                                             cgcTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCAcc                                                      >  2:2467454/1‑100 (MQ=255)
                                                tCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCa                                                   <  2:1298551/100‑1 (MQ=255)
                                                  ggCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTCTTAACTTCCCACCTTTACCGCTTTACGCCACCGCaa                                                  >  1:2557837/1‑99 (MQ=255)
                                                  ggCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAg                                                 <  1:2009527/100‑1 (MQ=255)
                                                  ggCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAg                                                 <  1:1104827/100‑1 (MQ=255)
                                                   gCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGc                                                <  1:1174905/100‑1 (MQ=255)
                                                    cGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTACCTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                               >  1:1619240/1‑100 (MQ=255)
                                                     gTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCa                                              >  1:1438038/1‑100 (MQ=255)
                                                      tGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCa                                              <  2:2232649/99‑1 (MQ=255)
                                                        cccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTATCTTCCCACCTTTTCCGCTTTACGCCACCGCAAGCCAAAt                                           <  1:3885381/100‑1 (MQ=255)
                                                        cccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTGTACCGCTTTACGCCACCGCAAGCCAAAt                                           <  2:2998156/100‑1 (MQ=255)
                                                         ccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                          <  1:2884543/100‑1 (MQ=255)
                                                           tattatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                        >  2:2473801/1‑100 (MQ=255)
                                                            attatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAt                                       <  2:3008345/100‑1 (MQ=255)
                                                              tataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTg                                     <  1:3874200/100‑1 (MQ=255)
                                                              tataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTg                                     >  2:1015350/1‑100 (MQ=255)
                                                              tataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTg                                     <  1:3084652/100‑1 (MQ=255)
                                                               ataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGa                                    <  1:3344366/100‑1 (MQ=255)
                                                                taATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGat                                   <  1:1270806/100‑1 (MQ=255)
                                                                 aaTCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGata                                  <  1:1100306/100‑1 (MQ=255)
                                                                  aTCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATGAACTTCCCACCTTTACCGCTTTACGCCACAGCAAGCCAAATACATTGatat                                 >  1:3684698/1‑100 (MQ=255)
                                                                     aaCACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACa                              >  1:2013323/1‑100 (MQ=255)
                                                                     aaCACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACa                              >  1:2682872/1‑100 (MQ=255)
                                                                     aaCACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACa                              >  2:1790284/1‑100 (MQ=255)
                                                                        acGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                           <  2:4127280/100‑1 (MQ=255)
                                                                           cTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCgg                        <  1:2126062/100‑1 (MQ=255)
                                                                           cTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCgg                        <  1:3194298/100‑1 (MQ=255)
                                                                           cTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCgg                        <  1:680213/100‑1 (MQ=255)
                                                                            tATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGt                       <  2:2682856/100‑1 (MQ=255)
                                                                             aTGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTc                      <  1:1817815/100‑1 (MQ=255)
                                                                             aTGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTc                      <  1:1213811/100‑1 (MQ=255)
                                                                             aTGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAAATTCCCACCTTTACCGCTTTACGCCACCGTAAGCCAAACACATTGATATATAGCACGGTc                      >  2:475716/1‑100 (MQ=255)
                                                                              tGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCa                     >  1:1132035/1‑100 (MQ=255)
                                                                               gTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCAt                    >  2:111501/1‑100 (MQ=255)
                                                                                   tttGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATg                <  1:2344572/100‑1 (MQ=255)
                                                                                    ttGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGa               >  2:434263/1‑100 (MQ=255)
                                                                                    gttaataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCACATCCATTGATATACACCCCGGTCATAATg                >  2:2627453/9‑99 (MQ=37)
                                                                                    gttaataaaccccGAGATGTCGGGGCTTATAAACCTCCCAACTTTACCGCTCTACGCCACCGCAAGCCAAAAACCTCACAATACAACCCGGTCAAAAAGa               >  2:1196131/10‑100 (MQ=255)
                                                                                      taataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgc             >  1:3512566/7‑100 (MQ=255)
                                                                                      taataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgc             >  1:645522/7‑100 (MQ=255)
                                                                                      gTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgc             >  1:3614156/1‑100 (MQ=255)
                                                                                       ttCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgca            >  2:2868399/1‑100 (MQ=255)
                                                                                        tCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcac           >  1:22229/1‑100 (MQ=255)
                                                                                        tCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcac           >  2:1722367/1‑100 (MQ=255)
                                                                                        tCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcac           <  2:350866/100‑1 (MQ=255)
                                                                                        ataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcac           >  1:1799754/5‑100 (MQ=255)
                                                                                         cTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcacc          >  1:2071043/1‑100 (MQ=255)
                                                                                            gCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgca       >  1:3831494/1‑100 (MQ=255)
                                                                                            gCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTCCCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgc        >  2:851921/1‑99 (MQ=255)
                                                                                             ccccGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgcac      <  2:281419/100‑1 (MQ=255)
                                                                                             ccccGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgcac      >  1:1683554/1‑100 (MQ=255)
                                                                                             ccccGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgcac      >  1:139587/1‑100 (MQ=255)
                                                                                             ccccGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccg         >  2:2620807/1‑97 (MQ=255)
                                                                                              cccGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgcacc     <  1:478060/100‑1 (MQ=255)
                                                                                              cccGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgcac      >  2:3461814/1‑99 (MQ=255)
                                                                                                 gACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTaa  >  1:3422911/1‑100 (MQ=255)
                                                                                                  aCATCTCGGGGCGTATTAACTTCCGACCTTTACAGCTTTACGCCACCTCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACAGCACATaa  >  2:1842040/1‑99 (MQ=255)
                                                                                                   |                                                                                                 
CGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAA  >  NC_000913/1620125‑1620321

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: