Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 721,604 G→A R44C (CGT→TGT)  kdpE ← DNA‑binding transcriptional activator KdpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913721,6040GA83.8% 164.6 / 24.2 69R44C (CGT→TGT) kdpEDNA‑binding transcriptional activator KdpE
Reads supporting (aligned to +/- strand):  ref base G (6/5);  new base A (22/35);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CGGGCGCGGTGGTGGCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCGT  >  NC_000913/721358‑721839
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           
cgggcctgttggtggCAGAGTGGCGGAGTAATGTGGCGCGCAGACGGGCCTGCAATTAGCCAATGCCAAACGGCTTACTCAGATTATCATCCGATCCTGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCCCAATCCCCGGCAACGCGCTCCACTGGGGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGACAGGCCGAGATAG‑AGAATAATCAAATCTGGCTTACAGGtt                                                                                                                                                                                                                                         <  1:252521/242‑1 (MQ=255)
 gggCGCGGTGGTGGCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTg                                                                                                                                                                                                                                        >  1:384286/1‑251 (MQ=255)
       ggtggaGGCAGAGTGGAGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCtt                                                                                                                                                                                                                                  >  1:144624/1‑251 (MQ=255)
           gtggCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACGGGTTGCCGCTTcc                                                                                                                                                                                                                                >  1:51350/1‑249 (MQ=255)
            tggCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGc                                                                                                                                                                                                                             >  1:345530/1‑251 (MQ=255)
                agagTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGc                                                                                                                                                                                                                         >  1:27703/1‑251 (MQ=255)
                          cgTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAg                                                                                                                                                                                                               <  1:209696/251‑1 (MQ=255)
                            tAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑CGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTTCAGCg                                                                                                                                                                                                             >  2:330093/1‑251 (MQ=255)
                             aaTGCGACGCCGAGACGGGCCTGCAATTCGCCAATGGCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGAGCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAATCCGCGTTTCAGCGt                                                                                                                                                                                                            <  2:335172/251‑1 (MQ=255)
                              aTGCGACGCCCAGACCGGCCTGCAATTCGCCAATGCCAAACTGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCCCCCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCTGGATCG‑CGAATAACCAAATCTTGCGTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGtt                                                                                                                                                                                                           <  2:191479/251‑1 (MQ=255)
                                  gACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCgg                                                                                                                                                                                                       <  2:39291/251‑1 (MQ=255)
                                     gcgcAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCt                                                                                                                                                                                                    >  2:90049/1‑251 (MQ=255)
                                           aCGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGa                                                                                                                                                                                              <  2:10813/251‑1 (MQ=255)
                                            cGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATGATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGAc                                                                                                                                                                                             <  2:156280/251‑1 (MQ=255)
                                              ggCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGAcgc                                                                                                                                                                                           >  2:417691/1‑251 (MQ=255)
                                                    cAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATcc                                                                                                                                                                                     <  1:124900/251‑1 (MQ=255)
                                                        tCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGt                                                                                                                                                                                  <  2:233022/250‑1 (MQ=255)
                                                          gCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGc                                                                                                                                                                               >  1:256483/1‑251 (MQ=255)
                                                           ccAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCTTCGcc                                                                                                                                                                              <  1:417691/251‑1 (MQ=255)
                                                             aaTGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGCAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑ACAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCTCATCCCGTCGcccc                                                                                                                                                                            >  2:355530/1‑250 (MQ=255)
                                                                  cAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAACA‑AATCACCGGCACCGCGCCCCCCTGGCGCCGGTCGCGGGTAAAATCAAACCCCTCCACAACGGGCCAGCCGAGAACGGAGAATA‑TCAAATTCGGCTTACCGGTTGCCGACTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACCCGCCTCCCGTCGCCCTCCaacc                                                                                                                                                                      >  2:292860/1‑248 (MQ=255)
                                                                   aaaCGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCa                                                                                                                                                                                                                <  1:161053/209‑1 (MQ=255)
                                                                   aaaCGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCa                                                                                                                                                                                                                >  2:161053/1‑209 (MQ=255)
                                                                    aaCGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAgcgc                                                                                                                                                                     >  2:339480/1‑251 (MQ=255)
                                                                        gCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTg                                                                                                                                                                 <  1:339480/251‑1 (MQ=255)
                                                                               cAGATAATCAGCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAg                                                                                                                                                            <  2:68584/249‑1 (MQ=255)
                                                                                    aaTCATCCGCTCCTGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGcg                                                                                                                                                      <  1:135144/250‑1 (MQ=255)
                                                                                    aaTCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGcg                                                                                                                                                      <  2:51350/250‑1 (MQ=255)
                                                                                    aaTCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGcg                                                                                                                                                      <  2:161839/250‑1 (MQ=255)
                                                                                                 ggCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGtt                                                                                                                                        >  1:46400/1‑251 (MQ=255)
                                                                                                    aTCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCAt                                                                                                                                     >  2:390799/1‑251 (MQ=255)
                                                                                                         gcgcGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTCCATCTTca                                                                                                                                <  2:394711/251‑1 (MQ=255)
                                                                                                              gcgATTTTGTCGCTCTCTTCGCTGCGGGCGGAAATCACAATCTCCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTGCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCACCTTCAACAAt                                                                                                                           <  2:342013/251‑1 (MQ=255)
                                                                                                                 aTTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATc                                                                                                                          >  1:306443/1‑249 (MQ=255)
                                                                                                                  ttntGTCGCTCTCTTCGTTGCGTGCGGAAAGCACACTCACCGGCACCGCGCTCCACTGGCGCATGTCTCGGATAAACTCAATCCCATCACCATCGTACTGGCCGAGATCG‑AGAATAATCAAATCTTGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCATCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGCGCGAAGAAAGCGACGAGTAGGCTGTTCAGCTTCAACAATCAg                                                                                                                        <  2:450875/250‑1 (MQ=255)
                                                                                                                     tGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAAc                                                                                                                     <  2:283363/250‑1 (MQ=255)
                                                                                                                      gTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGt                                                                                                                   <  2:78975/251‑1 (MQ=255)
                                                                                                                         gCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGATGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTg                                                                                                                >  1:472457/1‑251 (MQ=255)
                                                                                                                           tctctTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTc                                                                                                              >  1:487820/1‑251 (MQ=255)
                                                                                                                            ctctTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTc                                                                                                              >  1:377157/1‑250 (MQ=255)
                                                                                                                                  gCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATc                                                                                                       <  1:338676/251‑1 (MQ=255)
                                                                                                                                     gcgtgcgGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATcc                                                                                                      <  2:345530/249‑1 (MQ=255)
                                                                                                                                         gcgGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCCCCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGTGTTGCAGCGGTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATg                                                                                                 <  2:323249/250‑1 (MQ=255)
                                                                                                                                           ggAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGaa                                                                                               >  1:302395/1‑250 (MQ=255)
                                                                                                                                                     aTCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGt                                                                                    >  1:162265/1‑251 (MQ=255)
                                                                                                                                                       cACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTc                                                                                  >  1:64667/1‑251 (MQ=255)
                                                                                                                                                          cGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAgg                                                                               <  1:230627/251‑1 (MQ=255)
                                                                                                                                                            gCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAgggg                                                                             >  2:206619/1‑251 (MQ=255)
                                                                                                                                                              aCCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCa                                                                           <  2:384286/251‑1 (MQ=255)
                                                                                                                                                                  cgcTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGttt                                                                       <  2:27703/251‑1 (MQ=255)
                                                                                                                                                                      ccACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAg                                                                          <  2:413724/244‑1 (MQ=255)
                                                                                                                                                                      ccACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAg                                                                          >  1:413724/1‑244 (MQ=255)
                                                                                                                                                                       cACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCGCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAGTCGGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTgg                                                                  <  1:317959/251‑1 (MQ=255)
                                                                                                                                                                           ggCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTgggg                                                                <  1:124942/249‑1 (MQ=255)
                                                                                                                                                                              gcAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAgtg                                                           <  1:355282/251‑1 (MQ=255)
                                                                                                                                                                                         aTAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCCACAATCAGCACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGtt                                                                        >  2:285028/1‑227 (MQ=255)
                                                                                                                                                                                         aTAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGtt                                                                        <  1:285028/227‑1 (MQ=255)
                                                                                                                                                                                           aanCTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAAc                                               <  2:447475/250‑1 (MQ=255)
                                                                                                                                                                                                cAATCCCATCACCAGCGTGCAGGCCGAGATCG‑CGAATAATCAAATCTGAGTTACAGGTGGACGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGAATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGGTCATCTTCAACAATCAGAACGTTTGTCACATAGCCTCATGAAATTTTTCAAGTTCAGGGGCAGTTTTCTGGGGAAGTGTAACGCGAAAACAGGcac                                         <  2:148845/251‑1 (MQ=255)
                                                                                                                                                                                                       aTCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTc                                  <  1:163367/251‑1 (MQ=255)
                                                                                                                                                                                                            cATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCcggtc                             <  2:46400/251‑1 (MQ=255)
                                                                                                                                                                                                                    aGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAAc                     >  1:443097/1‑251 (MQ=255)
                                                                                                                                                                                                                     gtccGAGATCG‑AGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAAcg                    >  1:76791/3‑251 (MQ=255)
                                                                                                                                                                                                                        cGAGATCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCgg                 >  1:93867/1‑251 (MQ=255)
                                                                                                                                                                                                                            aTCG‑AGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCCCCACCTTCCGGTCGGTTGAACGCGGTAAt             <  2:192031/251‑1 (MQ=255)
                                                                                                                                                                                                                               g‑agaATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGt          <  1:206619/251‑1 (MQ=255)
                                                                                                                                                                                                                                 agaATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGtt                         >  2:420609/1‑235 (MQ=255)
                                                                                                                                                                                                                                 agaATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGtt                         <  1:420609/235‑1 (MQ=255)
                                                                                                                                                                                                                                         cAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACTCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCgt  <  1:448653/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           
CGGGCGCGGTGGTGGCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCG‑AGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCGT  >  NC_000913/721358‑721839

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑