Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,703,200 C→G 5.6% intergenic (+20/‑43) 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,703,2000CG 5.6% 670.9 / 22.8 196intergenic (+20/‑43)phaA/phaBpoly(3‑hydroxyalkanoate) polymerase/poly(3‑hydroxyalkanoate) depolymerase
Reads supporting (aligned to +/- strand):  ref base C (85/100);  new base G (4/7);  total (89/107)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 197 total aligned reads displayed.

TGGTGGTTGCACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATCTTTAACGG  >  NC_002947/5703052‑5703349
                                                                                                                                                    |                                                                                                                                                     
tggtggTTGCACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCt                                                                                                                                                      <  1:518566/150‑1 (MQ=255)
 ggtggtTGCACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTg                                                                                                                                                     <  1:1060008/150‑1 (MQ=255)
  gtggtTGCACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGc                                                                                                                                                    >  1:83380/1‑150 (MQ=255)
  gtggtTGCACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGc                                                                                                                                                    >  1:2701264/1‑150 (MQ=255)
    ggtTGCACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCAc                                                                                                                                                  >  1:1490961/1‑150 (MQ=255)
      ttGCACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACgg                                                                                                                                                <  1:1304654/150‑1 (MQ=255)
       tGCACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGgc                                                                                                                                               <  1:3944222/150‑1 (MQ=255)
        gCACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGgcg                                                                                                                                              >  1:4453011/1‑150 (MQ=255)
        gCACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGgcg                                                                                                                                              <  2:1335440/150‑1 (MQ=255)
          aCTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGCCGTGCAGgccgcc                                                                                                                                            <  1:2622307/150‑1 (MQ=255)
            tGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGCCGTGCAGGCCGCCAc                                                                                                                                          <  2:2607798/150‑1 (MQ=255)
            tGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGCCGTGCAGGCCGCCAc                                                                                                                                          <  2:349656/150‑1 (MQ=255)
              gCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACgg                                                                                                                                        >  2:522343/1‑150 (MQ=255)
              gCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACgg                                                                                                                                        >  1:238814/1‑150 (MQ=255)
                agagTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGgtg                                                                                                                                      <  1:2842974/150‑1 (MQ=255)
                 gagTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGgtgt                                                                                                                                     <  2:4287892/150‑1 (MQ=255)
                  aggtGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCGGGCGAAGCCTCCCCCGGCACCTACGGTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGtt                                                                                                                                    <  2:819751/147‑1 (MQ=255)
                  agTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGtt                                                                                                                                    <  2:534905/150‑1 (MQ=255)
                   gTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGttt                                                                                                                                   >  2:3318097/1‑150 (MQ=255)
                       gCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACtt                                                                                                                               >  2:4337374/1‑150 (MQ=255)
                          gggCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCAc                                                                                                                            <  1:3682957/150‑1 (MQ=255)
                           ggCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCAcc                                                                                                                           >  1:3743988/1‑150 (MQ=255)
                           ggCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCAcc                                                                                                                           <  2:412455/150‑1 (MQ=255)
                            gCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCAccc                                                                                                                          >  1:459761/1‑150 (MQ=255)
                            gCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCAccc                                                                                                                          <  2:2922881/150‑1 (MQ=255)
                            gCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCAccc                                                                                                                          <  2:4324076/150‑1 (MQ=255)
                             cGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCa                                                                                                                         >  1:3696719/1‑150 (MQ=255)
                                gCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAaga                                                                                                                      <  1:2319120/150‑1 (MQ=255)
                                gCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAaga                                                                                                                      <  1:1515643/150‑1 (MQ=255)
                                   tGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTc                                                                                                                   >  1:3270606/1‑150 (MQ=255)
                                       ggCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGTGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGt                                                                                                               <  1:3233331/150‑1 (MQ=255)
                                       ggCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGt                                                                                                               <  1:1136846/150‑1 (MQ=255)
                                       ggCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGt                                                                                                               <  1:2263839/150‑1 (MQ=255)
                                       ggCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGt                                                                                                               <  2:3200789/150‑1 (MQ=255)
                                          gcgcTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCa                                                                                                            <  2:1244827/150‑1 (MQ=255)
                                          gcgcTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCa                                                                                                            <  1:3032361/150‑1 (MQ=255)
                                            gcTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATg                                                                                                          >  1:1640122/1‑150 (MQ=255)
                                              tGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGcc                                                                                                        >  2:2575289/1‑150 (MQ=255)
                                               gAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCg                                                                                                       >  1:1357066/1‑150 (MQ=255)
                                               gAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCg                                                                                                       >  1:2081332/1‑150 (MQ=255)
                                                     ggCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAAccc                                                                                                 <  2:2243198/150‑1 (MQ=255)
                                                      gCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCt                                                                                                <  1:514914/150‑1 (MQ=255)
                                                      gCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCt                                                                                                >  1:4245992/1‑150 (MQ=255)
                                                          cGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTatat                                                                                            <  1:1376588/150‑1 (MQ=255)
                                                           gACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATAtt                                                                                           >  2:4218325/1‑150 (MQ=255)
                                                              ccGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTAc                                                                                                                                 >  2:406528/1‑109 (MQ=255)
                                                              ccGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTAc                                                                                                                                 <  1:406528/109‑1 (MQ=255)
                                                              ccGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTc                                                                                        >  1:4352759/1‑150 (MQ=255)
                                                              ccGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTc                                                                                        >  2:3939809/1‑150 (MQ=255)
                                                              ccGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGGGTTTACTTCACCCAAGAGTCACGTGCATGCCGCCACCCTATATTTTc                                                                                        >  2:3119616/1‑150 (MQ=255)
                                                               cGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCa                                                                                       <  1:362359/150‑1 (MQ=255)
                                                               cGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCa                                                                                       <  2:238814/150‑1 (MQ=255)
                                                               cGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCa                                                                                       <  1:2230937/150‑1 (MQ=255)
                                                               cGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCa                                                                                       <  2:1194844/150‑1 (MQ=255)
                                                               cGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCa                                                                                       <  2:1194917/150‑1 (MQ=255)
                                                               cGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCa                                                                                       <  1:2457937/150‑1 (MQ=255)
                                                               cGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCACCGTGGCGCCGTGCAGGCCGCCACggcgct                                                                                                                                    >  2:2205320/1‑101 (MQ=37)
                                                                gCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAg                                                                                      <  2:873233/150‑1 (MQ=255)
                                                                  cTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGa                                                                                    <  2:1903948/150‑1 (MQ=255)
                                                                   tGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGAc                                                                                   >  2:3260938/1‑150 (MQ=255)
                                                                    gggCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGAcc                                                                                  <  1:3110159/150‑1 (MQ=255)
                                                                     ggCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCg                                                                                 >  1:2681637/1‑150 (MQ=255)
                                                                      gCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAaga                                                                                                                      >  2:1088637/1‑112 (MQ=255)
                                                                         aCCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGa                                                                             <  2:4181894/150‑1 (MQ=255)
                                                                           cGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTGTATTTTCAGGACCGTCGAGc                                                                           >  2:2861500/1‑150 (MQ=255)
                                                                           cGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGc                                                                           >  2:2570479/1‑150 (MQ=255)
                                                                                cTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGAc                                                                      >  1:557470/1‑150 (MQ=255)
                                                                                cTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGAc                                                                      <  1:100823/150‑1 (MQ=255)
                                                                                   tGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAAc                                                                   <  2:1751223/150‑1 (MQ=255)
                                                                                   tGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAAc                                                                   <  1:665774/150‑1 (MQ=255)
                                                                                     ccgccgGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCa                                                                 <  2:3737219/150‑1 (MQ=255)
                                                                                     ccgccgGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCa                                                                 <  2:3737191/150‑1 (MQ=255)
                                                                                        ccgGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTc                                                              >  1:984362/1‑150 (MQ=255)
                                                                                              aaGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCg                                                        <  2:3782161/150‑1 (MQ=255)
                                                                                                   tCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCg                                                   >  1:1810646/1‑150 (MQ=255)
                                                                                                   tCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGGGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCg                                                   >  1:1242744/1‑150 (MQ=255)
                                                                                                      ccAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTg                                                >  2:761617/1‑150 (MQ=255)
                                                                                                      ccAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTg                                                >  1:2157857/1‑150 (MQ=255)
                                                                                                       cAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGt                                               <  1:3541640/150‑1 (MQ=255)
                                                                                                       cAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGt                                               >  1:3870638/1‑150 (MQ=255)
                                                                                                        aGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGtt                                              <  2:68578/150‑1 (MQ=255)
                                                                                                            aCCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGcc                                          <  2:4453011/150‑1 (MQ=255)
                                                                                                             ccTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGccc                                         <  2:3015804/150‑1 (MQ=255)
                                                                                                                 cGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGc                                     >  1:796683/1‑150 (MQ=255)
                                                                                                                   ttCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCaa                                   <  2:2117677/150‑1 (MQ=255)
                                                                                                                      aCGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGcc                                >  1:631766/1‑150 (MQ=255)
                                                                                                                        gAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGc                              >  2:2890121/1‑150 (MQ=255)
                                                                                                                             ttGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTg                         >  1:2804548/1‑150 (MQ=255)
                                                                                                                              tGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGa                        >  1:338726/1‑150 (MQ=255)
                                                                                                                                  ttACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGcc                    >  1:1838751/1‑150 (MQ=255)
                                                                                                                                   gtaaacaccGTGGCGCCGTGCAGGCCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCg                   >  2:1548669/8‑150 (MQ=255)
                                                                                                                                    aCAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGt                  >  1:1761890/1‑150 (MQ=255)
                                                                                                                                        cgcCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCt              <  1:2292746/150‑1 (MQ=255)
                                                                                                                                        caccGTGGCGCCGTGCAGGCCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCt              >  1:1331305/3‑150 (MQ=255)
                                                                                                                                         gcCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCa                                                                                       >  2:229596/1‑76 (MQ=255)
                                                                                                                                         gcCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTg             <  2:4352759/150‑1 (MQ=255)
                                                                                                                                         accGTGGCGCCGTGCAGGCCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTg             >  1:3282118/2‑150 (MQ=255)
                                                                                                                                            gTGGCGGCCTGTACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATc          <  1:4326665/150‑1 (MQ=255)
                                                                                                                                               gcggcCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCACCCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATCttt       >  1:4182730/1‑150 (MQ=255)
                                                                                                                                                    cTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATCTTTAACgg  >  1:1026527/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGTGGTTGCACTGGCAGAGTTGGCTGGGCGAGCGTGCCGGCGCGCTGAAAAAGGCACCGACCCGCCTGGGCAACCGTACCTATGCCGCCGGCGAAGCCTCCCCAGGCACCTACGTTCACGAGCGTTGAGTTACAGCGCCGTGGCGGCCTGCACGGCGCCACGGTGTTTACTTCACCCAAGAGTCACGTGCATGCCGCAACCCTATATTTTCAGGACCGTCGAGCTGGACAACCAGTCCATCCGCACCGCTGTTCGCCCCGGCAAGCCGCACCTGACGCCGTTGCTGATCTTTAACGG  >  NC_002947/5703052‑5703349

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: