Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 625,067 A→T 6.6% intergenic (‑158/+25) 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,0670AT 6.6% 803.1 / 31.3 248intergenic (‑158/+25)PP_0537/ppatranscriptional regulator/inorganic pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (127/103);  new base T (6/10);  total (135/113)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 250 total aligned reads displayed.

TAGCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGCTCGA  >  NC_002947/624921‑625216
                                                                                                                                                  |                                                                                                                                                     
tAGCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGc                                                                                                                                                    >  2:1036921/1‑150 (MQ=255)
   cggcggGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCggg                                                                                                                                                 <  1:2986270/150‑1 (MQ=255)
       gggTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGGCGGGGttt                                                                                                                                             <  1:2890045/150‑1 (MQ=255)
       gggTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGGCGGGGttt                                                                                                                                             <  1:412815/150‑1 (MQ=255)
       gggTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGttt                                                                                                                                             >  2:3864721/1‑150 (MQ=255)
         gTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGttttt                                                                                                                                           <  1:1336389/150‑3 (MQ=255)
          ttGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGGCGGGGtttttc                                                                                                                                          <  2:2591207/150‑4 (MQ=255)
          ttGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTccc                                                                                                                                          >  2:1843229/1‑150 (MQ=255)
           tGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCt                                                                                                                                         <  2:3496916/150‑1 (MQ=255)
               aGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTCGGATGCTGGTGTAGGACTGGGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCACAGCGGGGTTTCCCTGCAg                                                                                                                                     >  2:4559152/1‑150 (MQ=255)
                 tAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTc                                                                                                                                   <  1:486427/150‑1 (MQ=255)
                   aGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGa                                                                                                                                 >  1:967282/1‑150 (MQ=255)
                    ggACATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTTGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGaa                                                                                                                                >  2:4833008/1‑150 (MQ=255)
                     gAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGaaa                                                                                                                               >  1:46071/1‑150 (MQ=255)
                     gAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGaaa                                                                                                                               <  1:5113408/150‑1 (MQ=255)
                      aaaTGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAAt                                                                                                                              <  1:2103057/150‑1 (MQ=255)
                            aCCGGGTGTTGACGGGATTTGTACGAAAAGAGGCGGTTTGTCGGATGCTGGGGTAGGACGGGGTAGGATTTGTTCGCGGCAGTGAGCGCGTGCCATTCAGGCATGAAAAACCCCGCCAGTGCGGGGTTTCCCTGCAGTCGAAAGCCGccc                                                                                                                        >  2:807735/1‑150 (MQ=255)
                                     tGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCa                                                                                                               <  1:5436360/150‑1 (MQ=255)
                                      gACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAgg                                                                                                                 <  2:1247738/147‑1 (MQ=255)
                                         gggATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGa                                                                                                           >  2:1315859/1‑150 (MQ=255)
                                                gTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTTTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGAttt                                                                                                    <  1:3391876/150‑1 (MQ=255)
                                                 tACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTg                                                                                                   <  1:2317228/150‑1 (MQ=255)
                                                 tACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTg                                                                                                   >  2:3714684/1‑150 (MQ=255)
                                                 tACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTg                                                                                                   <  1:4074012/150‑1 (MQ=255)
                                                  aCGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTgg                                                                                                  <  1:4741201/150‑1 (MQ=255)
                                                   cGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGt                                                                                                 >  2:3084325/1‑150 (MQ=255)
                                                   cGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGGAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGt                                                                                                 <  2:5240150/150‑1 (MQ=255)
                                                   cGAAAAGAGGCTGATTGTAGGATGCGGGAGGAGGACTGTGTAGGAATAGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGAtttcgc                                                                                                 <  2:5459126/150‑4 (MQ=255)
                                                    gAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTg                                                                                                >  2:4488573/1‑150 (MQ=255)
                                                    gAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTg                                                                                                >  1:5370777/1‑150 (MQ=255)
                                                         gagGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATcgc                                                                                           <  2:2654329/150‑1 (MQ=255)
                                                          agGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATcgcg                                                                                          <  1:5448859/150‑1 (MQ=255)
                                                           ggCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAGACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCgg                                                                                         <  2:1823928/150‑1 (MQ=255)
                                                              tGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCAc                                                                                      >  2:4596115/1‑150 (MQ=255)
                                                              tGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCAc                                                                                      >  1:254064/1‑150 (MQ=255)
                                                              tGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCAc                                                                                      <  1:5273469/150‑1 (MQ=255)
                                                               gTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACg                                                                                     >  1:1795504/1‑150 (MQ=255)
                                                                  tGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGggc                                                                                  >  2:3846163/1‑150 (MQ=255)
                                                                   gTAGGATGCTGGTGTAGGACTGTGTAGGATTTGGTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGggcg                                                                                 <  2:4151547/150‑1 (MQ=255)
                                                                    tAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTGGTAGGCAGCGACCGATGTGGTGATCGCGGCACGggcgg                                                                                >  2:5337896/1‑150 (MQ=255)
                                                                     aGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCTCGTTCCTTTAAGGAATGAAAATCCCCGCCAAAGCGGTTTTTCCCTGCAGTCGAATTAAGCCCTTGTAGGCTGCGTCCGATTTGGTGATCGCGGCACGggcggc                                                                               <  2:2703513/150‑1 (MQ=255)
                                                                      ggATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTg                                                                                                   >  2:1889313/1‑129 (MQ=255)
                                                                      ggATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGgcggcg                                                                              >  2:2244187/1‑150 (MQ=255)
                                                                        aTGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTc                                                                            >  1:5245779/1‑150 (MQ=255)
                                                                            tGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGgcg                                                                        >  1:493170/1‑150 (MQ=255)
                                                                              gtgtAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGcc                                                                      >  1:1710592/1‑150 (MQ=255)
                                                                              gtgtAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGcc                                                                      >  1:5198804/1‑150 (MQ=255)
                                                                              gtgtAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGcc                                                                      >  1:4545932/1‑150 (MQ=255)
                                                                                 tAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTc                                                                   <  2:960704/150‑1 (MQ=255)
                                                                                  aGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGCCGACgccgcccgtg                                                                  <  1:1468025/150‑5 (MQ=255)
                                                                                         tgtAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATc                                                                                             >  1:2428992/1‑116 (MQ=255)
                                                                                         tgtAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCtt                                                           >  1:115675/1‑150 (MQ=255)
                                                                                          gtAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTc                                                          >  2:935499/1‑150 (MQ=255)
                                                                                           tAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCt                                                                                                                                         <  2:3990249/70‑1 (MQ=255)
                                                                                           tAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCg                                                         >  1:871897/1‑150 (MQ=255)
                                                                                           tAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCg                                                         >  2:2111395/1‑150 (MQ=255)
                                                                                           tAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCAAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCCGCCTTCg                                                         >  2:3669420/1‑150 (MQ=255)
                                                                                                   gTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGAc                                                 >  1:1963730/1‑150 (MQ=255)
                                                                                                    ttCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGAcc                                                >  1:1413649/1‑150 (MQ=255)
                                                                                                     tCCCAGCAGTGAGGGCGTTCCATTCTGGCATGAAAATCCCCGCCAATGCGGGGTTTCGCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGAGCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGTTCTTGAccc                                               >  2:4888037/1‑150 (MQ=255)
                                                                                                      cccAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCa                                              <  1:1005653/150‑1 (MQ=255)
                                                                                                       ccAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCAc                                             >  2:3144761/1‑150 (MQ=255)
                                                                                                        cagcagTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACt                                            >  1:2235315/1‑150 (MQ=255)
                                                                                                          gcagTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTg                                          <  2:4730684/150‑1 (MQ=255)
                                                                                                               gAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCtt                                     >  1:535100/1‑150 (MQ=255)
                                                                                                               gAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCtt                                     >  1:3056512/1‑150 (MQ=255)
                                                                                                                aGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTtc                                    <  1:4405351/150‑1 (MQ=255)
                                                                                                                aGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTtc                                    >  2:1119869/1‑150 (MQ=255)
                                                                                                                 gcgcGTTCCATTCAGGCCTGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCGTCCCAGCCTTCGATCTTGACCCACTTGCCCTtct                                   <  2:2895501/150‑1 (MQ=255)
                                                                                                                 gcgcGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTtct                                   >  2:1678076/1‑150 (MQ=255)
                                                                                                                  cgcgTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTtctc                                  >  1:649927/1‑150 (MQ=255)
                                                                                                                    cgTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGTCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCTCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCga                                <  1:3535209/150‑1 (MQ=255)
                                                                                                                      ttCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCgaga                              <  2:606956/150‑1 (MQ=255)
                                                                                                                      ttCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCgaga                              <  2:4087214/150‑1 (MQ=255)
                                                                                                                      ttCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCgaga                              <  1:3800751/150‑1 (MQ=255)
                                                                                                                           ttCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCTCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                         <  1:348640/150‑1 (MQ=255)
                                                                                                                           ttCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                         <  2:2943879/150‑1 (MQ=255)
                                                                                                                           ttCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                         <  2:4774300/150‑1 (MQ=255)
                                                                                                                           ttCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                         <  1:1843202/150‑1 (MQ=255)
                                                                                                                             cAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGt                       <  1:1245533/150‑1 (MQ=255)
                                                                                                                              aGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTa                      >  2:2335269/1‑150 (MQ=255)
                                                                                                                              aGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTa                      <  2:563909/150‑1 (MQ=255)
                                                                                                                               ggCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAg                     <  2:345971/150‑1 (MQ=255)
                                                                                                                                 cATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGtt                   <  2:1537228/150‑1 (MQ=255)
                                                                                                                                  aTGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTg                  >  1:360540/1‑150 (MQ=255)
                                                                                                                                  aTGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTg                  >  1:2147504/1‑150 (MQ=255)
                                                                                                                                  aTGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTg                  >  1:4795597/1‑150 (MQ=255)
                                                                                                                                  aTGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTg                  <  2:297182/150‑1 (MQ=255)
                                                                                                                                   tGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTgg                 <  1:4588622/150‑1 (MQ=255)
                                                                                                                                   tGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTgg                 >  1:1060593/1‑150 (MQ=255)
                                                                                                                                    gAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGc                >  2:4357115/1‑150 (MQ=255)
                                                                                                                                     ggaaaCCCCGCTTTGGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCa               >  2:133157/3‑150 (MQ=255)
                                                                                                                                     ggaaaCCCCGCTTTGGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCa               >  2:586025/3‑150 (MQ=255)
                                                                                                                                     aaaaaCCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCa               >  1:5550452/1‑150 (MQ=255)
                                                                                                                                          ccccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaagaa          <  2:5165085/150‑1 (MQ=255)
                                                                                                                                           cccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaagaag         >  2:4438717/1‑150 (MQ=255)
                                                                                                                                              gCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGc      >  1:3122829/1‑150 (MQ=255)
                                                                                                                                                cAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGCTc    >  2:1160315/1‑150 (MQ=255)
                                                                                                                                                cAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGCTc    <  1:4178732/150‑1 (MQ=255)
                                                                                                                                                  aaGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGCTCGa  >  1:370573/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TAGCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGCTCGA  >  NC_002947/624921‑625216

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