Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 625,068 A→G 7.3% intergenic (‑159/+24) 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,0680AG 7.3% 669.3 / 27.7 220intergenic (‑159/+24)PP_0537/ppatranscriptional regulator/inorganic pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (112/90);  new base G (7/9);  total (119/101)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

Only 100 of 224 total aligned reads displayed.

AGCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTG  >  NC_002947/624922‑625211
                                                                                                                                                  |                                                                                                                                               
 gcggcgGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCTTTCAGGCATGAAAAACCCCGCCAAAGCgg                                                                                                                                             >  2:4741456/1‑150 (MQ=255)
  cggcggGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCggg                                                                                                                                            >  1:803961/1‑150 (MQ=255)
  cggcggGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCggg                                                                                                                                            <  2:1613218/150‑1 (MQ=255)
  cggcggGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCggg                                                                                                                                            <  1:3515847/150‑1 (MQ=255)
  cggcggGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCggg                                                                                                                                            <  1:5108472/150‑1 (MQ=255)
      gggTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGGCGGGGttt                                                                                                                                        <  2:446404/150‑1 (MQ=255)
       ggTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGGCGGGGtttt                                                                                                                                       <  1:85196/150‑2 (MQ=255)
        gTTGCCAGTAAGGAAATGACCCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTcc                                                                                                                                      <  2:4535299/150‑1 (MQ=255)
        gTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTcc                                                                                                                                      <  1:3981668/150‑1 (MQ=255)
         ttGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGGCGGGGtttttc                                                                                                                                     <  1:3869552/150‑4 (MQ=255)
         ttGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGGCGGGGtttttc                                                                                                                                     <  1:3869478/150‑4 (MQ=255)
              aGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAg                                                                                                                                >  1:2458476/1‑150 (MQ=255)
                  aGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGa                                                                                                                            >  1:899592/1‑150 (MQ=255)
                    gAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGaaa                                                                                                                          >  1:3556485/1‑150 (MQ=255)
                      aaTGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCggg                                                                                                                                            >  2:979866/1‑130 (MQ=255)
                      aaTGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATc                                                                                                                        >  1:5223914/1‑150 (MQ=255)
                      aaTGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATc                                                                                                                        >  1:2996504/1‑150 (MQ=255)
                      aaTGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATc                                                                                                                        >  2:1306306/1‑150 (MQ=255)
                      aaTGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATc                                                                                                                        >  1:2996491/1‑150 (MQ=255)
                       aTGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCa                                                                                                                       >  1:4875767/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)
                               ggTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGt                                                                                                               <  1:1228073/150‑1 (MQ=255)
                                    tGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCa                                                                                                          <  2:2454896/150‑1 (MQ=255)
                                      aCGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCATCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCgg                                                                                                                                             <  2:1430180/113‑1 (MQ=255)
                                               gTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGAttt                                                                                               >  1:1193926/1‑150 (MQ=255)
                                                tACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTa                                                                                                              >  1:3383731/1‑134 (MQ=255)
                                                tACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTa                                                                                                              <  2:3383731/134‑1 (MQ=255)
                                                  cGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGt                                                                                            >  1:1804674/1‑150 (MQ=255)
                                                   gAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCGGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTg                                                                                           >  2:1809048/1‑150 (MQ=255)
                                                   gAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTg                                                                                                                <  1:1558478/129‑1 (MQ=255)
                                                            cTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCa                                                                                  >  2:294692/1‑150 (MQ=255)
                                                             tGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCTCGCCAAAGCGGGGTTTCCCTGCAGTCTAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCAc                                                                                 <  1:4746330/150‑1 (MQ=255)
                                                             tGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCAc                                                                                 >  1:3137176/1‑150 (MQ=255)
                                                              gTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACg                                                                                >  2:493841/1‑150 (MQ=255)
                                                               tttGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGAcc                                                                                                    >  2:1088930/1‑129 (MQ=255)
                                                               tttGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGAcc                                                                                                    <  1:1088930/129‑1 (MQ=255)
                                                               tttGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACgg                                                                               >  2:4878811/1‑150 (MQ=255)
                                                                   tAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGggcgg                                                                           >  2:312014/1‑150 (MQ=255)
                                                                      gATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGt                                                                        >  1:4317574/1‑150 (MQ=255)
                                                                      gATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGt                                                                        >  1:2612480/1‑150 (MQ=255)
                                                                           tGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGgcg                                                                   >  2:2701456/1‑150 (MQ=255)
                                                                            ggTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGgcgc                                                                  <  2:1238768/150‑1 (MQ=255)
                                                                             gtgtAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGcc                                                                 >  2:4038869/1‑150 (MQ=255)
                                                                             gtgtAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGcc                                                                 >  2:3164276/1‑150 (MQ=255)
                                                                                tAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTc                                                              >  2:2122595/1‑150 (MQ=255)
                                                                                tAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTc                                                              <  1:4368581/150‑1 (MQ=255)
                                                                                   gACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCa                                                           >  2:2842304/1‑150 (MQ=255)
                                                                                   gACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCa                                                           >  2:2414694/1‑150 (MQ=255)
                                                                                      tgtgTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATATGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGcc                                                        >  2:165175/1‑150 (MQ=255)
                                                                                        tgtAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCtt                                                      <  2:2612480/150‑1 (MQ=255)
                                                                                        tgtAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCtt                                                      >  2:2318594/1‑150 (MQ=255)
                                                                                          tAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCg                                                    >  1:3086986/1‑150 (MQ=255)
                                                                                           aGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGa                                                   <  1:2132436/150‑1 (MQ=255)
                                                                                            ggATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCa                                                           <  2:2278645/141‑1 (MQ=255)
                                                                                            ggATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCa                                                           >  1:2278645/1‑141 (MQ=255)
                                                                                            ggATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGAt                                                  >  2:4951690/1‑150 (MQ=255)
                                                                                             gATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATc                                                 <  1:2848055/150‑1 (MQ=255)
                                                                                             gATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATc                                                 <  2:4071182/150‑1 (MQ=255)
                                                                                              aTTTGTTCCCAGCAGTGAGCGCGCTCCATTCAGGCCTGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCt                                                <  1:517124/150‑1 (MQ=255)
                                                                                                 tGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGa                                             <  1:367568/150‑1 (MQ=255)
                                                                                                  gTTCCTAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGAc                                            <  1:4847482/150‑1 (MQ=255)
                                                                                                   ttCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGAcc                                           <  1:1016001/150‑1 (MQ=255)
                                                                                                      ccAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCAc                                        <  2:1193926/150‑1 (MQ=255)
                                                                                                       ctgcagTGAGCGCGTTCCATTCAGGCATGTAATACCCCGCCAATGCGGGGTTTCGCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCTGCCTTCGATCTTGACCctct                                       <  2:5448090/148‑1 (MQ=255)
                                                                                                       cagcagTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACt                                       >  1:4352029/1‑150 (MQ=255)
                                                                                                        agcagTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGGGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACtt                                      <  2:5287818/150‑1 (MQ=255)
                                                                                                          cagTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGc                                    <  2:4325194/150‑1 (MQ=255)
                                                                                                            gTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGccc                                  <  2:3402683/150‑1 (MQ=255)
                                                                                                             tGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCt                                 >  1:4626943/1‑150 (MQ=255)
                                                                                                                 cgcgTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTtctc                             <  2:4626943/150‑1 (MQ=255)
                                                                                                                  gcgGTCCAGTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCg                            >  2:5225100/1‑150 (MQ=255)
                                                                                                                        cATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCt                      <  2:1007516/150‑1 (MQ=255)
                                                                                                                          ttCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                    <  1:271232/150‑1 (MQ=255)
                                                                                                                          ttCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                    <  1:4864651/150‑1 (MQ=255)
                                                                                                                           tCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTg                   <  2:3736029/150‑1 (MQ=255)
                                                                                                                               gCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGgcgc                                                                  <  2:3862296/99‑1 (MQ=255)
                                                                                                                               gCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGt               <  1:1856349/150‑1 (MQ=255)
                                                                                                                                cATGAAACACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGtt              <  2:5233704/150‑1 (MQ=255)
                                                                                                                                cATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGtt              >  1:4736605/1‑150 (MQ=255)
                                                                                                                                 aTGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTg             >  2:3668769/1‑150 (MQ=255)
                                                                                                                                 aTGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTg             >  2:926230/1‑150 (MQ=255)
                                                                                                                                  tGAAAAACCCCGGATAAGCGGGGTTTCCCTTCTGTCGAAATCATCCCTTGTATGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTgg            <  2:1855961/150‑1 (MQ=255)
                                                                                                                                  tGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGgc                                                                    >  1:1724077/1‑94 (MQ=255)
                                                                                                                                  tGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTgg            <  1:2571180/150‑1 (MQ=255)
                                                                                                                                   gAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTa                 >  1:3786175/1‑144 (MQ=255)
                                                                                                                                   gAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTa                 <  2:3786175/144‑1 (MQ=255)
                                                                                                                                   gAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGc           <  2:296828/150‑1 (MQ=255)
                                                                                                                                   gAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACGTGCCCGTCTCGAGATCTTTGTAGTTGGc           >  1:3584099/1‑150 (MQ=255)
                                                                                                                                    ggaaaCCCCGCTTTGGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCa          >  1:597081/3‑150 (MQ=255)
                                                                                                                                    ggaaaCCCCGCTTTGGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCa          >  2:866321/3‑150 (MQ=255)
                                                                                                                                         ccccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaagaa     <  1:1416916/150‑1 (MQ=255)
                                                                                                                                         ccccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaagaa     <  2:3036529/150‑1 (MQ=255)
                                                                                                                                         ccccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaagaa     <  2:803961/150‑1 (MQ=255)
                                                                                                                                          cccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaagaag    >  2:1984195/1‑150 (MQ=255)
                                                                                                                                          cccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaagaag    <  1:3118530/150‑1 (MQ=255)
                                                                                                                                           ccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGt   <  1:2204488/150‑1 (MQ=255)
                                                                                                                                           ccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGt   >  1:3334562/1‑150 (MQ=255)
                                                                                                                                                  aGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCa                                                           >  2:1932666/1‑87 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
AGCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTG  >  NC_002947/624922‑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: