Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 191,912 Δ1 bp 5.3% coding (1606/2157 nt) paxB ← toxin secretion ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947191,9120C. 5.3% 1098.8 / 52.7 323coding (1606/2157 nt)paxBtoxin secretion ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (171/135);  new base . (11/6);  total (182/141)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 326 total aligned reads displayed.

GTAGCGGGCGCCGCTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGACCCCGCGG  >  NC_002947/191763‑192060
                                                                                                                                                     |                                                                                                                                                     
   gCGGGCGCCGCTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTa                                                                                                                                                    <  2:8792664/150‑1 (MQ=255)
   gCGGGCGCCGCTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTa                                                                                                                                                    >  2:842833/1‑150 (MQ=255)
      ggCGCCGCTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAg                                                                                                                                                <  1:7189545/150‑1 (MQ=255)
            gcTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTGTACGAGGCcgat                                                                                                                                          <  2:3361546/150‑2 (MQ=255)
             cTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCcgacg                                                                                                                                         >  1:3992489/1‑150 (MQ=255)
              tGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCcgacga                                                                                                                                        >  2:7706803/1‑150 (MQ=255)
                  cAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCGTCGGCCTGTACGAGGCcgatggtggc                                                                                                                                    <  2:2870970/150‑8 (MQ=255)
                  cAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAg                                                                                                                                    >  1:7912856/1‑150 (MQ=255)
                  cAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAg                                                                                                                                    >  2:7028158/1‑150 (MQ=255)
                   aGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGc                                                                                                                                   <  2:1363498/150‑1 (MQ=255)
                     gTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCtt                                                                                                                                 <  1:512953/150‑1 (MQ=255)
                     gTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCtt                                                                                                                                 <  1:103084/150‑1 (MQ=255)
                      ttGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTg                                                                                                                                <  1:3925159/150‑1 (MQ=255)
                           aCGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAg                                                                                                                           >  1:482735/1‑150 (MQ=255)
                             gCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCg                                                                                                                         >  2:4730711/1‑150 (MQ=255)
                             gCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCg                                                                                                                         >  1:7260146/1‑150 (MQ=255)
                             gCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCg                                                                                                                         >  2:125256/1‑150 (MQ=255)
                               aaGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAg                                                                                                                       >  2:3721792/1‑150 (MQ=255)
                                aGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGc                                                                                                                      <  2:761776/150‑1 (MQ=255)
                                 ggTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGct                                                                                                                     >  2:7177453/1‑150 (MQ=255)
                                    gCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCtt                                                                                                                  >  2:559310/1‑150 (MQ=255)
                                     ccGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTg                                                                                                                 >  2:1845191/1‑150 (MQ=255)
                                          cagcagCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGa                                                                                                            <  1:1490979/150‑1 (MQ=255)
                                          cagcagCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGa                                                                                                            <  2:6932405/150‑1 (MQ=255)
                                          cagcagCTGGATATACTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCCCTGACATCGATCTGGCGAATGTCGACGAAATCAACCAGCAGCGAGCCACAATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGa                                                                                                            >  2:3103265/1‑150 (MQ=255)
                                           agcagcTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAg                                                                                                           <  1:7152109/150‑1 (MQ=255)
                                           agcagcTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAg                                                                                                           <  1:4396148/150‑1 (MQ=255)
                                           agcagcTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAg                                                                                                           <  1:2765531/150‑1 (MQ=255)
                                           agcagcTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAg                                                                                                           <  1:1118437/150‑1 (MQ=255)
                                              agcTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCg                                                                                                        <  2:7952556/150‑1 (MQ=255)
                                               gcTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGc                                                                                                       <  2:917379/150‑1 (MQ=255)
                                                  ggATATCCTGGGGGACGGAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGc                                                                                                    <  2:5207339/150‑1 (MQ=255)
                                                   gATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCg                                                                                                   >  2:7448978/1‑150 (MQ=255)
                                                   gATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCg                                                                                                   <  2:6147592/150‑1 (MQ=255)
                                                       tCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCg                                                                                               >  1:5737627/1‑150 (MQ=255)
                                                        ccTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCga                                                                                              >  2:1532805/1‑150 (MQ=255)
                                                          tGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCgatg                                                                                            >  1:1343906/1‑150 (MQ=255)
                                                           gggggACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCgatga                                                                                           <  2:6660169/150‑1 (MQ=255)
                                                            ggggACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCgatgat                                                                                          <  2:2710794/150‑1 (MQ=255)
                                                              ggACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGc                                                                                        >  1:1530383/1‑150 (MQ=255)
                                                                  gTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGAc                                                                                    <  2:3757831/150‑1 (MQ=255)
                                                                  gTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGAc                                                                                    >  2:2052350/1‑150 (MQ=255)
                                                                     gCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCtt                                                                                 >  1:3646097/1‑150 (MQ=255)
                                                                       cGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGACGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTtct                                                                               >  1:4092933/1‑150 (MQ=255)
                                                                        gATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTtctc                                                                              >  1:8213959/1‑150 (MQ=255)
                                                                             tcagccgcAGATCGCAGACATCGAACTGGCGAATGTCGACGCCAACAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCACGCCAg                                                                         >  1:8366899/6‑150 (MQ=255)
                                                                                   gcAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGa                                                                   >  2:1975914/1‑150 (MQ=255)
                                                                                   gcAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGa                                                                   >  2:5135823/1‑150 (MQ=255)
                                                                                    cAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGAt                                                                  <  2:3898396/150‑1 (MQ=255)
                                                                                    cAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGAt                                                                  <  1:8832754/150‑1 (MQ=255)
                                                                                     aGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATg                                                                 <  2:7887469/150‑1 (MQ=255)
                                                                                     aGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATg                                                                 <  2:586046/150‑1 (MQ=255)
                                                                                     aGTTCGCGGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATg                                                                 <  1:6404998/150‑1 (MQ=255)
                                                                                      gTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATgg                                                                <  1:7083792/150‑1 (MQ=255)
                                                                                            tGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGt                                                          <  2:439872/150‑1 (MQ=255)
                                                                                            tGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGt                                                          >  1:3390243/1‑150 (MQ=255)
                                                                                              aCATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTg                                                        <  1:6671640/150‑1 (MQ=255)
                                                                                               cATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGa                                                       <  2:9044808/150‑1 (MQ=255)
                                                                                                aTCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGat                                                      <  1:1765179/150‑1 (MQ=255)
                                                                                                   gATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATAtt                                                   >  1:5883852/1‑150 (MQ=255)
                                                                                                    aTCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATAttc                                                  <  1:3328316/150‑1 (MQ=255)
                                                                                                    aTCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATAttc                                                  <  1:5483634/150‑1 (MQ=255)
                                                                                                        ggCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCa                                              <  2:4112708/150‑1 (MQ=255)
                                                                                                        ggCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCa                                              <  1:4649677/150‑1 (MQ=255)
                                                                                                        ggCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCa                                              <  2:5421506/150‑1 (MQ=255)
                                                                                                         gCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAg                                             <  2:1076063/150‑1 (MQ=255)
                                                                                                          cGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAgg                                            <  1:2231139/150‑1 (MQ=255)
                                                                                                          cGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAgg                                            <  1:6688092/150‑1 (MQ=255)
                                                                                                            aaTGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGc                                          <  1:6226896/150‑1 (MQ=255)
                                                                                                                tCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAgg                                      >  1:7394271/1‑150 (MQ=255)
                                                                                                                 cGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGt                                     >  2:4991212/1‑150 (MQ=255)
                                                                                                                   aCGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTc                                   <  1:2716079/150‑1 (MQ=255)
                                                                                                                     gCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTctg                                 >  1:6601210/1‑150 (MQ=255)
                                                                                                                      ccATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTctgc                                <  1:184810/150‑1 (MQ=255)
                                                                                                                       cATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTctgct                               >  2:1172574/1‑150 (MQ=255)
                                                                                                                       cATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTctgct                               >  2:657390/1‑150 (MQ=255)
                                                                                                                        aTCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTctgctg                              <  1:2305815/150‑1 (MQ=255)
                                                                                                                         tCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTgg                             >  1:6009740/1‑150 (MQ=255)
                                                                                                                            aCCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTc                          <  1:7706803/150‑1 (MQ=255)
                                                                                                                            aCCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTc                          <  1:559003/150‑1 (MQ=255)
                                                                                                                              cagcagCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCGTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTGCTCGCCAGGGCGGATGTCAGGGTGAGCTTCTTCAGGGCCAGGTTCTGCGGGGGCGGa                       >  2:7659569/1‑150 (MQ=255)
                                                                                                                                 cagcGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATa                     >  1:9163061/1‑150 (MQ=255)
                                                                                                                                  agcGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATaa                    >  1:8607475/1‑150 (MQ=255)
                                                                                                                                    cGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGt                  <  2:5559147/150‑1 (MQ=255)
                                                                                                                                     gAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTg                 >  2:1675260/1‑150 (MQ=255)
                                                                                                                                        tcatcgtCGGCCTGTACGAGGCCGATGAAGAGCTTGGCCAGAGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCTCCAGGGCGGATGGTCAGGATGATATTCTTCAGGGCCAGGTTCTGCAGGTACGGATAAGTGAAg              >  1:2439645/7‑150 (MQ=255)
                                                                                                                                        ccaccaTCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAg              >  2:3172274/1‑150 (MQ=255)
                                                                                                                                        ccaccaTCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAg              <  1:7542069/150‑1 (MQ=255)
                                                                                                                                        ccaccaTCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAg              <  1:7760467/150‑1 (MQ=255)
                                                                                                                                          accaTCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTc            >  2:2464476/1‑150 (MQ=255)
                                                                                                                                           tcgtCGGCCTGTACGAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCg           >  1:8117849/4‑150 (MQ=255)
                                                                                                                                           tcgtCGGCCTGTACGAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCg           >  1:5319378/4‑150 (MQ=255)
                                                                                                                                            cgtCGGCCTGTACGAGGCCGATGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGa          >  2:4267295/3‑150 (MQ=39)
                                                                                                                                            caTCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGa          >  1:2513561/1‑150 (MQ=255)
                                                                                                                                              tCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGAcc        >  2:457438/1‑150 (MQ=255)
                                                                                                                                               cGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGAccc       <  2:7295513/150‑1 (MQ=255)
                                                                                                                                                 gCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGACCCcg     >  2:8067713/1‑150 (MQ=255)
                                                                                                                                                    tCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGACCCCGCgg  >  1:1901453/1‑150 (MQ=255)
                                                                                                                                                    tCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGACCCCGCgg  >  2:3482356/1‑150 (MQ=255)
                                                                                                                                                    tCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGACCCCGCgg  >  1:8543372/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTAGCGGGCGCCGCTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGACCCCGCGG  >  NC_002947/191763‑192060

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