Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 625,065 C→T 7.2% intergenic (‑156/+27) PP_0537 ← / ← ppa transcriptional regulator/inorganic pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947625,0650CT 7.2% 667.6 / 27.6 221intergenic (‑156/+27)PP_0537/ppatranscriptional regulator/inorganic pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (113/91);  new base T (7/9);  total (120/101)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 222 total aligned reads displayed.

GAGTGTAGCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTG  >  NC_002947/624916‑625211
                                                                                                                                                     |                                                                                                                                                  
gAGTGTAGCGGCGGGTTGCCAGTAAGGAAATGAAACGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGcc                                                                                                                                                    >  2:1878859/1‑150 (MQ=255)
  gtgtAGCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCaa                                                                                                                                                  >  2:1131344/1‑150 (MQ=255)
       gcggcgGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCTTTCAGGCATGAAAAACCCCGCCAAAGCgg                                                                                                                                             >  2:4741456/1‑150 (MQ=255)
        cggcggGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCggg                                                                                                                                            <  2:1613218/150‑1 (MQ=255)
            gggTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGGCGGGGttt                                                                                                                                        <  1:3877957/150‑1 (MQ=255)
            gggTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGGCGGGGttt                                                                                                                                        <  2:446404/150‑1 (MQ=255)
            gggTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGttt                                                                                                                                        >  1:2386257/1‑150 (MQ=255)
              gTTGCCAGTAAGGAAATGACCCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTcc                                                                                                                                      <  2:4535299/150‑1 (MQ=255)
              gTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAGAAACCCCGCTTTGGCGGGGttttt                                                                                                                                      <  1:3226058/150‑3 (MQ=255)
               ttGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGGCGGGGtttttc                                                                                                                                     <  1:3869489/150‑4 (MQ=255)
                   cAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCa                                                                                                                                 >  2:2495642/1‑150 (MQ=255)
                      tAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTc                                                                                                                              >  2:4666159/1‑150 (MQ=255)
                         ggAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGaa                                                                                                                           <  1:880028/150‑1 (MQ=255)
                         ggAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGaa                                                                                                                           >  1:4851545/1‑150 (MQ=255)
                          gAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGaaa                                                                                                                          >  1:3643451/1‑150 (MQ=255)
                          gAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGaaa                                                                                                                          >  1:3556349/1‑150 (MQ=255)
                            aaTGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCggg                                                                                                                                            >  2:979866/1‑130 (MQ=255)
                            aaTGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATc                                                                                                                        >  1:2996498/1‑150 (MQ=255)
                            aaTGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATc                                                                                                                        >  2:1306306/1‑150 (MQ=255)
                               gAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGc                                                                                                                     <  1:2046362/150‑1 (MQ=255)
                                aaCCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGcc                                                                                                                    <  1:404194/150‑1 (MQ=255)
                                  ccGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCt                                                                                                                  >  2:2613816/1‑150 (MQ=255)
                                    gggTGTTGACGGGATTTGTACGAAAAGAGGCTGTTGGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTg                                                                                                                <  2:3713913/150‑1 (MQ=255)
                                     ggTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGt                                                                                                               <  1:1228073/150‑1 (MQ=255)
                                     ggTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGt                                                                                                               >  2:2848055/1‑150 (MQ=255)
                                            aCGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCATCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCgg                                                                                                                                             <  2:1430180/113‑1 (MQ=255)
                                            aCGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGc                                                                                                        <  1:2981818/150‑1 (MQ=255)
                                                gATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGAcc                                                                                                    >  2:1833204/1‑150 (MQ=255)
                                                     gTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGAttt                                                                                               >  1:1193926/1‑150 (MQ=255)
                                                      tACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTa                                                                                                              <  2:3383731/134‑1 (MQ=255)
                                                      tACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTa                                                                                                              >  1:3383731/1‑134 (MQ=255)
                                                       aCGAAAAGATGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCTGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGTCAGCGACCGAttttg                                                                                             <  2:3740360/150‑3 (MQ=255)
                                                        cGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGt                                                                                            >  1:1804674/1‑150 (MQ=255)
                                                         gAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCGGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTg                                                                                           >  2:1809048/1‑150 (MQ=255)
                                                         gAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTg                                                                                                                <  1:1558478/129‑1 (MQ=255)
                                                               agGCTGTTTGTAGGATGCTGGTGTAGGACGGTGTAGGATTTGGTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAACGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATcgcg                                                                                     <  2:3155723/150‑1 (MQ=255)
                                                                  cTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCa                                                                                  >  2:294692/1‑150 (MQ=255)
                                                                    gTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACg                                                                                >  2:493841/1‑150 (MQ=255)
                                                                    gTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACg                                                                                >  1:4324690/1‑150 (MQ=255)
                                                                     tttGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGAcc                                                                                                    <  1:1088930/129‑1 (MQ=255)
                                                                     tttGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACgg                                                                               >  2:4878811/1‑150 (MQ=255)
                                                                       tGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGggc                                                                             <  1:4747038/150‑1 (MQ=255)
                                                                          aGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGggcggc                                                                          >  2:2204488/1‑150 (MQ=255)
                                                                            gATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGt                                                                        >  1:2612480/1‑150 (MQ=255)
                                                                                 tGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGgcg                                                                   >  1:5034041/1‑150 (MQ=255)
                                                                                  ggTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGgcgc                                                                  <  2:1238768/150‑1 (MQ=255)
                                                                                   gtgtAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGcc                                                                 >  2:3164276/1‑150 (MQ=255)
                                                                                      tAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTc                                                              >  2:2122595/1‑150 (MQ=255)
                                                                                       aGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTcc                                                             <  1:312014/150‑1 (MQ=255)
                                                                                         gACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCa                                                           >  2:2414694/1‑150 (MQ=255)
                                                                                            tgtgTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATATGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGcc                                                        >  2:165175/1‑150 (MQ=255)
                                                                                              tgtAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCtt                                                      >  2:2318594/1‑150 (MQ=255)
                                                                                              tgtAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCtt                                                      <  2:2612480/150‑1 (MQ=255)
                                                                                                tAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCg                                                    >  1:3086986/1‑150 (MQ=255)
                                                                                                 aGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGa                                                   <  1:4868110/150‑1 (MQ=255)
                                                                                                 aGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGa                                                   <  1:844033/150‑1 (MQ=255)
                                                                                                 aGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGa                                                   >  1:5140927/1‑150 (MQ=255)
                                                                                                   gATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATc                                                 >  1:436045/1‑150 (MQ=255)
                                                                                                     tttGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCtt                                               <  1:2780985/150‑1 (MQ=255)
                                                                                                         ttCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGAcc                                           >  2:5111481/1‑150 (MQ=255)
                                                                                                         ttCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGAcc                                           <  1:1016001/150‑1 (MQ=255)
                                                                                                         ttCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGAcc                                           <  1:4519916/150‑1 (MQ=255)
                                                                                                            ccAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCAc                                        <  2:1193926/150‑1 (MQ=255)
                                                                                                            ccAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGCGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCAc                                        >  2:3099551/1‑150 (MQ=255)
                                                                                                             ctgcagTGAGCGCGTTCCATTCAGGCATGTAATACCCCGCCAATGCGGGGTTTCGCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCTGCCTTCGATCTTGACCctct                                       <  2:5448090/148‑1 (MQ=255)
                                                                                                              agcagTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGGGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACtt                                      <  2:5287818/150‑1 (MQ=255)
                                                                                                               gcagTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTg                                     <  1:530571/150‑1 (MQ=255)
                                                                                                                cagTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGc                                    <  2:4325194/150‑1 (MQ=255)
                                                                                                                cagTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGc                                    <  2:4059652/150‑1 (MQ=255)
                                                                                                                  gTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGccc                                  <  2:3402683/150‑1 (MQ=255)
                                                                                                                       cgcgTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTtctc                             <  1:3438539/150‑1 (MQ=255)
                                                                                                                        gcgTTCCTTTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCg                            <  1:4741456/150‑1 (MQ=255)
                                                                                                                        gcgGTCCAGTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCg                            >  2:5225100/1‑150 (MQ=255)
                                                                                                                           ttCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCgaga                         >  2:4120829/1‑150 (MQ=255)
                                                                                                                              cATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCt                      >  1:2167466/1‑150 (MQ=255)
                                                                                                                              cATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCt                      <  2:1007604/150‑1 (MQ=255)
                                                                                                                               aTTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCtt                     >  1:465200/1‑150 (MQ=255)
                                                                                                                                ttCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                    >  1:782032/1‑150 (MQ=255)
                                                                                                                                ttCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                    <  1:271232/150‑1 (MQ=255)
                                                                                                                                ttCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                    <  1:904795/150‑1 (MQ=255)
                                                                                                                                   aGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTa                 >  1:2569276/1‑150 (MQ=255)
                                                                                                                                     gCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGgcgc                                                                  <  2:3862296/99‑1 (MQ=255)
                                                                                                                                     gCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGgcgc                                                                  >  1:3862296/1‑99 (MQ=255)
                                                                                                                                     gCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGt               >  1:4307245/1‑150 (MQ=255)
                                                                                                                                     gCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGt               <  1:1856349/150‑1 (MQ=255)
                                                                                                                                      cATGAAACACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGtt              <  2:5233704/150‑1 (MQ=255)
                                                                                                                                       aTGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTg             >  2:926230/1‑150 (MQ=255)
                                                                                                                                        tGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGgc                                                                    >  1:1724077/1‑94 (MQ=255)
                                                                                                                                         gAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTa                 <  2:3786175/144‑1 (MQ=255)
                                                                                                                                         gAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTa                 >  1:3786175/1‑144 (MQ=255)
                                                                                                                                         gAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGc           >  1:3584095/1‑150 (MQ=255)
                                                                                                                                         gAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGc           <  2:296828/150‑1 (MQ=255)
                                                                                                                                          ggaaaCCCCGCTTTGGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCtt                                               >  1:5219463/3‑113 (MQ=255)
                                                                                                                                          ggaaaCCCCGCTTTGGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCa          >  2:866321/3‑150 (MQ=255)
                                                                                                                                          aaaaaCCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCa          >  1:2211241/1‑150 (MQ=255)
                                                                                                                                               ccccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaagaa     <  1:1416916/150‑1 (MQ=255)
                                                                                                                                               ccccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaagaa     <  2:2725891/150‑1 (MQ=255)
                                                                                                                                               ccccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaagaa     <  2:803961/150‑1 (MQ=255)
                                                                                                                                                cccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaagaag    <  1:3118347/150‑1 (MQ=255)
                                                                                                                                                cccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaagaag    >  2:1984195/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GAGTGTAGCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTG  >  NC_002947/624916‑625211

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: