Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,684,019 C→T 16.9% R380W (CGG→TGG)  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,0190CT16.9% 251.6 / 39.0 108R380W (CGG→TGG) rstBsensory histidine kinase of RstAB two‑component system
Reads supporting (aligned to +/- strand):  ref base C (41/48);  new base T (8/10);  total (50/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

Only 100 of 109 total aligned reads displayed.

GGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGG  >  NC_000913/1683921‑1684118
                                                                                                  |                                                                                                   
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)
     agagCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAACCCGCGACCATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGAt                                                                                                <  1:2758206/99‑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:1670898/100‑4 (MQ=255)
         cGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgagg                                                                                           <  2:2039659/100‑4 (MQ=255)
         cGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgagg                                                                                           <  2:1624994/100‑4 (MQ=255)
         cGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTc                                                                                           >  1:1446187/1‑100 (MQ=255)
          gACATTAATTGTTGGGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCa                                                                                          <  1:2052085/100‑1 (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:369632/100‑6 (MQ=255)
           aCATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgaggcg                                                                                         <  1:436897/100‑6 (MQ=255)
           aCATTAATTGTTGAATATGAAGGCCCAGGGATTGCCCTAGTAAATCGCGAACATATCTTTGAACCTTTTGTTTTCCTTGATCTCTGCGGGTATCGCTCa                                                                                          >  1:1992279/1‑99 (MQ=255)
             aTTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAAcc                                                                                       <  2:2641255/100‑1 (MQ=255)
             aTTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAAcc                                                                                       >  1:2575047/1‑100 (MQ=255)
             aTTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAAcc                                                                                       >  2:341566/1‑100 (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                                                                         >  1:2223834/1‑100 (MQ=255)
                           atgatgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggct                                                                         >  2:1421069/1‑100 (MQ=255)
                            tgatgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctg                                                                        <  1:1668135/100‑1 (MQ=255)
                              atgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctgg                                                                       <  1:1765392/99‑1 (MQ=255)
                                ggCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctggggc                                                                    >  1:2333044/1‑100 (MQ=255)
                                ggCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctggggc                                                                    <  1:821848/100‑1 (MQ=255)
                                 gCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctggggct                                                                   >  2:2878391/1‑100 (MQ=255)
                                  cccAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctggggctg                                                                  >  2:2553964/1‑100 (MQ=255)
                                  cccAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctggggctg                                                                  >  1:40013/1‑100 (MQ=255)
                                    ttccgatctccccAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGc                                                                <  2:2505169/91‑1 (MQ=255)
                                    cAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGggctggggctgg                                                                 <  2:1653115/99‑1 (MQ=255)
                                    cAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGc                                                                >  1:811074/1‑100 (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)
                                        gATTGCCCCAGAAAACCGCGAACATATCTTTAAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAAt                                                             <  1:1478087/99‑1 (MQ=255)
                                          ttGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGt                                                          <  1:628578/100‑1 (MQ=255)
                                           tGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTc                                                         <  2:2493713/100‑1 (MQ=255)
                                            gCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTcc                                                        >  2:1579/1‑100 (MQ=255)
                                            gCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTcc                                                        <  2:290411/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)
                                             ccccAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCa                                                       >  2:1205198/1‑100 (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:1724615/1‑100 (MQ=255)
                                                        cgcgAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCAc                                            <  2:1587737/100‑1 (MQ=255)
                                                            aaCATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGAAATTGTCCACTCTATAGCACTGGc                                        <  2:2439922/100‑1 (MQ=255)
                                                                  tCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATggg                                  >  1:1276664/1‑100 (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)
                                                                           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)
                                                                                   cGCCTCGATCCCAGCCGGGATCGCTCAACCGGCTGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATtgtg                 <  2:2823663/100‑1 (MQ=255)
                                                                                   cGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATtgtg                 <  2:33819/100‑1 (MQ=255)
                                                                                    gCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGa                <  2:1675045/100‑1 (MQ=255)
                                                                                    gCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGa                >  2:2008648/1‑100 (MQ=255)
                                                                                     nagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGac               >  2:1316293/4‑100 (MQ=255)
                                                                                     gagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGac               >  2:2142589/4‑100 (MQ=255)
                                                                                     gagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGac               >  2:108549/4‑100 (MQ=255)
                                                                                     gagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGac               >  1:1836977/4‑100 (MQ=255)
                                                                                     gagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGac               >  2:1750559/4‑100 (MQ=255)
                                                                                     gagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGac               >  2:17304/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:1308304/1‑100 (MQ=255)
                                                                                      cTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGaca              >  1:2708580/1‑100 (MQ=255)
                                                                                        cGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACAcc            <  2:400863/100‑1 (MQ=255)
                                                                                        cGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACAcc            <  2:1187856/100‑1 (MQ=255)
                                                                                        cGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACAcc            <  2:1357145/100‑1 (MQ=255)
                                                                                         gATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGAAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACAcc            <  1:2552972/99‑1 (MQ=255)
                                                                                                  cGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTgg  <  2:2314585/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGG  >  NC_000913/1683921‑1684118

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: