Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,706,232 C→G 5.1% intergenic (+25/‑38) phaA → / → phaB poly(3‑hydroxyalkanoate) polymerase/poly(3‑hydroxyalkanoate) depolymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,706,2320CG 5.1% 1199.5 / 37.4 351intergenic (+25/‑38)phaA/phaBpoly(3‑hydroxyalkanoate) polymerase/poly(3‑hydroxyalkanoate) depolymerase
Reads supporting (aligned to +/- strand):  ref base C (147/186);  new base G (14/4);  total (161/190)
Fisher's exact test for biased strand distribution p-value = 6.63e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 351 total aligned reads displayed.

TTGCACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATCTTTAACGGCATC  >  NC_002947/5706085‑5706380
                                                                                                                                                   |                                                                                                                                                    
ttGCACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACgg                                                                                                                                                    >  2:2146485/1‑150 (MQ=255)
ttGCACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACgg                                                                                                                                                    >  1:380752/1‑150 (MQ=255)
    aCTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGcc                                                                                                                                                >  2:3540565/1‑150 (MQ=255)
    aCTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGCCGTGCAGgccgcc                                                                                                                                                <  2:245583/150‑1 (MQ=255)
     cTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCa                                                                                                                                               <  1:5641837/150‑1 (MQ=255)
       ggCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACg                                                                                                                                             >  1:8278465/1‑150 (MQ=255)
         cAGAGTTGGCTGGGCGAGCGTCCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGgt                                                                                                                                           <  1:3463983/150‑1 (MQ=255)
           gagTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCACCGTGGCGCCGTGCAGGCCGCCACggcgc                                                                                                                                         <  2:5973581/150‑4 (MQ=255)
            agTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGtt                                                                                                                                        <  1:8607745/150‑1 (MQ=255)
            agTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGtt                                                                                                                                        <  1:922314/150‑1 (MQ=255)
            agTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGtt                                                                                                                                        <  2:3638559/150‑1 (MQ=255)
              ttGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTa                                                                                                                                      <  1:6294504/150‑1 (MQ=255)
              tgggCTGGGCGAGCGTGCCGGCGCGCTGAAACAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTa                                                                                                                                      <  2:5420472/148‑1 (MQ=255)
                 gCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACtt                                                                                                                                   <  2:5125522/150‑1 (MQ=255)
                 gCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACtt                                                                                                                                   <  2:1713188/150‑1 (MQ=255)
                 gCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACtt                                                                                                                                   <  2:2615523/150‑1 (MQ=255)
                   tGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCa                                                                                                                                 <  1:3256319/150‑1 (MQ=255)
                   tGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCa                                                                                                                                 >  1:6363165/1‑150 (MQ=255)
                   tGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCa                                                                                                                                 >  1:5103723/1‑150 (MQ=255)
                    gggCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTGTACTTCAc                                                                                                                                >  2:474566/1‑150 (MQ=255)
                     ggCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCAcc                                                                                                                               >  2:1921065/1‑150 (MQ=255)
                      gCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCAccc                                                                                                                              >  2:6878772/1‑150 (MQ=255)
                      gCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCAccc                                                                                                                              >  1:7562069/1‑150 (MQ=255)
                       cGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTCTGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCa                                                                                                                             >  2:7675749/1‑150 (MQ=255)
                        gAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCaa                                                                                                                            <  2:4153048/150‑1 (MQ=255)
                          gCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAaga                                                                                                                          <  2:6814548/150‑1 (MQ=255)
                          gCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAaga                                                                                                                          <  2:5784273/150‑1 (MQ=255)
                               ccGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCAc                                                                                                                     <  1:6738555/150‑1 (MQ=255)
                               ccGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCAc                                                                                                                     <  2:8135111/150‑1 (MQ=255)
                                  gcgcgcTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTg                                                                                                                  <  2:3423638/150‑1 (MQ=255)
                                       cTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGc                                                                                                             >  2:5005038/1‑150 (MQ=255)
                                       cTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGc                                                                                                             >  1:6641926/1‑150 (MQ=255)
                                        tGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCAGGTGCATGcc                                                                                                            >  1:3133546/1‑150 (MQ=255)
                                         gAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCg                                                                                                           >  2:220088/1‑150 (MQ=255)
                                          aaaaaGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGc                                                                                                          >  1:4590133/1‑150 (MQ=255)
                                          aaaaaGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGc                                                                                                          >  2:444174/1‑150 (MQ=255)
                                            aaaGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCaa                                                                                                        >  2:2298067/1‑150 (MQ=255)
                                              aGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAAcc                                                                                                      <  1:4770379/150‑1 (MQ=255)
                                                gCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCt                                                                                                    >  1:7764747/1‑150 (MQ=255)
                                                gCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCt                                                                                                    >  1:4079174/1‑150 (MQ=255)
                                                  aCCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCtat                                                                                                  >  1:592168/1‑150 (MQ=255)
                                                    cGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTatat                                                                                                >  2:6190359/1‑150 (MQ=255)
                                                        ccGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCCTGCCGCAACCCTATATTTTc                                                                                            >  1:3679370/1‑150 (MQ=255)
                                                        ccGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTc                                                                                            >  1:7198453/1‑150 (MQ=255)
                                                         cGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCa                                                                                           <  2:3046965/150‑1 (MQ=255)
                                                         cGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCa                                                                                           <  2:4291137/150‑1 (MQ=255)
                                                         cGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCa                                                                                           <  1:1196498/150‑1 (MQ=255)
                                                          gCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAg                                                                                          <  1:2351854/150‑1 (MQ=255)
                                                          gCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTACAg                                                                                          >  2:4905292/1‑150 (MQ=255)
                                                             tGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGAc                                                                                       >  1:2635205/1‑150 (MQ=255)
                                                             tGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGAc                                                                                       >  1:7924215/1‑150 (MQ=255)
                                                                 cAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTc                                                                                   <  1:434003/150‑1 (MQ=255)
                                                                      gTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCt                                                                              <  1:2750066/150‑1 (MQ=255)
                                                                      gTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCt                                                                              <  1:7649564/150‑1 (MQ=255)
                                                                      gTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCt                                                                              <  2:2953339/150‑1 (MQ=255)
                                                                       tACCTATGCCGCCGGCGTAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTg                                                                             <  2:3300874/150‑1 (MQ=255)
                                                                          cTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGAc                                                                          <  2:7389889/150‑1 (MQ=255)
                                                                           tATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACa                                                                         >  2:2384502/1‑150 (MQ=255)
                                                                               ccgccgGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCa                                                                     <  2:3730150/150‑1 (MQ=255)
                                                                               ccgccgGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCa                                                                     <  1:249448/150‑1 (MQ=255)
                                                                                cgccgGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAg                                                                    >  2:5642011/1‑150 (MQ=255)
                                                                                cgccgGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAg                                                                    >  2:6618632/1‑150 (MQ=255)
                                                                                 gccgGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGt                                                                   <  1:6223816/150‑1 (MQ=255)
                                                                                  ccgGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTc                                                                  >  1:3354354/1‑150 (MQ=255)
                                                                                     gCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCAt                                                               >  2:3894715/1‑150 (MQ=255)
                                                                                      cGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATc                                                              <  2:6769277/150‑1 (MQ=255)
                                                                                       gAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATcc                                                             <  2:8164775/150‑1 (MQ=255)
                                                                                               cccAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCt                                                     <  1:2129503/150‑1 (MQ=255)
                                                                                               cccAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCt                                                     <  1:291371/150‑1 (MQ=255)
                                                                                                   ggCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTc                                                 <  1:1835650/150‑1 (MQ=255)
                                                                                                     cACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGc                                               <  1:5351040/150‑1 (MQ=255)
                                                                                                       ccTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCAc                                                         >  1:6849115/1‑138 (MQ=255)
                                                                                                        cTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGcccc                                            >  1:8413641/1‑150 (MQ=255)
                                                                                                        cTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGcccc                                            <  2:8323006/150‑1 (MQ=255)
                                                                                                        cTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGcccc                                            >  1:6153933/1‑150 (MQ=255)
                                                                                                         tACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCg                                           <  2:2864485/150‑1 (MQ=255)
                                                                                                              tCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAg                                      <  1:492256/150‑1 (MQ=255)
                                                                                                                aCGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGcc                                    >  2:6493017/1‑150 (MQ=255)
                                                                                                                  gAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGc                                  >  2:4203267/1‑150 (MQ=255)
                                                                                                                    gCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCAc                                <  1:223118/150‑1 (MQ=255)
                                                                                                                       ttGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTg                             <  1:1088739/150‑1 (MQ=255)
                                                                                                                            ttACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGcc                        >  1:875867/1‑150 (MQ=255)
                                                                                                                             gtaaacaccGTGGCGCCGTGCAGGCCGCCACGGCGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCg                       >  2:8723869/8‑150 (MQ=255)
                                                                                                                              aCAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGt                      >  1:7333953/1‑150 (MQ=255)
                                                                                                                                aacaccGTGGCGCCGTGCAGGCCGCCACGGCGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTg                    >  1:1447160/5‑150 (MQ=255)
                                                                                                                                 acaccGTGGCGCCGTGCAGGCCGCCACGGCGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGc                   >  1:7468962/4‑150 (MQ=255)
                                                                                                                                   accGTGGCGCCGTGCAGGCCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTg                 >  2:8221541/2‑150 (MQ=255)
                                                                                                                                   accGTGGCGCCGTGCAGGCCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTg                 >  2:1281481/2‑150 (MQ=255)
                                                                                                                                   accGTGGCGCCGTGCAGGCCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTg                 >  2:1145647/2‑150 (MQ=255)
                                                                                                                                    ccGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGa                <  2:8413641/150‑1 (MQ=255)
                                                                                                                                      gTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATc              <  1:8163627/150‑1 (MQ=255)
                                                                                                                                       tGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATCt             <  2:2438749/150‑1 (MQ=255)
                                                                                                                                         gcggcCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCACCCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGGGGCTGATCttt           >  2:6112186/1‑150 (MQ=255)
                                                                                                                                           ggcCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATCTTTaa         <  1:7436105/150‑1 (MQ=255)
                                                                                                                                           ggcCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATCTTTaa         <  1:4246482/150‑1 (MQ=255)
                                                                                                                                           ggcCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCATATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATCTTTaa         <  2:7863814/150‑1 (MQ=255)
                                                                                                                                            gcCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATCTTTAAc        <  1:598283/150‑1 (MQ=255)
                                                                                                                                               tGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATCTTTAACGGc     <  2:1425693/150‑1 (MQ=255)
                                                                                                                                                gCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATCTTTAACGGCa    >  1:695248/1‑150 (MQ=255)
                                                                                                                                                  aCGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATCTTTAACGGCATc  <  2:6363165/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTGCACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATCTTTAACGGCATC  >  NC_002947/5706085‑5706380

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: