Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,972,572 G→A 6.3% intergenic (‑363/+31) fliC ← / ← PP_4379 flagellin/3‑oxoacyl‑ACP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,972,5720GA 6.3% 539.0 / 17.5 175intergenic (‑363/+31)fliC/PP_4379flagellin/3‑oxoacyl‑ACP synthase
Reads supporting (aligned to +/- strand):  ref base G (72/92);  new base A (7/4);  total (79/96)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 186 total aligned reads displayed.

CATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAACACTGGAGG  >  NC_002947/4972423‑4972721
                                                                                                                                                      |                                                                                                                                                     
    ccGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGgccgc                                                                                                                                                   >  1:442150/1‑150 (MQ=255)
    ccGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGgccgc                                                                                                                                                   >  2:1594978/1‑150 (MQ=255)
      gTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGccgccg                                                                                                                                                 <  2:3806442/150‑1 (MQ=255)
         ggTGGCGCGGGTGAGGCTGCGAAGACGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCGGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGgc                                                                                                                                              <  1:1772173/150‑1 (MQ=255)
         ggTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGgc                                                                                                                                              <  1:5432493/150‑1 (MQ=255)
         ggTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCGCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGgc                                                                                                                                              <  2:5018570/150‑1 (MQ=255)
            ggcgggggTGAGGCTGCGAAGATGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGATGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCt                                                                                                                                           <  2:3619986/145‑1 (MQ=255)
               gggggTGCGGCTGCTAAGAAGGAGACGCGCGGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCACCTACAAGGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCACAAAAAAAAGCGCCCGGCGGCCACCGCCCGGCGCtttt                                                                                                                                        <  1:1632770/148‑1 (MQ=37)
                 gggTGAGGCTGCGAAGACGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTGCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAg                                                                                                                                      <  2:1809137/150‑1 (MQ=255)
                 gggTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAg                                                                                                                                      <  2:4438754/150‑1 (MQ=255)
                 gggTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAg                                                                                                                                      <  1:819600/150‑1 (MQ=255)
                        gCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGa                                                                                                                               <  2:3347687/150‑1 (MQ=255)
                           gCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTg                                                                                                                            >  2:4971317/1‑150 (MQ=255)
                             gaagaaGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTggg                                                                                                                          >  1:676279/1‑150 (MQ=255)
                             gaagaaGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTggg                                                                                                                          >  1:2623693/1‑150 (MQ=255)
                             gaagaaGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTggg                                                                                                                          >  1:3322906/1‑150 (MQ=255)
                              atgatggAGCCGCGAGGGGTGGCACGGGCTTTGCGGGTTTTCGCGGGTTAACCCGCTCCTACAAGGGGCTGGGTTCGGCAGAAGCTTATTTCCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGa                                                                                                                         <  1:4660018/145‑1 (MQ=255)
                                gacggAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATc                                                                                                                       <  1:2061718/147‑1 (MQ=255)
                                  cggAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAg                                                                                                                     <  1:3112060/149‑1 (MQ=255)
                                   ggAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGt                                                                                                                    >  1:4662241/1‑150 (MQ=255)
                                   ggAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGt                                                                                                                    >  1:3137206/1‑150 (MQ=255)
                                   ggAGACGCGAGGGTTGGCACCGGCTTTGCCGGTGTTCGCGGGTTAACCCGCTCCTACAAGGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGt                                                                                                                    <  1:4297997/150‑1 (MQ=255)
                                    gagaCGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTc                                                                                                                   >  1:551720/1‑150 (MQ=255)
                                      gaCGCGATGGGTGGCACGGGATTTGCCGGTGTTCGCGGGTCAACCCGCTCCAAAAATGGGCTGCTTTAGCCATAAGCTTATTTTCGGCACAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTcgc                                                                                                                 <  1:21871/150‑1 (MQ=255)
                                      gaCGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTcgc                                                                                                                 >  1:457289/1‑150 (MQ=255)
                                       aCGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGAAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTcgcg                                                                                                                <  1:4706755/150‑1 (MQ=255)
                                        cgcgATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACCATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTcgcg                                                                                                                <  2:2225731/150‑1 (MQ=255)
                                         gcgATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGgc                                                                                                                                              <  2:1178741/118‑1 (MQ=255)
                                         gcgATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGt                                                                                                              >  2:4789167/1‑150 (MQ=255)
                                            aTGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACa                                                                                                           <  2:332427/150‑1 (MQ=255)
                                            aTGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACa                                                                                                           >  1:2361536/1‑150 (MQ=255)
                                             ccgcgaacACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAg                                                                                                          >  1:4760621/8‑150 (MQ=255)
                                             ccgcgaacACCGGCAAAGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAgg                                                                                                         >  2:1746830/8‑150 (MQ=39)
                                              gggTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAgg                                                                                                         <  2:2988468/150‑1 (MQ=255)
                                              gggTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACAGGGATCAGTCGCGGTACAgg                                                                                                         >  2:1166210/1‑150 (MQ=255)
                                                gTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGAt                                                                                                       <  1:2211470/150‑1 (MQ=255)
                                                 tGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATg                                                                                                      <  2:1482115/150‑1 (MQ=255)
                                                     aCCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCg                                                                                                  >  1:4205836/1‑150 (MQ=255)
                                                     aCCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCg                                                                                                  >  1:1459060/1‑150 (MQ=255)
                                                        ggCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGc                                                                                               <  2:1859277/150‑1 (MQ=255)
                                                        ggCTTGGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGc                                                                                               <  1:4633780/150‑1 (MQ=255)
                                                            ttGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCa                                                                                           <  1:4115314/150‑1 (MQ=255)
                                                             tGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCAc                                                                                          <  1:3605715/150‑1 (MQ=255)
                                                                cGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGAc                                                                                       >  1:369065/1‑150 (MQ=255)
                                                                cGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCCGGCGGCCACCGCCCGGCGCtttt                                                                                                                                        >  2:3954601/1‑101 (MQ=38)
                                                                  gtgtTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAg                                                                                     >  1:2949403/1‑150 (MQ=255)
                                                                        gcgGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCAc                                                                               <  1:1188500/150‑1 (MQ=255)
                                                                          gggTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGc                                                                             <  1:1594978/150‑1 (MQ=255)
                                                                            gTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCa                                                                           >  1:885556/1‑150 (MQ=255)
                                                                              gAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCaaa                                                                         <  1:3675223/150‑1 (MQ=255)
                                                                              gAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCaaa                                                                         >  1:2464826/1‑150 (MQ=255)
                                                                              gAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTAGTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCaaa                                                                         <  2:2037133/150‑1 (MQ=255)
                                                                               aaCCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAg                                                                        <  1:1223518/150‑1 (MQ=255)
                                                                               aaCCAGCTCCTACAATGGGCTGAGTTAGTCAGAAGCTTATTTTATGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCTGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCAGAGCCCCACGACAGGCCCACGCAAAAg                                                                        >  2:5199326/1‑150 (MQ=255)
                                                                                 cccGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGTCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGcc                                                                      >  2:3183924/1‑150 (MQ=255)
                                                                                  ccGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCa                                                                     >  1:1326085/1‑150 (MQ=255)
                                                                                   cGCTCCTACACTGGGCTGCGTTAGGCAGAAGATTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCAc                                                                    <  2:3649515/150‑1 (MQ=255)
                                                                                   cGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCAc                                                                    <  1:4153702/150‑1 (MQ=255)
                                                                                        cTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATc                                                               <  1:2455557/150‑1 (MQ=255)
                                                                                               gggCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCAcgc                                                        <  2:1459060/150‑1 (MQ=255)
                                                                                               gggCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCAcgc                                                        <  2:4205836/150‑1 (MQ=255)
                                                                                                 gCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGAt                                                                >  1:2138123/1‑140 (MQ=255)
                                                                                                 gCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACgcgc                                                      >  2:2525458/1‑150 (MQ=255)
                                                                                                 gCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACgcgc                                                      >  2:5252053/1‑150 (MQ=255)
                                                                                                   tGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCtt                                                    >  2:924826/1‑150 (MQ=255)
                                                                                                   tGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCtt                                                    >  2:4665796/1‑150 (MQ=255)
                                                                                                       ttAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAg                                                <  1:622786/150‑1 (MQ=255)
                                                                                                         aGGCAGACGCTTATTTTCGGCAGAAAAAAAAGCGGCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCt                                              <  2:4069355/150‑1 (MQ=255)
                                                                                                         aGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCt                                              <  1:5100817/150‑1 (MQ=255)
                                                                                                             aGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTAc                                                                                                            <  2:4463619/84‑1 (MQ=255)
                                                                                                             aGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCa                                          >  2:606153/1‑150 (MQ=255)
                                                                                                             aGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCa                                          >  1:1536998/1‑150 (MQ=255)
                                                                                                               aaGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCaaa                                                                         >  2:2562108/1‑117 (MQ=255)
                                                                                                                aGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCg                                       <  2:3424068/150‑1 (MQ=255)
                                                                                                                 gCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAg                                                                                     <  2:81767/103‑1 (MQ=255)
                                                                                                                 gCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGa                                      <  1:1082444/150‑1 (MQ=255)
                                                                                                                 gCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGa                                      <  1:1091469/150‑1 (MQ=255)
                                                                                                                      ttttCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACCGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCAtg                                 <  2:4342071/150‑1 (MQ=255)
                                                                                                                       tttCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAgg                                                                                                         <  1:2945929/77‑1 (MQ=255)
                                                                                                                             cAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCt                          <  2:676279/150‑1 (MQ=255)
                                                                                                                              aGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTc                           >  2:3847239/1‑148 (MQ=255)
                                                                                                                              aGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTg                         <  1:3830682/150‑1 (MQ=255)
                                                                                                                               gAAAAAAATGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTgc                        >  2:3812393/1‑150 (MQ=255)
                                                                                                                                  cctaaacGCGCCCGGCGGCCACCGCCCGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCAc                                                                    <  1:4972356/101‑1 (MQ=37)
                                                                                                                                  cctaaaaGCGCCCGGCGGCCACCGCCCGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCAc                                                                    >  2:4972356/4‑104 (MQ=37)
                                                                                                                                   aaaaaaGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTgcagc                     >  1:1976965/1‑150 (MQ=255)
                                                                                                                                      aaaGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTgcagcagc                  <  2:2623693/150‑1 (MQ=255)
                                                                                                                                        aGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGcagcagcag                >  1:3167124/1‑150 (MQ=255)
                                                                                                                                        aGCGCCCGGCGGCCACCGCCCGGCGCTTTTTGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGcagcagcag                >  1:189416/1‑150 (MQ=255)
                                                                                                                                           gcCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCgg             >  1:3934803/1‑150 (MQ=255)
                                                                                                                                           gcCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCgg             >  1:5279093/1‑150 (MQ=255)
                                                                                                                                           gcCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTAGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCgg             >  1:2874547/1‑150 (MQ=255)
                                                                                                                                             cGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGaa           <  2:1046503/150‑1 (MQ=255)
                                                                                                                                             cGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGaa           >  1:3526722/1‑150 (MQ=255)
                                                                                                                                             cGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGaa           >  1:3491612/1‑150 (MQ=255)
                                                                                                                                                 cGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAACCCt       <  2:2911617/150‑1 (MQ=255)
                                                                                                                                                  ggTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAACACTg      <  2:521774/150‑1 (MQ=255)
                                                                                                                                                    tGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAACACTgga    <  2:457793/150‑1 (MQ=255)
                                                                                                                                                    tGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAACACTgga    <  1:4971317/150‑1 (MQ=255)
                                                                                                                                                      gccgccGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAACACTgga    >  1:3961354/1‑148 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAACACTGGAGG  >  NC_002947/4972423‑4972721

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