Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,684,016 A→G 16.6% S379G (AGC→GGC)  rstB → sensory histidine kinase of RstAB two‑component system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,684,0160AG16.6% 238.3 / 39.5 110S379G (AGC→GGC) rstBsensory histidine kinase of RstAB two‑component system
Reads supporting (aligned to +/- strand):  ref base A (42/49);  new base G (8/10);  total (51/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.76e-01

Only 100 of 110 total aligned reads displayed.

GGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACC  >  NC_000913/1683919‑1684108
                                                                                                 |                                                                                            
ggggAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAg                                                                                             <  2:530042/99‑1 (MQ=255)
ggggAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGc                                                                                            >  2:1999746/1‑100 (MQ=255)
  ggAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCg                                                                                          >  2:146870/1‑100 (MQ=255)
   gAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCgg                                                                                         <  2:1570711/100‑1 (MQ=255)
     aTAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGa                                                                                       >  2:412749/1‑100 (MQ=255)
      tatagCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAACACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGAt                                                                                      <  1:2103721/97‑1 (MQ=255)
      tAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCgggtt                                                                                      >  1:1767585/1‑98 (MQ=255)
      tAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGAt                                                                                      <  2:2882498/100‑1 (MQ=255)
      tAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCCGGAt                                                                                      >  1:923651/1‑100 (MQ=255)
       agagCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATc                                                                                     <  2:2416154/100‑1 (MQ=255)
        gagCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCg                                                                                    <  2:1245377/100‑1 (MQ=255)
          gCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCt                                                                                  >  2:2017810/1‑100 (MQ=255)
           cGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgagg                                                                                 <  2:1153563/100‑4 (MQ=255)
           cGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgagg                                                                                 <  2:1624994/100‑4 (MQ=255)
           cGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgagg                                                                                 <  2:1670898/100‑4 (MQ=255)
           cGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgagg                                                                                 <  2:2039659/100‑4 (MQ=255)
           cGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTc                                                                                 >  1:1446187/1‑100 (MQ=255)
            gACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgaggc                                                                                <  2:1387185/100‑5 (MQ=255)
            gACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgaggc                                                                                <  2:1172784/100‑5 (MQ=255)
            gACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgaggc                                                                                <  1:209467/100‑5 (MQ=255)
            gACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCa                                                                                >  2:637871/1‑100 (MQ=255)
             aCATTAATTGTTGAGGATGATGGTCCAGGGATTGCCCCAGTAAACCGCGACCAGATCTTTGAGCCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCaa                                                                               <  1:1035031/100‑1 (MQ=255)
             aCATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgaggcg                                                                               <  1:436897/100‑6 (MQ=255)
             aCATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgaggcg                                                                               <  1:369632/100‑6 (MQ=255)
             aCATTAATTGTTGAATATGAAGGCCCAGGGATTGCCCTAGTAAATCGCGAACATATCTTTGAACCTTTTGTTTTCCTTGATCTCTGCGGGTATCGCTCa                                                                                >  1:1992279/1‑99 (MQ=255)
               aTTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAAcc                                                                             >  2:341566/1‑100 (MQ=255)
               aTTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAAcc                                                                             >  1:2575047/1‑100 (MQ=255)
               aTTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAAcc                                                                             <  2:2641255/100‑1 (MQ=255)
                 tAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAACCCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgaggcgaaca                                                                           <  2:1198160/100‑10 (MQ=38)
                  aaTTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACcggc                                                                          >  1:364435/1‑100 (MQ=255)
                     tgttgAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCggcggc                                                                       <  2:1795275/100‑1 (MQ=255)
                      gttgAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCt                                                                      >  1:1226509/1‑100 (MQ=255)
                       ttgAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTg                                                                     >  2:1527958/1‑100 (MQ=255)
                         gAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCg                                                                   >  1:1241544/1‑100 (MQ=255)
                           ggATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCggg                                                                 >  2:1619413/1‑100 (MQ=255)
                             atgatgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggct                                                               >  2:1421069/1‑100 (MQ=255)
                             atgatgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggct                                                               >  1:2223834/1‑100 (MQ=255)
                              tgatgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctg                                                              <  1:1668135/100‑1 (MQ=255)
                                atgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctgg                                                             <  1:1765392/99‑1 (MQ=255)
                                  ggCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctggggc                                                          <  1:821848/100‑1 (MQ=255)
                                  ggCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctggggc                                                          >  1:2333044/1‑100 (MQ=255)
                                   gCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctggggct                                                         >  2:2878391/1‑100 (MQ=255)
                                    cccAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctggggctg                                                        >  1:40013/1‑100 (MQ=255)
                                    cccAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctggggctg                                                        >  2:2553964/1‑100 (MQ=255)
                                      ttccgatctccccAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGc                                                      <  2:2505169/91‑1 (MQ=255)
                                      cAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGggctggggctgg                                                       <  2:1653115/99‑1 (MQ=255)
                                      cAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTCCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGc                                                      <  1:1270697/100‑1 (MQ=255)
                                       aGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCa                                                     <  2:763311/100‑1 (MQ=255)
                                         ggATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAAt                                                   <  2:2838362/100‑1 (MQ=255)
                                          gATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAAtt                                                  <  2:2707225/100‑1 (MQ=255)
                                          gATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAAtt                                                  >  2:2811888/1‑100 (MQ=255)
                                            ttGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGt                                                <  1:628578/100‑1 (MQ=255)
                                              gCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTcc                                              >  2:1005269/1‑100 (MQ=255)
                                               ncccAGAAAACCGCGAACACATCTTTGAACCTTTTGTTCGCCCCGATCCCAGCCGGCATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTc                                               >  1:1504017/2‑98 (MQ=255)
                                               ccccAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCagatcggaa                                             >  1:2505158/1‑92 (MQ=255)
                                                  cAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGTCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCAc                                            <  1:1213318/98‑1 (MQ=255)
                                                   aGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCActct                                         >  1:1727582/1‑100 (MQ=255)
                                                      aaaCCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCtgta                                      <  1:1952228/100‑1 (MQ=255)
                                                      aaaCCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCtata                                      >  1:2770273/1‑100 (MQ=255)
                                                       aaCCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAg                                     <  1:2811878/100‑1 (MQ=255)
                                                       aaCCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAg                                     <  1:637870/100‑1 (MQ=255)
                                                         ccGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCa                                   >  2:1070749/1‑100 (MQ=255)
                                                         ccGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCa                                   >  2:1724615/1‑100 (MQ=255)
                                                          cgcgAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCAc                                  <  2:1587737/100‑1 (MQ=255)
                                                              aaCATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGAAATTGTCCACTCTATAGCACTGGc                              <  2:2439922/100‑1 (MQ=255)
                                                                   atCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATgg                         <  2:2880223/100‑1 (MQ=255)
                                                                    tCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATggg                        >  1:1276664/1‑100 (MQ=255)
                                                                     cTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGc                       <  1:339065/100‑1 (MQ=255)
                                                                     cTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGc                       <  1:1188899/100‑1 (MQ=255)
                                                                      tttGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCg                      <  2:1695984/100‑1 (MQ=255)
                                                                        tGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGt                    <  2:626717/100‑1 (MQ=255)
                                                                          aaCCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTAc                  <  2:567138/100‑1 (MQ=255)
                                                                            ccTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGACTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACgg                <  1:1516936/100‑1 (MQ=255)
                                                                             cTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGt               <  1:9825/100‑1 (MQ=255)
                                                                             cTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGt               <  1:2142780/100‑1 (MQ=255)
                                                                             cTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACAGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGt               <  2:1334257/100‑1 (MQ=255)
                                                                                 tGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAAt           <  2:1338257/100‑1 (MQ=255)
                                                                                   ttCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGAACGGTACGGTTCAt           >  1:1068197/1‑98 (MQ=255)
                                                                                     ttgagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATtgtg       >  1:2126828/6‑100 (MQ=255)
                                                                                     cGCCTCGATCCCAGCCGGGATCGCTCAACCGGCTGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATtgtg       <  2:2823663/100‑1 (MQ=255)
                                                                                     cGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATtgtg       <  2:33819/100‑1 (MQ=255)
                                                                                      gCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGa      >  2:2008648/1‑100 (MQ=255)
                                                                                      gCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGa      <  2:1675045/100‑1 (MQ=255)
                                                                                       nagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGac     >  2:1316293/4‑100 (MQ=255)
                                                                                       gagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGac     >  2:2142589/4‑100 (MQ=255)
                                                                                       gagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGac     >  1:1836977/4‑100 (MQ=255)
                                                                                       gagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGac     >  2:108549/4‑100 (MQ=255)
                                                                                       gagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGac     >  2:17304/4‑100 (MQ=255)
                                                                                       gagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGac     >  2:1687930/4‑100 (MQ=255)
                                                                                       ccTCGCTCCCAGCCGAGATAGCAAAAAAAAAGGCTGCGGGCTGGGGCTGGCAATTGTACAATCTATAGCACTGGCAATGGGCGGTACGGAGAATTGTGa      >  1:1867575/1‑99 (MQ=255)
                                                                                       ccTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGTCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGac     >  2:366832/1‑100 (MQ=255)
                                                                                       ccTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGac     <  1:2783774/100‑1 (MQ=255)
                                                                                        cccGATCCCAGCCGGGATCGCTCAACCGGCGCCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGaca    >  2:1965383/3‑100 (MQ=255)
                                                                                        cTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGaca    >  2:2462464/1‑100 (MQ=255)
                                                                                        cTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGaca    >  1:2708580/1‑100 (MQ=255)
                                                                                        cTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGaca    >  1:1308304/1‑100 (MQ=255)
                                                                                          cGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACAcc  <  2:400863/100‑1 (MQ=255)
                                                                                          cGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACAcc  <  2:1357145/100‑1 (MQ=255)
                                                                                          cGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACAcc  <  2:1187856/100‑1 (MQ=255)
                                                                                           gATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGAAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACAcc  <  1:2552972/99‑1 (MQ=255)
                                                                                                 |                                                                                            
GGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACC  >  NC_000913/1683919‑1684108

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: