Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,535,548 G→T 7.3% F107L (TTC→TTA envZ ← sensory histidine kinase EnvZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,535,5480GT 7.3% 307.3 / 12.8 96F107L (TTC→TTAenvZsensory histidine kinase EnvZ
Reads supporting (aligned to +/- strand):  ref base G (47/42);  new base T (2/5);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01

ACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAA  >  NC_000913/3535451‑3535644
                                                                                                 |                                                                                                
aCAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGa                                                                                                 <  1:1565714/99‑1 (MQ=255)
 cAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGaa                                                                                                <  1:1259518/99‑1 (MQ=255)
  aGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAAt                                                                                               <  2:2814976/99‑1 (MQ=255)
  aGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTAACCTCAACGCGCACTTCCGACGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAAt                                                                                               <  2:1604833/99‑1 (MQ=255)
   gCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAAtt                                                                                              >  2:424087/1‑99 (MQ=255)
   gCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAAtt                                                                                              >  1:1492049/1‑99 (MQ=255)
     cAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCa                                                                                            >  1:1430650/1‑99 (MQ=255)
     cAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCa                                                                                            <  1:1396829/99‑1 (MQ=255)
     cAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCa                                                                                            >  2:1386180/1‑99 (MQ=255)
      aGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCAt                                                                                           >  1:253801/1‑99 (MQ=255)
       ggTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATa                                                                                          <  1:3006343/99‑1 (MQ=255)
        gTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAg                                                                                         <  2:497159/99‑1 (MQ=255)
        gTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAg                                                                                         <  2:1384604/99‑1 (MQ=255)
        gTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAg                                                                                         >  2:916939/1‑99 (MQ=255)
          tttCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAgtg                                                                                       <  2:1926857/99‑1 (MQ=255)
          tttCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAgtg                                                                                       <  2:1845910/99‑1 (MQ=255)
             cAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTg                                                                                    <  2:1346792/99‑1 (MQ=255)
              aGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGGATTCATAGTGTTgc                                                                                   <  2:1852146/99‑1 (MQ=255)
                 cAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAATAATTCATAGTGTTGCGcc                                                                                <  1:2252176/99‑1 (MQ=255)
                 cAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGcc                                                                                <  1:1554106/99‑1 (MQ=255)
                   gacgacAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCa                                                                              >  1:761699/1‑99 (MQ=255)
                    acgacAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAATAATTCATAGTGTTGCGCCCaa                                                                             <  1:2808570/99‑1 (MQ=255)
                       acAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGc                                                                          >  1:2847224/1‑99 (MQ=255)
                       acAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGc                                                                          >  2:118291/1‑99 (MQ=255)
                        cAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCa                                                                         <  2:283001/99‑1 (MQ=255)
                        cAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCa                                                                         <  1:254075/99‑1 (MQ=255)
                          ggCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGa                                                                       >  2:343444/1‑99 (MQ=255)
                          ggCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGa                                                                       >  1:825360/1‑99 (MQ=255)
                           gCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGAc                                                                      <  1:337868/99‑1 (MQ=255)
                           gCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGAc                                                                      <  1:1714864/99‑1 (MQ=255)
                           gCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGAc                                                                      <  1:1204180/99‑1 (MQ=255)
                            cGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAATAATTCATAGTGTTGCGCCCAACGCAGAcc                                                                     <  1:2563978/99‑1 (MQ=255)
                              aaCTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTg                                                                   >  1:2144644/1‑99 (MQ=255)
                               aCTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAATAATTCATAGTGTTGCGCCCAACGCAGACCTGc                                                                  >  1:2825541/1‑99 (MQ=255)
                               aCTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGc                                                                  <  1:368593/99‑1 (MQ=255)
                                 ttttGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCct                                                                <  1:736371/99‑1 (MQ=255)
                                  tttGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAATAATTCATAGTGTTGCGCCCAACGCAGACCTGCctc                                                               <  2:513059/99‑1 (MQ=255)
                                     gttgACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAATAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTc                                                            <  1:869810/99‑1 (MQ=255)
                                     gttgACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTc                                                            >  2:2055442/1‑99 (MQ=255)
                                      ttgACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCg                                                           >  1:769609/1‑99 (MQ=255)
                                       tgACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCgg                                                          <  1:260967/99‑1 (MQ=255)
                                        gACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGc                                                         >  2:1961819/1‑99 (MQ=255)
                                           cTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGc                                                      >  1:2882425/1‑99 (MQ=255)
                                             cAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCt                                                    <  1:1628577/99‑1 (MQ=255)
                                             cAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCt                                                    <  2:465962/99‑1 (MQ=255)
                                               aCGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCACAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCg                                                  >  2:833800/1‑99 (MQ=255)
                                                  cgcACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTg                                               <  2:2190419/99‑1 (MQ=255)
                                                      cTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGt                                           <  2:1248660/99‑1 (MQ=255)
                                                      cTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGt                                           >  2:1086172/1‑99 (MQ=255)
                                                      cTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGt                                           >  1:1609844/1‑99 (MQ=255)
                                                      cTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGt                                           <  1:1029680/99‑1 (MQ=255)
                                                        tCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTag                                         <  1:392621/99‑1 (MQ=255)
                                                        tCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTag                                         >  1:2157691/1‑99 (MQ=255)
                                                        tCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTag                                         >  1:1859144/1‑99 (MQ=255)
                                                          cGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTagag                                       >  2:1154764/1‑99 (MQ=255)
                                                           gTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTagaga                                      >  1:1276693/1‑99 (MQ=255)
                                                            tCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTagagag                                     >  2:2964652/1‑99 (MQ=255)
                                                              gggCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTagagagag                                   >  2:975915/1‑99 (MQ=255)
                                                              gggCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTagagagag                                   <  1:1234868/99‑1 (MQ=255)
                                                              gggCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTagagagag                                   >  2:1057785/1‑99 (MQ=255)
                                                               ggCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAgagagaga                                  <  1:1266959/99‑1 (MQ=255)
                                                               ggCCGCCCAGTTGCTGCGCCATCCGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAgagagaga                                  >  1:1080183/1‑99 (MQ=255)
                                                                 ccgccCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATc                                <  2:2980100/99‑1 (MQ=255)
                                                                  cgccCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATcc                               <  1:1026309/99‑1 (MQ=255)
                                                                  cgccCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATcc                               >  2:1946931/1‑99 (MQ=255)
                                                                   gccCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATccc                              <  2:2809981/99‑1 (MQ=255)
                                                                     ccAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCa                            <  1:2317962/99‑1 (MQ=255)
                                                                      cAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAg                           >  1:685334/1‑99 (MQ=255)
                                                                       aGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGc                          <  2:522729/99‑1 (MQ=255)
                                                                       aGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGc                          <  1:455287/99‑1 (MQ=255)
                                                                         ttGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTc                        >  1:2117021/1‑99 (MQ=255)
                                                                         ttGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTc                        >  1:2443617/1‑99 (MQ=255)
                                                                         ttGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTc                        <  1:1394486/99‑1 (MQ=255)
                                                                           gctgcGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCAc                      >  2:1705572/1‑99 (MQ=255)
                                                                           gctgcGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCAc                      <  2:1352070/99‑1 (MQ=255)
                                                                            ctgcGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACg                     >  1:2083347/1‑99 (MQ=255)
                                                                            ctgcGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACg                     >  1:2916064/1‑99 (MQ=255)
                                                                              gcgcCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGt                   <  1:1940482/99‑1 (MQ=255)
                                                                                 ccATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGa                >  2:2734679/1‑99 (MQ=255)
                                                                                 ccATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGa                >  1:2067563/1‑99 (MQ=255)
                                                                                  cATCTGATGGCTTAATAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGAt               >  1:716228/1‑99 (MQ=255)
                                                                                  cATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGAt               <  2:1034840/99‑1 (MQ=255)
                                                                                  cATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGAt               >  2:2122411/1‑99 (MQ=255)
                                                                                  cATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGAt               <  1:407338/99‑1 (MQ=255)
                                                                                    tCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGAtct             >  1:1500526/1‑99 (MQ=255)
                                                                                    tCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGAtct             >  2:2798598/1‑99 (MQ=255)
                                                                                     cTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGAtctc            >  1:1748310/1‑99 (MQ=255)
                                                                                     cTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGAtctc            <  2:168736/99‑1 (MQ=255)
                                                                                       gATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTccc          >  1:553089/1‑99 (MQ=255)
                                                                                       gATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTccc          >  2:852873/1‑99 (MQ=255)
                                                                                       gATGGATTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTccc          >  2:2350401/1‑99 (MQ=255)
                                                                                        aTGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCg         >  1:1653416/1‑99 (MQ=255)
                                                                                        aTGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCg         <  1:3028883/99‑1 (MQ=255)
                                                                                          ggCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGAc       >  2:160387/1‑99 (MQ=255)
                                                                                            cTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACgg     >  1:86727/1‑99 (MQ=255)
                                                                                               aaGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGaaa  >  1:1086515/1‑99 (MQ=255)
                                                                                                 |                                                                                                
ACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAA  >  NC_000913/3535451‑3535644

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