Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,657,901 A→T 15.9% intergenic (+31/‑168) 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,9010AT15.9% 339.3 / 54.4 139intergenic (+31/‑168)ynfD/ynfEDUF1161 family periplasmic protein/putative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (58/58);  new base T (11/11);  total (69/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-01

Only 100 of 144 total aligned reads displayed.

CCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCT  >  NC_000913/1657803‑1658000
                                                                                                  |                                                                                                   
ccaccaGCGGTAACGTCTCTGCTCCCCCGCAGTCTAGCCAGGATCGCCCGCCTGCCCAACCCCAGTAACAACTTATCGACTCACCCCCGGCTCTGACaaa                                                                                                    >  2:3487936/1‑100 (MQ=255)
      gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGc                                                                                              <  2:2058815/100‑1 (MQ=255)
      gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc                                                                                              <  1:2930046/100‑1 (MQ=255)
       cGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCg                                                                                             <  1:377980/100‑1 (MQ=255)
       cGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCg                                                                                             <  1:2507102/100‑1 (MQ=255)
        ggTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGGAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgt                                                                                            <  2:1022843/100‑2 (MQ=255)
         gTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgt                                                                                            <  1:851203/99‑2 (MQ=255)
           aaCGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAttt                                                                                         <  1:1120699/100‑1 (MQ=255)
            aCGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgttcag                                                                                        <  1:273976/100‑6 (MQ=38)
            aCGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgttcag                                                                                        <  1:1080406/100‑6 (MQ=38)
             cGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGa                                                                                       >  2:3119414/1‑100 (MQ=255)
              gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgttcagtc                                                                                      <  2:1570884/100‑8 (MQ=38)
              gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGaa                                                                                      <  1:136664/100‑1 (MQ=255)
                ctctGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACa                                                                                    >  1:636979/1‑100 (MQ=255)
                    gCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATaca                                                                                <  1:4118676/100‑1 (MQ=255)
                     cacccGCGCAGTCAAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaa                                                                               >  1:2849855/3‑100 (MQ=255)
                     cTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGCACGCCGGATATGACAAATCCGGCGATTTGAACATacaa                                                                               >  2:1980690/1‑100 (MQ=255)
                     cTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaa                                                                               >  2:3262023/1‑100 (MQ=255)
                     cTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaa                                                                               <  2:3915711/100‑1 (MQ=255)
                     cTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTCGAACATac                                                                                 >  2:3274193/1‑98 (MQ=255)
                     cTCCCGAGCAGTCTAGCCAGGATGGCGCGCCTGAAGGGCCGCAGTAAAACATTATCGACTGAACGCCGGCTATGACAAATCCGGCGATTTGACCATacaa                                                                               >  2:3687635/1‑100 (MQ=255)
                       cccGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaac                                                                              <  1:1023912/99‑1 (MQ=255)
                        ccGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACAt                                                                            >  2:3716519/1‑100 (MQ=255)
                            gcAGTCTAGCCAGGATGGCGCGCCTGCCAAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATc                                                                        <  2:2324876/100‑1 (MQ=255)
                             cAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATcc                                                                       <  2:3982038/100‑1 (MQ=255)
                              aGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATccc                                                                      <  2:3698572/100‑1 (MQ=255)
                               gTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTACCACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATccc                                                                      <  2:4197622/99‑1 (MQ=255)
                               gTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCTGGCGATTTGAACATACAACATAATCCCa                                                                     >  1:1829943/1‑100 (MQ=255)
                                tCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCAc                                                                    <  1:2533558/100‑1 (MQ=255)
                                tCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCAc                                                                    <  2:195754/100‑1 (MQ=255)
                                 cTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCAcc                                                                   >  2:2374725/1‑100 (MQ=255)
                                 cTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCAcc                                                                   <  1:1050541/100‑1 (MQ=255)
                                 cTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCAcc                                                                   >  1:1038695/1‑100 (MQ=255)
                                 cTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCAcc                                                                   >  1:2166882/1‑100 (MQ=255)
                                  tAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCt                                                                  <  2:4080022/100‑1 (MQ=255)
                                    gCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCtta                                                                >  1:1789646/1‑100 (MQ=255)
                                      cAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatt                                                              <  1:2141215/100‑1 (MQ=255)
                                      cAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatt                                                              >  2:427416/1‑100 (MQ=255)
                                        ggATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTAc                                                            >  1:1276657/1‑100 (MQ=255)
                                        ggATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTAc                                                            >  1:3928697/1‑100 (MQ=255)
                                         gATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACt                                                           <  1:4055300/100‑1 (MQ=255)
                                            ggCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCAt                                                        <  1:3998771/100‑1 (MQ=255)
                                            ggCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCAt                                                        <  1:2043502/100‑1 (MQ=255)
                                             gcgcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATa                                                       >  1:2345999/1‑100 (MQ=255)
                                             gcgcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATa                                                       >  2:1166678/1‑100 (MQ=255)
                                              cgcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATAc                                                      >  1:2392501/1‑100 (MQ=255)
                                              cgcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATAc                                                      >  2:1053053/1‑100 (MQ=255)
                                               gcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATAcc                                                     <  2:2356739/100‑1 (MQ=255)
                                                 gcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCt                                                   <  1:324867/100‑1 (MQ=255)
                                                 gcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCt                                                   >  2:3340044/1‑100 (MQ=255)
                                                  ccTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCtt                                                  >  2:1117281/1‑100 (MQ=255)
                                                  ccTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCt                                                   >  2:3353311/1‑99 (MQ=255)
                                                   cTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTc                                                 <  1:1244791/100‑1 (MQ=255)
                                                   cTGCCGAACCCCAGTAACACATTATCTACTGAACGCCCGATATCACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTc                                                 <  2:2814279/100‑1 (MQ=255)
                                                         aaCCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGa                                           <  1:794860/100‑1 (MQ=255)
                                                          aCCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGat                                          >  1:3237460/1‑100 (MQ=255)
                                                           ccGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGata                                         >  2:3032655/1‑100 (MQ=255)
                                                             gCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATg                                       >  1:584652/1‑100 (MQ=255)
                                                             gCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATg                                       >  1:3990971/1‑100 (MQ=255)
                                                              cAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATgg                                      <  2:1783701/100‑1 (MQ=255)
                                                              cAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATgg                                      >  1:1380212/1‑100 (MQ=255)
                                                                 tAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAtt                                   <  2:1124435/100‑1 (MQ=255)
                                                                  aaCACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAtt                                   <  1:2013511/99‑1 (MQ=255)
                                                                     acaTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaata                               <  2:4163933/100‑1 (MQ=255)
                                                                          aTCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCt                          <  2:1380210/100‑1 (MQ=255)
                                                                           tCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCtt                         <  1:495/100‑1 (MQ=255)
                                                                            cGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTa                        >  1:1561420/1‑100 (MQ=255)
                                                                             ggcTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTaa                       >  1:2835255/3‑100 (MQ=255)
                                                                             gACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTaa                       <  1:2656073/100‑1 (MQ=255)
                                                                             gACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTaa                       >  1:3441418/1‑100 (MQ=255)
                                                                              aCTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAc                      <  2:3525427/100‑1 (MQ=255)
                                                                               caaatcGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAcc                     >  2:2455025/6‑100 (MQ=38)
                                                                               cTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAcc                     >  1:1528300/1‑100 (MQ=255)
                                                                                tGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAcc                     <  1:2873672/99‑1 (MQ=255)
                                                                                aaatcGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAccc                    >  2:4042367/5‑100 (MQ=38)
                                                                                 gAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCa                   >  1:3051547/1‑100 (MQ=255)
                                                                                 aatcGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTTATATGGATTAATAATTCTTAACCCa                   >  1:3342676/4‑100 (MQ=38)
                                                                                 aatcGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCa                   >  1:3021966/4‑100 (MQ=38)
                                                                                  atcGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCaa                  >  1:3644877/3‑100 (MQ=255)
                                                                                     gCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAAt               >  1:2050193/1‑100 (MQ=255)
                                                                                     gCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAAt               <  1:847192/100‑1 (MQ=255)
                                                                                     gCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTAATATGGATTAATAATTCTTAACCCAAAAt               >  2:742268/1‑100 (MQ=255)
                                                                                      ccGGATATGACACATCCGGCGATTTGAACATACAACATAATACCACCTTAATACTCATACCCTTCTACTGATATGGATTAATAATTCTTAACCCaaa                 <  2:3537656/97‑1 (MQ=255)
                                                                                      ccGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATg              >  2:2271150/1‑100 (MQ=255)
                                                                                       cGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATgg             <  1:1598071/100‑1 (MQ=255)
                                                                                        ggATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATggg            <  1:1389548/100‑1 (MQ=255)
                                                                                        ggATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATggg            <  2:458936/100‑1 (MQ=255)
                                                                                        ggATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATggg            >  2:74475/1‑100 (MQ=255)
                                                                                        ggATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATggg            >  2:1692912/1‑100 (MQ=255)
                                                                                           tctGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATCCCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTa          <  1:2639547/97‑1 (MQ=255)
                                                                                           tatGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGTATCAATAATTCTTAACCCAAAATGGGTAg         >  2:3819249/1‑100 (MQ=255)
                                                                                           tatGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAg         <  2:295006/100‑1 (MQ=255)
                                                                                             tGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGAc       >  1:978005/1‑100 (MQ=255)
                                                                                              gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTaagat      >  2:2639540/1‑96 (MQ=255)
                                                                                              gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACt      >  1:434673/1‑100 (MQ=255)
                                                                                              gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACt      <  2:1051766/100‑1 (MQ=255)
                                                                                                cAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTcc    >  1:1625807/1‑100 (MQ=255)
                                                                                                 aaaTCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTccc   <  1:302049/100‑1 (MQ=255)
                                                                                                  aaTCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTccc   <  2:552754/99‑1 (MQ=255)
                                                                                                  aaTCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCct  >  2:2903331/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCT  >  NC_000913/1657803‑1658000

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: