Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,657,898 A→T 15.8% intergenic (+28/‑171) ynfD → / → ynfE DUF1161 family periplasmic protein/putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,657,8980AT15.8% 292.4 / 43.3 114intergenic (+28/‑171)ynfD/ynfEDUF1161 family periplasmic protein/putative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (47/49);  new base T (5/13);  total (52/62)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-01

Only 100 of 114 total aligned reads displayed.

GCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACT  >  NC_000913/1657800‑1657996
                                                                                                  |                                                                                                  
gCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGAc                                                                                                   <  2:1607542/100‑1 (MQ=255)
gCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGAc                                                                                                   >  1:47408/1‑100 (MQ=255)
 caccacCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCa                                                                                                  <  1:2890487/100‑1 (MQ=255)
 caccacCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGAc                                                                                                   <  1:2604662/99‑1 (MQ=255)
    caccaGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAAt                                                                                               >  1:145104/1‑100 (MQ=255)
      ccaGCGGTAACGTCCCTGCTCCCGCGCAGCCTAGCCCGGATGCCGCTCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATcc                                                                                             >  1:2052338/1‑100 (MQ=255)
        aGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCgg                                                                                           <  2:1352102/100‑1 (MQ=255)
         gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTATCACATTATCGACTGAACGCCGGATCTGTCATATCCGGc                                                                                          <  1:1799479/100‑1 (MQ=255)
         gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGc                                                                                          <  1:2201437/100‑1 (MQ=255)
         gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGc                                                                                          <  2:1167930/100‑1 (MQ=255)
         gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc                                                                                          <  1:119942/100‑1 (MQ=255)
         gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc                                                                                          >  1:2794782/1‑100 (MQ=255)
         gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc                                                                                          <  1:757511/100‑1 (MQ=255)
         gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc                                                                                          <  2:1236912/100‑1 (MQ=255)
         gCGGTAACGTCGCTGCTCCCGCGCAGTCGAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc                                                                                          <  1:2367937/100‑1 (MQ=255)
          cGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCg                                                                                         <  2:1174659/100‑1 (MQ=255)
          cGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCg                                                                                         <  1:911594/100‑1 (MQ=255)
          cGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCg                                                                                         <  2:1215754/100‑1 (MQ=255)
          cGGTAACGTCTCTGCCCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCg                                                                                         <  2:2171333/100‑1 (MQ=255)
           ggTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGTGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgt                                                                                        <  2:2063370/100‑2 (MQ=255)
           ggTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgt                                                                                        <  2:1048728/100‑2 (MQ=255)
             tAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAtt                                                                                      <  2:2136557/100‑1 (MQ=255)
             tAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAtt                                                                                      <  2:809844/100‑1 (MQ=255)
               aCGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgttcag                                                                                    <  2:1372485/100‑6 (MQ=38)
               aCGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgttcag                                                                                    <  1:530218/100‑6 (MQ=38)
                cGTCTCTTCTCCCGCGCATTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGa                                                                                   <  1:1918418/100‑1 (MQ=255)
                 gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGaa                                                                                  <  2:1299109/100‑1 (MQ=255)
                   ctctGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACa                                                                                <  2:1887800/100‑1 (MQ=255)
                   ctctGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACa                                                                                <  2:795399/100‑1 (MQ=255)
                     ctGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATa                                                                              >  2:1111266/1‑100 (MQ=255)
                      tGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATac                                                                             <  2:1170907/100‑1 (MQ=255)
                       gCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATaca                                                                            >  2:1421683/1‑100 (MQ=255)
                       gCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATaca                                                                            >  1:2331912/1‑100 (MQ=255)
                       gCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATaca                                                                            <  2:391928/100‑1 (MQ=255)
                        cTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaa                                                                           >  2:2135458/1‑100 (MQ=255)
                        cTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaa                                                                           <  2:1821771/100‑1 (MQ=255)
                          cccGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaaca                                                                         <  1:1507281/100‑1 (MQ=255)
                           ccGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACAt                                                                        >  2:243811/1‑100 (MQ=255)
                             gcgcAGTATAGCAAGGATGGCGCGCCTGCCGAACCGTAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATa                                                                              <  1:1209397/92‑1 (MQ=255)
                                  gTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCa                                                                 >  2:2704182/1‑100 (MQ=255)
                                    cTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCAcc                                                               <  2:1759660/100‑1 (MQ=255)
                                    cTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCAcc                                                               <  2:2348320/100‑1 (MQ=255)
                                    cTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCAcc                                                               >  1:350414/1‑100 (MQ=255)
                                      aGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCtt                                                             <  1:191391/100‑1 (MQ=255)
                                        ccAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttat                                                           >  2:2048388/1‑100 (MQ=255)
                                        ccAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttat                                                           >  2:272190/1‑100 (MQ=255)
                                           ggATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTAc                                                        >  1:2332886/1‑100 (MQ=255)
                                              tGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCa                                                     <  1:1577079/100‑1 (MQ=255)
                                              tGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCa                                                     <  1:772353/100‑1 (MQ=255)
                                              tGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATACGGCGATTTGAACATACAACATAATCCCACCTTATTACTCa                                                     >  1:954841/1‑100 (MQ=255)
                                                gcgcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCTGCGATTTGAACATACAACATAATCCCACCTTATTACTCATa                                                   <  1:1587353/100‑1 (MQ=255)
                                                 cgcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATAc                                                  >  1:407448/1‑100 (MQ=255)
                                                 cgcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATAc                                                  >  1:1351328/1‑100 (MQ=255)
                                                  gcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATAcc                                                 >  2:1714099/1‑100 (MQ=255)
                                                   cgacTGCCGAACCGCAGTAACAAATTATCGACTGAACCCCGGCTATGACAAATCCATCGACTTGAAAATAAAACAAAATCCCACCTTATTACTAATAcc                                                 >  1:2019865/4‑99 (MQ=255)
                                                    gcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCCGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCt                                               <  1:1243923/100‑1 (MQ=255)
                                                     ccTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGTGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCtt                                              >  1:2218328/1‑100 (MQ=255)
                                                     ccTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCtt                                              >  2:1938579/1‑100 (MQ=255)
                                                     ccTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCtt                                              <  2:523622/100‑1 (MQ=255)
                                                       tGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCt                                            >  2:1065872/1‑100 (MQ=255)
                                                        gCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTa                                           <  2:2228683/100‑1 (MQ=255)
                                                        gCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTa                                           >  2:1452914/1‑100 (MQ=255)
                                                        gCCGAACCGCAGTAACACATTATCGACTGAACGCCACATATGACAAATCCGGCGATTTGAACCTACAACATAATCCCACCTTATTACTCATACCCTTc                                             <  1:1221080/98‑1 (MQ=255)
                                                           gAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTg                                        >  2:2293095/1‑100 (MQ=255)
                                                           gAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTg                                        >  2:1055227/1‑100 (MQ=255)
                                                           gAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTg                                        <  2:672205/100‑1 (MQ=255)
                                                              ccGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGata                                     <  2:2340718/100‑1 (MQ=255)
                                                              ccGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGata                                     >  1:50500/1‑100 (MQ=255)
                                                                gCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATg                                   >  1:2372570/1‑100 (MQ=255)
                                                                gCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATg                                   >  1:2280346/1‑100 (MQ=255)
                                                                gCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATg                                   <  1:1065868/100‑1 (MQ=255)
                                                                   gTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAt                                <  1:822009/100‑1 (MQ=255)
                                                                    tAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACGTACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAtt                               <  1:2072100/100‑1 (MQ=255)
                                                                       cacaTTATCGACTGAACGCCGGATATGACAAATCTGGCGATTTGAACCTACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGGTtaat                            >  2:2178325/1‑100 (MQ=255)
                                                                         caTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaataa                          >  1:915431/1‑100 (MQ=255)
                                                                          aTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATCCCCTTCTATTGATATGGATtaataa                          <  1:1087886/99‑1 (MQ=255)
                                                                          aTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTaataat                         <  2:2753545/100‑1 (MQ=255)
                                                                             aTCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCt                      <  2:1348933/100‑1 (MQ=255)
                                                                               gttcaaatcGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTa                    >  2:1263553/9‑100 (MQ=37)
                                                                               cGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTa                    >  1:2473487/1‑100 (MQ=255)
                                                                                ttcaaatcGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTaa                   >  2:2466714/8‑100 (MQ=37)
                                                                                gACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTaa                   >  2:754210/1‑100 (MQ=255)
                                                                                 aCTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAc                  <  2:211749/100‑1 (MQ=255)
                                                                                 aCTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAc                  >  2:2892174/1‑100 (MQ=255)
                                                                                  tctaaCGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAcc                 <  1:1385130/97‑1 (MQ=255)
                                                                                    gAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCa               >  1:854997/1‑100 (MQ=255)
                                                                                     aaCGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAccaga              >  2:1385135/1‑97 (MQ=255)
                                                                                      tcGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCaaa             >  1:613727/2‑100 (MQ=255)
                                                                                      aCGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCaaa             >  1:1139917/1‑100 (MQ=255)
                                                                                       cGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCaaaa            >  1:1458110/1‑100 (MQ=255)
                                                                                        gCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAAt           >  1:1076595/1‑100 (MQ=255)
                                                                                         ccGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATg          >  2:2217334/1‑100 (MQ=255)
                                                                                           ggATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATggg        >  2:1677033/1‑100 (MQ=255)
                                                                                             atatGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTa      <  2:145105/100‑1 (MQ=255)
                                                                                              tatGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAg     >  1:2895310/1‑100 (MQ=255)
                                                                                                tGACAAATCCGGCGATTTGAACATACAATATAAGCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGAc   <  1:1186415/100‑1 (MQ=255)
                                                                                                tGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGAc   <  2:2765900/100‑1 (MQ=255)
                                                                                                 gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACt  >  2:2781811/1‑100 (MQ=255)
                                                                                                 gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACt  >  2:2509771/1‑100 (MQ=255)
                                                                                                 gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACt  >  2:417718/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
GCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACT  >  NC_000913/1657800‑1657996

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