Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,657,901 A→T 15.5% 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.5% 277.4 / 42.3 110intergenic (+31/‑168)ynfD/ynfEDUF1161 family periplasmic protein/putative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (46/47);  new base T (5/12);  total (51/59)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-02

Only 100 of 113 total aligned reads displayed.

CACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCT  >  NC_000913/1657804‑1658000
                                                                                                 |                                                                                                   
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                                                                                              <  2:1167930/100‑1 (MQ=255)
     gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGc                                                                                              <  1:2201437/100‑1 (MQ=255)
     gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc                                                                                              <  2:1236912/100‑1 (MQ=255)
     gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc                                                                                              >  1:2794782/1‑100 (MQ=255)
     gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc                                                                                              <  1:119942/100‑1 (MQ=255)
     gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc                                                                                              <  1:757511/100‑1 (MQ=255)
     gCGGTAACGTCGCTGCTCCCGCGCAGTCGAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc                                                                                              <  1:2367937/100‑1 (MQ=255)
      cGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCg                                                                                             <  1:911594/100‑1 (MQ=255)
      cGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCg                                                                                             <  1:1227840/100‑1 (MQ=255)
      cGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCg                                                                                             <  2:1174659/100‑1 (MQ=255)
      cGGTAACGTCTCTGCCCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCg                                                                                             <  2:2171333/100‑1 (MQ=255)
       ggTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGTGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgt                                                                                            <  2:2063370/100‑2 (MQ=255)
         tAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAtt                                                                                          <  2:809844/100‑1 (MQ=255)
         tAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAtt                                                                                          <  2:2136557/100‑1 (MQ=255)
           aCGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgttcag                                                                                        <  2:1372485/100‑6 (MQ=38)
           aCGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgttcag                                                                                        <  1:530218/100‑6 (MQ=38)
           aCGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTg                                                                                        <  2:196802/100‑1 (MQ=255)
            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)
                   gCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACCAATCCGGCGATTTGAACATaca                                                                                <  1:427416/100‑1 (MQ=255)
                   gCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATaca                                                                                >  1:2331912/1‑100 (MQ=255)
                   gCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATaca                                                                                <  2:391928/100‑1 (MQ=255)
                    cTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaa                                                                               <  2:1821771/100‑1 (MQ=255)
                    cTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaa                                                                               >  2:2135458/1‑100 (MQ=255)
                      cccGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaaca                                                                             <  1:1507281/100‑1 (MQ=255)
                       ccGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACAt                                                                            >  2:243811/1‑100 (MQ=255)
                         gcgcAGTATAGCAAGGATGGCGCGCCTGCCGAACCGTAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATa                                                                                  <  1:1209397/92‑1 (MQ=255)
                            cAGTCTAGCCAGGCTGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATcc                                                                       <  2:2473140/100‑1 (MQ=255)
                              gTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCa                                                                     >  2:2704182/1‑100 (MQ=255)
                                cTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCAcc                                                                   <  2:1759660/100‑1 (MQ=255)
                                cTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCAcc                                                                   >  1:350414/1‑100 (MQ=255)
                                cTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCAcc                                                                   <  2:2348320/100‑1 (MQ=255)
                                  aGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCtt                                                                 <  1:191391/100‑1 (MQ=255)
                                    ccAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttat                                                               >  2:272190/1‑100 (MQ=255)
                                    ccAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttat                                                               >  2:2048388/1‑100 (MQ=255)
                                       ggATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTAc                                                            >  1:2332886/1‑100 (MQ=255)
                                       ggATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTAc                                                            <  2:1099390/100‑1 (MQ=255)
                                        gATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACt                                                           <  2:2189688/100‑1 (MQ=255)
                                          tGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCa                                                         <  1:772353/100‑1 (MQ=255)
                                          tGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCa                                                         <  1:1577079/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:523622/100‑1 (MQ=255)
                                                   tGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCt                                                >  2:1065872/1‑100 (MQ=255)
                                                    gCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTa                                               >  2:1452914/1‑100 (MQ=255)
                                                    gCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTa                                               <  2:2228683/100‑1 (MQ=255)
                                                      cGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTAtt                                             <  1:824530/100‑1 (MQ=255)
                                                       gAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTg                                            >  2:1055227/1‑100 (MQ=255)
                                                       gAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTg                                            >  2:2293095/1‑100 (MQ=255)
                                                       gAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTg                                            <  2:672205/100‑1 (MQ=255)
                                                         aCCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGat                                          <  2:1583892/100‑1 (MQ=255)
                                                          ccGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGata                                         >  1:50500/1‑100 (MQ=255)
                                                          ccGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGata                                         <  2:2340718/100‑1 (MQ=255)
                                                            gCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATg                                       >  1:2280346/1‑100 (MQ=255)
                                                            gCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATg                                       <  1:1065868/100‑1 (MQ=255)
                                                            gCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATg                                       >  1:2372570/1‑100 (MQ=255)
                                                               gTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAt                                    <  1:822009/100‑1 (MQ=255)
                                                                tAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACGTACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAtt                                   <  1:2072100/100‑1 (MQ=255)
                                                                   cacaTTATCGACTGAACGCCGGATATGACAAATCTGGCGATTTGAACCTACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGGTtaat                                >  2:2178325/1‑100 (MQ=255)
                                                                    acaTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaata                               <  2:1304639/100‑1 (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)
                                                                            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)
                                                                                  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:2509771/1‑100 (MQ=255)
                                                                                             gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACt      >  2:2781811/1‑100 (MQ=255)
                                                                                             gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACt      >  2:417718/1‑100 (MQ=255)
                                                                                                 aaTCCGGCGATTTGAACATACAACATTATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTccc   <  1:1122904/99‑1 (MQ=255)
                                                                                                 aaTCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCct  <  1:205174/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
CACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCT  >  NC_000913/1657804‑1658000

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: