Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,171,134 G→A R217R (CGG→CGA qseC → sensory histidine kinase QseC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,171,1340GA99.7% 1378.4 / ‑4.5 331R217R (CGG→CGAqseCsensory histidine kinase QseC
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (158/172);  total (158/173)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

Only 100 of 331 total aligned reads displayed.

GAACTCGCGCCGCTGAACAAACTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGGCCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAACTTCGTA  >  NC_000913/3171036‑3171231
                                                                                                  |                                                                                                 
gAACTCGCGCCGCTGAACAAACTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGAcc                                                                                                 >  1:5889560/1‑101 (MQ=255)
   cTCGCGCCGCTGAACAAACTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACt                                                                                              >  1:5889565/1‑101 (MQ=255)
   cTCGCGCCGCTGAACAAACTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACt                                                                                              >  2:5889567/1‑101 (MQ=255)
     cgcgCCGCTGAACAAACTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTgg                                                                                            <  1:5889144/101‑1 (MQ=255)
     cgcgCCGATGAACAAACTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTgg                                                                                            >  2:5889570/1‑101 (MQ=255)
       cgccgcTGAACAAACTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGtt                                                                                          >  1:5889572/1‑101 (MQ=255)
         ccgcTGAACAAACTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGa                                                                                        >  2:5889578/1‑101 (MQ=255)
          cgcTGAACAAACTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAg                                                                                       >  1:5889579/1‑101 (MQ=255)
           gcTGAACAAACTGGCGCTGGCACTACGTATGCGTGCCCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTg                                                                                             <  2:5889021/94‑1 (MQ=255)
             tGAACAAACTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCg                                                                                    <  1:5888969/101‑1 (MQ=255)
              gAACAAACTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGc                                                                                   >  2:5889583/1‑101 (MQ=255)
               aaCAAACTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCt                                                                                  >  1:5889586/1‑101 (MQ=255)
                  aaaCTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTaaa                                                                               >  2:5889591/1‑101 (MQ=255)
                   aaCTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAAt                                                                              <  1:5889123/101‑1 (MQ=255)
                    aCTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATc                                                                             <  1:5889061/101‑1 (MQ=255)
                        gcgcTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACt                                                                         <  1:5888950/101‑1 (MQ=255)
                        gcgcTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACt                                                                         >  1:5889594/1‑101 (MQ=255)
                        gcgcTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACt                                                                         <  1:5889042/101‑1 (MQ=255)
                          gcTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGt                                                                       >  2:5889597/1‑101 (MQ=255)
                           cTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGtt                                                                      >  1:5889601/1‑101 (MQ=255)
                              gCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGc                                                                   >  1:5889605/1‑101 (MQ=255)
                                aCTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGccc                                                                 <  2:5888975/101‑1 (MQ=255)
                                 cTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCg                                                                <  2:5889189/101‑1 (MQ=255)
                                  tACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGc                                                               >  2:5889609/1‑101 (MQ=255)
                                  tACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGc                                                               >  1:5889610/1‑101 (MQ=255)
                                   aCGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTc                                                                     <  1:5888513/94‑1 (MQ=255)
                                    cGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGcac                                                             >  2:5889615/1‑101 (MQ=255)
                                      tATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGcacac                                                           <  1:5889050/101‑1 (MQ=255)
                                      tATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGcacac                                                           >  1:5889617/1‑101 (MQ=255)
                                       aTGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGcacaca                                                          <  1:5889078/101‑1 (MQ=255)
                                        tGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACAt                                                         <  2:5889491/101‑1 (MQ=255)
                                         gCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATg                                                        <  1:5888854/101‑1 (MQ=255)
                                          cGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGGCCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGc                                                       <  2:5889135/101‑1 (MQ=255)
                                          cGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGc                                                       >  1:5889622/1‑101 (MQ=255)
                                          cGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGc                                                       <  2:5889127/101‑1 (MQ=255)
                                          cGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGc                                                       <  2:5888808/101‑1 (MQ=255)
                                          cGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGc                                                       <  1:5888825/101‑1 (MQ=255)
                                           gTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCg                                                      <  1:5889017/101‑1 (MQ=255)
                                             gACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCTGCACACATGCgat                                                    >  2:5889627/1‑101 (MQ=255)
                                             gACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCgat                                                    >  1:5889626/1‑101 (MQ=255)
                                             gACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCgat                                                    <  2:5889116/101‑1 (MQ=255)
                                              aCCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCgatg                                                   >  1:5889630/1‑101 (MQ=255)
                                              aCCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCgatg                                                   <  2:5888694/101‑1 (MQ=255)
                                                ccTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCgatgat                                                 <  1:5889227/101‑1 (MQ=255)
                                                 cTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGatgatg                                                <  1:5888885/101‑1 (MQ=255)
                                                  tGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATgg                                               >  1:5889636/1‑101 (MQ=255)
                                                   gACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGt                                              >  2:5889637/1‑101 (MQ=255)
                                                    aCTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGtt                                             <  2:5889156/101‑1 (MQ=255)
                                                    aCTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGtt                                             <  2:5889095/101‑1 (MQ=255)
                                                    aCTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGtt                                             <  1:5889257/101‑1 (MQ=255)
                                                     cTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCgat                                                    >  1:5889641/1‑93 (MQ=255)
                                                     cTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTc                                            <  2:5889029/101‑1 (MQ=255)
                                                     cTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTc                                            <  2:5889281/101‑1 (MQ=255)
                                                       cGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGt                                          >  2:5889647/1‑101 (MQ=255)
                                                       cGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGt                                          >  1:5889646/1‑101 (MQ=255)
                                                        ggAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTg                                         <  2:5889067/101‑1 (MQ=255)
                                                         gAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGa                                        <  2:5888926/101‑1 (MQ=255)
                                                          aaaaaCCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGaa                                       >  1:5889652/1‑101 (MQ=255)
                                                           aaaaCCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAac                                      <  1:5889092/101‑1 (MQ=255)
                                                           aaaaCCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAac                                      >  2:5889656/1‑101 (MQ=255)
                                                             aaCCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAacga                                    <  1:5889284/101‑1 (MQ=255)
                                                             aaCCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAacga                                    >  2:5889661/1‑101 (MQ=255)
                                                             aaCCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAacga                                    <  1:5889250/101‑1 (MQ=255)
                                                               ccACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAacgacg                                  >  1:5889662/1‑101 (MQ=255)
                                                                cACTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGc                                 >  2:5889666/1‑101 (MQ=255)
                                                                 aCTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCt                                >  1:5889668/1‑101 (MQ=255)
                                                                  cTAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACACGCGATGATGGTTCGTGAACGACGCtt                               <  1:5889379/101‑1 (MQ=255)
                                                                   tAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCttt                              <  1:5889172/101‑1 (MQ=255)
                                                                   tAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCttt                              >  2:5889673/1‑101 (MQ=255)
                                                                   tAAACGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCttt                              <  1:5889109/101‑1 (MQ=255)
                                                                    aaaCGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGc                                 >  2:5889677/1‑97 (MQ=255)
                                                                    aaaCGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTa                             <  1:5888864/101‑1 (MQ=255)
                                                                     aaCGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTAc                            >  1:5889678/1‑101 (MQ=255)
                                                                      aCGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTAcc                           <  2:5889103/101‑1 (MQ=255)
                                                                      aCGCGACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTAcc                           <  2:5888994/101‑1 (MQ=255)
                                                                       cgcgACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTTAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCt                          <  2:5889084/101‑1 (MQ=255)
                                                                       cgcgACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCt                          >  1:5889682/1‑101 (MQ=255)
                                                                          gACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCg                       <  2:5889374/101‑1 (MQ=255)
                                                                          gACTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCg                       >  1:5889687/1‑101 (MQ=255)
                                                                           aCTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGa                      <  2:5889226/101‑1 (MQ=255)
                                                                           aCTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGa                      >  2:5889689/1‑101 (MQ=255)
                                                                            cTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGAc                     <  1:5889200/101‑1 (MQ=255)
                                                                            cTGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGAc                     >  2:5889693/1‑101 (MQ=255)
                                                                             tGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACg                    >  1:5889694/1‑101 (MQ=255)
                                                                             tGGCGTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACg                    >  1:5889698/1‑101 (MQ=255)
                                                                                 gTACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGc                >  2:5889704/1‑101 (MQ=255)
                                                                                  tACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCt               <  1:5889165/101‑1 (MQ=255)
                                                                                  tACCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCt               <  2:5888952/101‑1 (MQ=255)
                                                                                   aCCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTc              >  1:5889705/1‑101 (MQ=255)
                                                                                   aCCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTc              <  1:5889191/101‑1 (MQ=255)
                                                                                   aCCCAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTc              <  1:5889154/101‑1 (MQ=255)
                                                                                    cccAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCa             <  2:5889053/101‑1 (MQ=255)
                                                                                     ccAGCGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCAc            >  1:5889708/1‑101 (MQ=255)
                                                                                         cGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAAc        <  2:5889039/101‑1 (MQ=255)
                                                                                         cGAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAAc        >  2:5889712/1‑101 (MQ=255)
                                                                                          gAAGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAACt       >  2:5889717/1‑101 (MQ=255)
                                                                                           aaGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAACtt      >  1:5889718/1‑101 (MQ=255)
                                                                                            aGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAACTTc     <  2:5889310/101‑1 (MQ=255)
                                                                                            aGTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAACTTc     <  2:5889167/101‑1 (MQ=255)
                                                                                             gTGCGACCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTg                                         <  2:5888903/64‑1 (MQ=255)
                                                                                                  |                                                                                                 
GAACTCGCGCCGCTGAACAAACTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGGCCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAACTTCGTA  >  NC_000913/3171036‑3171231

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