Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,000,180 A→G 6.8% Q274R (CAG→CGG) ‡ PP_3528 → ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,000,1800AG 6.8% 344.5 / 10.9 117Q274R (CAG→CGG) ‡PP_3528ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (55/54);  new base G (4/4);  total (59/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 = 2.77e-01

Only 100 of 119 total aligned reads displayed.

CTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGAAGCGGTCTACCCGG  >  NC_002947/4000031‑4000329
                                                                                                                                                     |                                                                                                                                                     
 ttCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAg                                                                                                                                                      >  1:1196368/1‑150 (MQ=255)
   cTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGt                                                                                                                                                    <  2:1952643/150‑1 (MQ=255)
     agcTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCATGCCGAAGTGGCAGCGGTGGTGGCAACGCAGCTGGGGATAAGCAGACAGGTgg                                                                                                                                                  >  2:1169052/3‑149 (MQ=255)
     aCCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCCGTGGGCAAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGc                                                                                                                                                  >  2:685326/1‑150 (MQ=255)
       cTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTg                                                                                                                                                <  2:1506753/150‑1 (MQ=255)
         gTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGa                                                                                                                                              >  2:654248/1‑150 (MQ=255)
          tCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCt                                                                                                                                                 >  1:771952/1‑146 (MQ=255)
          tCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCt                                                                                                                                                 <  2:771952/146‑1 (MQ=255)
            cacaGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGt                                                                                                                                           <  2:1171830/150‑1 (MQ=255)
             acaGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAg                                                                                                                                                      <  2:2315491/138‑1 (MQ=255)
             acaGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAg                                                                                                                                                      >  1:2315491/1‑138 (MQ=255)
              cagccagcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGGGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTgg                                                                                                                                         <  1:654248/150‑1 (MQ=255)
                gccagcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGc                                                                                                                                                  <  1:1603071/139‑1 (MQ=255)
                gccagcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGc                                                                                                                                                  >  2:1603071/1‑139 (MQ=255)
                 ccagcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCGCCCAGGTGCTGGAAGTGGc                                                                                                                                        >  2:1237942/1‑148 (MQ=255)
                 ccagcTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGc                                                                                                                                        <  1:1237942/148‑1 (MQ=255)
                         gCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACgg                                                                                                                              <  2:2288917/150‑1 (MQ=255)
                          ccGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGgc                                                                                                                             >  1:480913/1‑150 (MQ=255)
                            gACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTgg                                                                                                                                               >  2:2186595/1‑130 (MQ=255)
                            gACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTgg                                                                                                                                               <  1:2186595/130‑1 (MQ=255)
                               cAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAAc                                                                                                                        >  2:1763699/1‑150 (MQ=255)
                               cAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAAc                                                                                                                        <  2:862245/150‑1 (MQ=255)
                                         ctctGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCg                                                                                                              >  2:2130200/1‑150 (MQ=255)
                                            tgctgAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCg                                                                                                           <  2:661094/150‑1 (MQ=255)
                                            tgctgAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCg                                                                                                           >  2:151507/1‑150 (MQ=255)
                                            tgctgAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCg                                                                                                           <  1:558842/150‑1 (MQ=255)
                                                gAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTAt                                                                                                       >  2:1941969/1‑150 (MQ=255)
                                                 aaaaCCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCTGGGTGCTgat                                                                                                                                              >  1:1753328/1‑108 (MQ=255)
                                                 aaaaCCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCTGGGTGCTgat                                                                                                                                              <  2:1753328/110‑3 (MQ=255)
                                                  aaaCCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCTGGGTGCTga                                                                                                                                               >  2:1976003/1‑107 (MQ=255)
                                                  aaaCCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCTGGGTGCTga                                                                                                                                               <  1:1976003/108‑2 (MQ=255)
                                                  aaaCCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCa                                                                                                     >  2:1838975/1‑150 (MQ=255)
                                                   aaCCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCt                                                                                                                                                 >  1:144026/1‑105 (MQ=255)
                                                      cttcttTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCAcccc                                                                                                 >  1:482049/1‑150 (MQ=255)
                                                      cttcttTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCAcccc                                                                                                 >  2:232969/1‑150 (MQ=255)
                                                       ttcttTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCg                                                                                                >  1:2366997/1‑150 (MQ=255)
                                                       ttcttTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCg                                                                                                >  1:573689/1‑150 (MQ=255)
                                                       ttcttTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGAGCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCg                                                                                                >  1:956572/1‑150 (MQ=255)
                                                          tttGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTg                                                                                                                                          <  1:591755/105‑1 (MQ=255)
                                                          tttGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTg                                                                                                                                          >  2:591755/1‑105 (MQ=255)
                                                               cGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTg                                                                                                                                          <  2:507341/100‑1 (MQ=255)
                                                               cGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTg                                                                                                                                          >  1:507341/1‑100 (MQ=255)
                                                               cGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCa                                                                                                                                      >  2:2123711/1‑104 (MQ=255)
                                                                  gCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCTTCGCTATCACCCCGCAGATCGTTGc                                                                                     <  1:2130200/150‑1 (MQ=255)
                                                                      ggCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGaa                                                                                 >  2:514906/1‑150 (MQ=255)
                                                                      ggCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCCGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGaa                                                                                 <  1:2106117/150‑1 (MQ=255)
                                                                       gCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGaac                                                                                >  2:1041622/1‑150 (MQ=255)
                                                                           ggTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGCCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGAGCAACTAATGCGCCGTCACTATCACCCCGCAGATCGTTGCCGAacaaca                                                                            >  1:634446/1‑150 (MQ=255)
                                                                             tCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACACCAgc                                                                          >  2:579951/1‑149 (MQ=255)
                                                                             tCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGa                                                                          >  1:338573/1‑150 (MQ=255)
                                                                               aGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGtt                                                                                       >  1:2116079/1‑135 (MQ=255)
                                                                               aGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGtt                                                                                       <  2:2116079/135‑1 (MQ=255)
                                                                                gCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAgc                                                                       >  2:1098709/1‑150 (MQ=255)
                                                                                   aGTGGCCCAACGCCAAGCCTGCCGAAGCGGCCGCGCTGCTTGCACCGCAGCTGGGGATCCGCCCCCAGGTGCTGGAAGTGGCCAGCCAACGCCGCAACTACAACGCGGTCGCTATCCCCCCGCAGATCGTTGCCGAACAACAGAAgctgg                                                                    <  2:1124800/150‑1 (MQ=255)
                                                                                    gTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAgctggc                                                                   <  1:352981/150‑1 (MQ=255)
                                                                                    gTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAgctggc                                                                   <  2:482049/150‑1 (MQ=255)
                                                                                       ggCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGa                                                                <  2:1196368/150‑1 (MQ=255)
                                                                                        gCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGac                                                               >  2:127593/1‑150 (MQ=255)
                                                                                         ccAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGaca                                                              <  1:1838975/150‑1 (MQ=255)
                                                                                         ccAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGaca                                                              >  1:2031546/1‑150 (MQ=255)
                                                                                          cAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTAc                                                                                                                     >  1:1752036/1‑94 (MQ=255)
                                                                                          cAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTAc                                                                                                                     <  2:1752036/94‑1 (MQ=255)
                                                                                              gCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCGAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGt                                                                                                             <  2:1747108/98‑1 (MQ=255)
                                                                                              gCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGt                                                                                                             >  1:1747108/1‑98 (MQ=255)
                                                                                              gCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTAGCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGGTGACACCTTc                                                         >  1:940105/1‑150 (MQ=255)
                                                                                                cAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAg                                                                                                                                                      >  2:831175/1‑55 (MQ=255)
                                                                                                cAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAACCTGGCTGACACCTTCCa                                                       <  2:1237160/150‑1 (MQ=255)
                                                                                                   gCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAAc                                                                                                                  >  2:1561869/1‑88 (MQ=255)
                                                                                                   gCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAAc                                                                                                                  <  1:1561869/88‑1 (MQ=255)
                                                                                                     cTGCCGAAGTGGCCGCGCTGTTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAgggcc                                                  >  1:1169278/1‑150 (MQ=255)
                                                                                                      tGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTAc                                                                                                                     <  1:1633058/82‑1 (MQ=255)
                                                                                                        ccGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACGCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAgggcctgg                                               <  1:1906338/150‑1 (MQ=255)
                                                                                                        ccGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCaa                                                                                                                         <  1:2034776/76‑1 (MQ=255)
                                                                                                        ccGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCaa                                                                                                                         >  2:2034776/1‑76 (MQ=255)
                                                                                                        ccGAAGTGGCCGCGCTCCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGCAGCTGGCTGACACCTTCCAgggcctgg                                               <  1:1171224/150‑1 (MQ=255)
                                                                                                         cGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAgggcctggg                                              >  2:42821/1‑150 (MQ=255)
                                                                                                          gAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAgggcctgggc                                             <  1:1177003/150‑1 (MQ=255)
                                                                                                          gAAGGGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAgggcctgggc                                             <  1:1991174/150‑1 (MQ=255)
                                                                                                            aGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAgggcctgggcct                                           <  1:22194/150‑1 (MQ=255)
                                                                                                              tGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGa                                         >  2:763965/1‑150 (MQ=255)
                                                                                                              tGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGa                                         <  1:232969/150‑1 (MQ=255)
                                                                                                                  cgcgCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTcc                                     >  2:1595666/1‑150 (MQ=255)
                                                                                                                    cgcTGTTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCAc                                   <  2:1169278/150‑1 (MQ=255)
                                                                                                                          tgGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGa                                         >  1:732527/1‑138 (MQ=255)
                                                                                                                          tgGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGa                                         <  2:732527/138‑1 (MQ=255)
                                                                                                                            gCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTa                           >  2:435529/1‑150 (MQ=255)
                                                                                                                                 gCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGt                      >  2:1879349/1‑150 (MQ=255)
                                                                                                                                  cAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTg                     <  1:984888/150‑1 (MQ=255)
                                                                                                                                       cacttccAGCACCTGGGTGCTAGAAGTAGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAAGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGa                >  2:1727034/7‑150 (MQ=255)
                                                                                                                                        gggATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATc                                                                                          <  2:853784/75‑1 (MQ=255)
                                                                                                                                        gggATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATc                                                                                          >  1:853784/1‑75 (MQ=255)
                                                                                                                                        gggATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGaa               <  2:956572/150‑1 (MQ=255)
                                                                                                                                          ttccAGCACCTGGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTAc                          >  1:1532125/4‑137 (MQ=255)
                                                                                                                                          ttccAGCACCTGGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTAc                          <  2:1532125/134‑1 (MQ=255)
                                                                                                                                          gATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGAAGc             <  2:2366997/150‑1 (MQ=255)
                                                                                                                                           aTCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGAAGCg            <  1:2196184/150‑1 (MQ=255)
                                                                                                                                               gCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAgctgg                                                                    >  1:612830/1‑90 (MQ=255)
                                                                                                                                               gCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGAAGCGGTCt        >  1:2173429/1‑150 (MQ=255)
                                                                                                                                                   ccAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAgggcctgg                                               <  2:1422956/107‑1 (MQ=255)
                                                                                                                                                     aGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGAAGCGGTCTACCCgg  >  1:119698/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTCTACCTGTCCACAGCCAGCTTTGCCGACCAGCACCGTGCTCTGCTGAAAACCTTCTTTGCCGAGCTGGGCAAGGTCAGCCAGTGGGCCAACGCCAAGCCTGCCGAAGTGGCCGCGCTGCTGGCACCGCAGCTGGGGATCAGCACCCAGGTGCTGGAAGTGGCCAGCGAACGGCGCAACTACAACGCCGTCGCTATCACCCCGCAGATCGTTGCCGAACAACAGAAGCTGGCTGACACCTTCCAGGGCCTGGGCCTGATTCCACGCAAGCTACAGGTGGCCGAAGCGGTCTACCCGG  >  NC_002947/4000031‑4000329

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