Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 191,916:1 +G 5.2% coding (1602/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,9161.G 5.2% 1462.6 / 67.0 325coding (1602/2157 nt)paxBtoxin secretion ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (173/135);  new base G (11/6);  total (184/141)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 326 total aligned reads displayed.

GGCGCCGCTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGACCCCGCGGAACTC  >  NC_002947/191769‑192065
                                                                                                                                                    |                                                                                                                                                     
ggCGCCGCTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAg                                                                                                                                                     <  1:7189545/150‑1 (MQ=255)
 gcgcCGCTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAgg                                                                                                                                                    >  1:8026952/1‑150 (MQ=255)
 gcgcCGCTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAgg                                                                                                                                                    <  2:3431695/150‑1 (MQ=255)
  cgccgcTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGc                                                                                                                                                   >  2:1885587/1‑150 (MQ=255)
     cgcTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCcga                                                                                                                                                >  1:8602940/1‑150 (MQ=255)
      gcTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGTCACCATCGGCCTCGTACAGGCcgac                                                                                                                                               >  1:3862030/1‑150 (MQ=255)
      gcTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTGTACGAGGCcgat                                                                                                                                               <  1:6026347/150‑2 (MQ=255)
           ccAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGa                                                                                                                                          >  2:8227942/1‑150 (MQ=255)
               gTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCtt                                                                                                                                      <  1:103084/150‑1 (MQ=255)
                     aCGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAg                                                                                                                                >  2:4103103/1‑150 (MQ=255)
                      cGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGc                                                                                                                               >  2:1987158/1‑150 (MQ=255)
                       gCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCg                                                                                                                              >  1:7260146/1‑150 (MQ=255)
                           ggTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGct                                                                                                                          >  2:7177620/1‑150 (MQ=255)
                             tGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGctct                                                                                                                        >  1:877403/1‑150 (MQ=255)
                                  gCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTGGCCg                                                                                                                   >  2:7032713/1‑150 (MQ=255)
                                    cagcagCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGa                                                                                                                 >  2:134968/1‑150 (MQ=255)
                                    cagcagCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGa                                                                                                                 <  1:1490979/150‑1 (MQ=255)
                                     agcagcTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAg                                                                                                                <  1:2765531/150‑1 (MQ=255)
                                     agcagcTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAg                                                                                                                <  1:1118437/150‑1 (MQ=255)
                                     agcagcTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAg                                                                                                                <  1:7152109/150‑1 (MQ=255)
                                         gcTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGc                                                                                                            <  2:917379/150‑1 (MQ=255)
                                          cTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCt                                                                                                           >  2:1069156/1‑150 (MQ=255)
                                             gATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCg                                                                                                        >  2:6629012/1‑150 (MQ=255)
                                             gATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCg                                                                                                        >  2:6940158/1‑150 (MQ=255)
                                              atatCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCgg                                                                                                       <  1:5074991/150‑1 (MQ=255)
                                               tatCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGc                                                                                                      >  2:7371208/1‑150 (MQ=255)
                                                 tCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCg                                                                                                    >  1:5737627/1‑150 (MQ=255)
                                                 tCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCg                                                                                                    >  2:7453457/1‑150 (MQ=255)
                                                  ccTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCga                                                                                                   <  2:7622739/150‑1 (MQ=255)
                                                    tGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCgatg                                                                                                 >  1:1343906/1‑150 (MQ=255)
                                                    tGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCgatg                                                                                                 >  1:5944459/1‑150 (MQ=255)
                                                        ggACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGc                                                                                             >  1:1530383/1‑150 (MQ=255)
                                                              aGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCt                                                                                       >  2:3827948/1‑150 (MQ=255)
                                                                 cGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTtct                                                                                    >  1:4092870/1‑150 (MQ=255)
                                                                  gATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTtctc                                                                                   >  2:4167955/1‑150 (MQ=255)
                                                                  gATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTtctc                                                                                   <  2:6246105/150‑1 (MQ=255)
                                                                       tgtgGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCCCCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAg                                                                              >  1:8366817/1‑150 (MQ=255)
                                                                         tgGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAggg                                                                            >  1:4394047/1‑150 (MQ=255)
                                                                           gcgcAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCg                                                                          >  2:50303/1‑150 (MQ=255)
                                                                             gcAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGa                                                                        >  2:8671213/1‑150 (MQ=255)
                                                                              cAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGAt                                                                       <  2:6868510/150‑1 (MQ=255)
                                                                               aGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATg                                                                      <  2:7887468/150‑1 (MQ=255)
                                                                                gTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATgg                                                                     <  1:7083792/150‑1 (MQ=255)
                                                                                    gCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAg                                                                 >  1:7708536/1‑150 (MQ=255)
                                                                                    gCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAg                                                                 <  2:7915157/150‑1 (MQ=255)
                                                                                    gCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAg                                                                 <  2:5659300/150‑1 (MQ=255)
                                                                                         cATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGa                                                            <  2:4442456/150‑1 (MQ=255)
                                                                                         cATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGa                                                            <  2:9044808/150‑1 (MQ=255)
                                                                                          aTCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGat                                                           <  1:1765179/150‑1 (MQ=255)
                                                                                          aTCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGat                                                           <  2:3753220/150‑1 (MQ=255)
                                                                                             gATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATAtt                                                        >  1:5201417/1‑150 (MQ=255)
                                                                                              aTCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATAttc                                                       <  2:542435/150‑1 (MQ=255)
                                                                                              aTCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATAttc                                                       <  1:5483634/150‑1 (MQ=255)
                                                                                              aTCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATAttc                                                       <  1:3327832/150‑1 (MQ=255)
                                                                                                 tGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATAttcttc                                                    >  2:2756784/1‑150 (MQ=255)
                                                                                                 tGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATAttcttc                                                    <  2:3229764/150‑1 (MQ=255)
                                                                                                  ggCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCa                                                   <  1:9072142/150‑1 (MQ=255)
                                                                                                  ggCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCa                                                   <  1:6463862/150‑1 (MQ=255)
                                                                                                   gCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAg                                                  >  1:3588268/1‑150 (MQ=255)
                                                                                                    cGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAgg                                                 <  1:2231139/150‑1 (MQ=255)
                                                                                                      aaTGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGc                                               <  1:6227186/150‑1 (MQ=255)
                                                                                                          tCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAgg                                           >  1:7394271/1‑150 (MQ=255)
                                                                                                           cGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGt                                          >  1:3966493/1‑150 (MQ=255)
                                                                                                             aCGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTc                                        <  2:117400/150‑1 (MQ=255)
                                                                                                             aCGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTc                                        >  2:528354/1‑150 (MQ=255)
                                                                                                             aCGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTc                                        <  2:117242/150‑1 (MQ=255)
                                                                                                             aCGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTc                                        <  1:2716079/150‑1 (MQ=255)
                                                                                                               gCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTctg                                      >  1:6601222/1‑150 (MQ=255)
                                                                                                               gCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTctg                                      <  2:8026952/150‑1 (MQ=255)
                                                                                                                ccATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTctgc                                     <  1:184810/150‑1 (MQ=255)
                                                                                                                 cATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGAc                                                                                         >  2:2284842/1‑97 (MQ=255)
                                                                                                                 cATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTctgct                                    <  1:8121410/150‑1 (MQ=255)
                                                                                                                  aTCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTctgctg                                   <  1:2305815/150‑1 (MQ=255)
                                                                                                                   tCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGcc                                                                                                                    >  1:7676393/1‑68 (MQ=255)
                                                                                                                   tCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGcc                                                                                                                    <  2:7676393/68‑1 (MQ=255)
                                                                                                                      aCCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTc                               <  1:559003/150‑1 (MQ=255)
                                                                                                                       ccAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCg                              >  2:1536980/1‑150 (MQ=255)
                                                                                                                       ccAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCg                              <  2:5960020/150‑1 (MQ=255)
                                                                                                                        cagcagCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCgg                             >  2:1427416/1‑150 (MQ=255)
                                                                                                                        cagcagCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCgg                             >  2:3129855/1‑150 (MQ=255)
                                                                                                                         agcagcGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGa                            >  1:3374430/1‑150 (MQ=255)
                                                                                                                             gcGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAg                        >  2:2990036/1‑150 (MQ=255)
                                                                                                                              cGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGt                       <  2:5559222/150‑1 (MQ=255)
                                                                                                                               gAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTg                      >  2:438021/1‑150 (MQ=255)
                                                                                                                                 gCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGaa                    <  1:7880853/150‑1 (MQ=255)
                                                                                                                                 gCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGaa                    <  2:4654284/150‑1 (MQ=255)
                                                                                                                                 gCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGaa                    <  1:8227942/150‑1 (MQ=255)
                                                                                                                                  tcatcgtCGGCCTGTACGAGGCCGATGAAGAGCTTGGCCAGAGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCTCCAGGGCGGATGGTCAGGATGATATTCTTCAGGGCCAGGTTCTGCAGGTACGGATAAGTGAAg                   >  1:2439645/7‑150 (MQ=255)
                                                                                                                                   caccaTCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGt                  >  1:4537561/1‑150 (MQ=255)
                                                                                                                                    accaTCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTc                 <  2:1752972/150‑1 (MQ=255)
                                                                                                                                    accaTCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTc                 >  2:2554935/1‑150 (MQ=255)
                                                                                                                                      cgtCGGCCTGTACGAGGCCGATGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGa               >  2:4788980/3‑150 (MQ=39)
                                                                                                                                      caTCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGa               >  1:2513561/1‑150 (MQ=255)
                                                                                                                                      caTCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGa               <  2:3519153/150‑1 (MQ=255)
                                                                                                                                       gtCGGCCTGTACGAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGAc              >  1:6831114/2‑150 (MQ=255)
                                                                                                                                          ggCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGAcccc           >  2:5159152/1‑150 (MQ=255)
                                                                                                                                              tCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGACCCCGCgg       >  1:8543362/1‑150 (MQ=255)
                                                                                                                                              tCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGACCCCGCgg       >  1:1901453/1‑150 (MQ=255)
                                                                                                                                              tCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGACCCCGCgg       >  1:6676294/1‑150 (MQ=255)
                                                                                                                                                gTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGACCCCGCGGaa     >  1:471402/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCGCCGCTGACCAGGTTGTCACGCAAGGTGCCGGCCAGCAGCTGGATATCCTGGGGGACGTAGCCGATGTTGTGGCGCAGTTCGCTGACATCGATCTGGCGAATGTCGACGCCATCAACCAGCAGCGAGCCACCATCGGCCTCGTACAGGCCGACGATGAGCTTGGCCAGCGAGCTCTTGCCGGAGCCGCTGCGGCCGATGATGCCGACCTTCTCGCCAGGGCGGATGGTCAGGTTGATATTCTTCAGGGCCAGGTTCTGCTGGTTCGGATAAGTGAAGTCGACCCCGCGGAACTC  >  NC_002947/191769‑192065

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: