Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 625,066 A→T 6.4% intergenic (‑157/+26) 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,0660AT 6.4% 830.9 / 31.1 248intergenic (‑157/+26)PP_0537/ppatranscriptional regulator/inorganic pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (129/103);  new base T (6/10);  total (135/113)
Fisher's exact test for biased strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 249 total aligned reads displayed.

TAGCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGCTCG  >  NC_002947/624921‑625215
                                                                                                                                                 |                                                                                                                                                     
tAGCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGc                                                                                                                                                   >  2:1036921/1‑150 (MQ=255)
  gcggcgGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCgg                                                                                                                                                 >  1:5224710/1‑150 (MQ=255)
       gggTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGGCGGGGttt                                                                                                                                            <  1:4569286/150‑1 (MQ=255)
       gggTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGGCGGGGttt                                                                                                                                            <  1:412815/150‑1 (MQ=255)
       gggTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGGCGGGGttt                                                                                                                                            <  1:2890045/150‑1 (MQ=255)
       gggTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGttt                                                                                                                                            >  2:3864721/1‑150 (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:967639/1‑150 (MQ=255)
                    ggACATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTTGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGaa                                                                                                                               >  2:4833008/1‑150 (MQ=255)
                     gAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGaaa                                                                                                                              >  1:371901/1‑150 (MQ=255)
                     gAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGaaa                                                                                                                              <  1:5113408/150‑1 (MQ=255)
                      aaaTGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAAt                                                                                                                             <  1:2103057/150‑1 (MQ=255)
                           aaCCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGcc                                                                                                                        >  1:2959081/1‑150 (MQ=255)
                            aCCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGccc                                                                                                                       >  1:5298142/1‑150 (MQ=255)
                            aCCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGccc                                                                                                                       >  1:5424/1‑150 (MQ=255)
                            aCCGGGTGTTGACGGGATTTGTACGAAAAGAGGCGGTTTGTCGGATGCTGGGGTAGGACGGGGTAGGATTTGTTCGCGGCAGTGAGCGCGTGCCATTCAGGCATGAAAAACCCCGCCAGTGCGGGGTTTCCCTGCAGTCGAAAGCCGccc                                                                                                                       >  2:807735/1‑150 (MQ=255)
                                 gtgtTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTa                                                                                                                  <  1:1254469/150‑1 (MQ=255)
                                     tGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCa                                                                                                              <  1:306981/150‑1 (MQ=255)
                                     tGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCa                                                                                                              <  1:5436360/150‑1 (MQ=255)
                                      gACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAgg                                                                                                                <  2:1247738/147‑1 (MQ=255)
                                      gACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAg                                                                                                             >  1:1769999/1‑150 (MQ=255)
                                         gggATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGa                                                                                                          >  2:1315859/1‑150 (MQ=255)
                                             tttGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGa                                                                                                      >  1:4625685/1‑150 (MQ=255)
                                                gTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTTTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGAttt                                                                                                   <  1:3391876/150‑1 (MQ=255)
                                                 tACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTg                                                                                                  >  2:3714684/1‑150 (MQ=255)
                                                 tACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTg                                                                                                  <  1:2316548/150‑1 (MQ=255)
                                                 tACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTg                                                                                                  <  1:4074018/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)
                                                     aaaaGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGa                                                                                              <  1:1432212/150‑1 (MQ=255)
                                                         gagGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATcgc                                                                                          <  2:2654329/150‑1 (MQ=255)
                                                           ggCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAGACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCgg                                                                                        <  2:1823928/150‑1 (MQ=255)
                                                              tGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCAc                                                                                     >  2:4596115/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)
                                                                      ggATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTg                                                                                                  <  1:1889313/129‑1 (MQ=255)
                                                                      ggATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGgcggcg                                                                             >  2:2244187/1‑150 (MQ=255)
                                                                        aTGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTc                                                                           >  1:5245779/1‑150 (MQ=255)
                                                                              gtgtAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGcc                                                                     >  1:1710592/1‑150 (MQ=255)
                                                                                 tAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTc                                                                  <  2:960704/150‑1 (MQ=255)
                                                                                        atgtAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCt                                                           <  1:3538070/149‑1 (MQ=255)
                                                                                         tgtAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATc                                                                                            >  1:2428992/1‑116 (MQ=255)
                                                                                         tgtAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCtt                                                          >  1:3441374/1‑150 (MQ=255)
                                                                                         tgtAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCtt                                                          >  1:2501698/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)
                                                                                            aGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGa                                                       <  1:4833008/150‑1 (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)
                                                                                                            agTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGcc                                       <  1:734090/150‑1 (MQ=255)
                                                                                                               gAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCtt                                    >  1:535100/1‑150 (MQ=255)
                                                                                                               gAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCtt                                    >  1:3056512/1‑150 (MQ=255)
                                                                                                                aGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTtc                                   >  1:113654/1‑150 (MQ=255)
                                                                                                                aGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTtc                                   >  2:1119869/1‑150 (MQ=255)
                                                                                                                aGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTtc                                   <  1:4405351/150‑1 (MQ=255)
                                                                                                                 gcgcGTTCCATTCAGGCCTGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCGTCCCAGCCTTCGATCTTGACCCACTTGCCCTtct                                  <  2:2895501/150‑1 (MQ=255)
                                                                                                                 gcgcGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTtct                                  >  2:1678076/1‑150 (MQ=255)
                                                                                                                    cgTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGTCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCTCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCga                               <  1:3535209/150‑1 (MQ=255)
                                                                                                                    cgTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCga                               <  1:4971623/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)
                                                                                                                           ttCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                        <  1:1843202/150‑1 (MQ=255)
                                                                                                                           ttCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                        <  2:2943879/150‑1 (MQ=255)
                                                                                                                           ttCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                        <  2:4774300/150‑1 (MQ=255)
                                                                                                                            tCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTg                       <  1:4262313/150‑1 (MQ=255)
                                                                                                                             cAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGt                      <  1:1245533/150‑1 (MQ=255)
                                                                                                                              aGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTa                     <  2:563909/150‑1 (MQ=255)
                                                                                                                              aGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTa                     >  2:2335269/1‑150 (MQ=255)
                                                                                                                               ggCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAg                    <  2:345971/150‑1 (MQ=255)
                                                                                                                                 cATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGtt                  <  2:1537228/150‑1 (MQ=255)
                                                                                                                                  aTGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTg                 >  1:2147498/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:5538694/1‑150 (MQ=255)
                                                                                                                                     aaaaaCCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCa              >  1:5550452/1‑150 (MQ=255)
                                                                                                                                          ccccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaagaa         <  2:5165085/150‑1 (MQ=255)
                                                                                                                                           cccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaagaag        >  2:4438717/1‑150 (MQ=255)
                                                                                                                                                cAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGCTc   >  2:1160315/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TAGCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCCCTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGCTCG  >  NC_002947/624921‑625215

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: