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,6040GA98.2% 190.0 / NA 54R44C (CGT→TGT) kdpEDNA‑binding transcriptional activator KdpE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/34);  minor base C (1/0);  total (20/34)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGCGCGGTGGTGGCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCGTGTACATCCA  >  NC_000913/721358‑721848
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
cGGGCGCGGTGGTGGCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGtt                                                                                                                                                                                                                                                  >  1:126049/1‑251 (MQ=255)
 gggCGCGGTGGTGGCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTg                                                                                                                                                                                                                                                 <  2:493035/251‑1 (MQ=255)
          ggtggCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCa                                                                                                                                                                                                                                        <  1:506337/251‑1 (MQ=255)
           gtggCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAACCGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTTCCGCTTCCAg                                                                                                                                                                                                                                       <  2:239382/251‑1 (MQ=255)
               cAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAg                                                                                                                                                                                                                                   <  2:431707/251‑1 (MQ=255)
                       cggcgTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGtt                                                                                                                                                                                                                                                  >  1:255302/1‑228 (MQ=255)
                       cggcgTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGtt                                                                                                                                                                                                                                                  <  2:255302/228‑1 (MQ=255)
                                 cGACGCGCAGAAGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCg                                                                                                                                                                                                                 >  2:277497/1‑251 (MQ=255)
                                  gACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACCGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCgg                                                                                                                                                                                                                >  2:507749/1‑251 (MQ=255)
                                          gccGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCTCCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAg                                                                                                                                                                                                        <  1:224827/249‑1 (MQ=255)
                                            cGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGAc                                                                                                                                                                                                      <  1:425915/251‑1 (MQ=255)
                                                   gCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCGCAAAGACGCGCAt                                                                                                                                                                                                <  1:396707/250‑1 (MQ=255)
                                                         cGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTc                                                                                                                                                                                          <  1:422612/250‑1 (MQ=255)
                                                          gCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCCAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGc                                                                                                                                                                                        >  1:320516/1‑251 (MQ=255)
                                                          gCCAATGCCAAACGGCTTACTAAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCCGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGc                                                                                                                                                                                        >  1:385358/1‑251 (MQ=255)
                                                                        gCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTg                                                                                                                                                                          <  2:141776/251‑1 (MQ=255)
                                                                                 gATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGACGTGCGCAGAAAg                                                                                                                                                                 <  1:466023/251‑1 (MQ=255)
                                                                                    aaTCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGcg                                                                                                                                                               <  1:282669/250‑1 (MQ=255)
                                                                                      tcatcCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCCCCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGCCAGGCCGAGATCGAGAATAATCAACTCTGGCTTACAGGTTGCCGCTTCCAGCACGCCGCGTTCCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCACAAAAAGcgacc                                                                                                                                                            >  2:530679/1‑250 (MQ=255)
                                                                                          ccGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATa                                                                                                                                                        <  2:235760/251‑1 (MQ=255)
                                                                                           cGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAg                                                                                                                                                       <  2:487676/251‑1 (MQ=255)
                                                                                            gCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGc                                                                                                                                                      >  1:512092/1‑251 (MQ=255)
                                                                                                cGGCAGCCAGCGCGGCGAGTGTGTCGCTCTCTCCGCTGCTTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAACCCCATCGCCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGGTGCCGCGTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAAGAGcccgt                                                                                                                                                  <  2:205731/251‑4 (MQ=255)
                                                                                                      ccAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCt                                                                                                                                            >  2:121307/1‑251 (MQ=255)
                                                                                                        aGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTc                                                                                                                                          <  1:172778/251‑1 (MQ=255)
                                                                                                             ggcgATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTcaaca                                                                                                                                      >  2:386289/1‑250 (MQ=255)
                                                                                                                 aTTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAg                                                                                                                                 <  1:157148/251‑1 (MQ=255)
                                                                                                                    ttGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGaa                                                                                                                               <  1:458169/250‑1 (MQ=255)
                                                                                                                            ctctTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCTAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCAGCTTCAACAATCAGAACGTTTGt                                                                                                                        <  1:383363/249‑1 (MQ=255)
                                                                                                                              ctTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCg                                                                                                                                                                                                                 <  1:16534/158‑1 (MQ=255)
                                                                                                                              ctTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCg                                                                                                                                                                                                                 >  2:16534/1‑158 (MQ=255)
                                                                                                                                     gcgtgcgGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTGCATCTTCAACAATCAGAACGTTTGTCACATATCCt                                                                                                              <  2:284495/250‑1 (MQ=255)
                                                                                                                                      cgtgcgGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGGCGCGAATAAACTCAATCCCATCACCATCTGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGGCGCCCTCCAGCGCCGTTCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAAGCAGAACGTTGGTCACATAGCCTCa                                                                                                            <  2:52311/251‑1 (MQ=255)
                                                                                                                                           ggAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGaaa                                                                                                       <  1:389716/251‑1 (MQ=255)
                                                                                                                                             aaaGCACAATCACCGGCACCGCGCTCCCCTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGCCGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAAtt                                                                                                     >  1:452468/1‑251 (MQ=255)
                                                                                                                                                              aCCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGc                                                                                     >  1:501486/1‑250 (MQ=255)
                                                                                                                                                               ccGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGc                                                                                     >  1:419745/1‑249 (MQ=255)
                                                                                                                                                                cgcgCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGt                                                                                  <  2:83606/251‑1 (MQ=255)
                                                                                                                                                                      ccACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTg                                                                            <  2:45275/251‑1 (MQ=255)
                                                                                                                                                                           ggCGCAGGTCGCGGATAAACTCAATCCCCTCTCCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAGCGTGTGTCTCATATCCTCATGAAATTCTTCAAGTGCAGTGGCAGTTTGCTgggg                                                                         <  2:371121/249‑1 (MQ=255)
                                                                                                                                                                                 ggTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCGGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTaa                                                                 <  2:126049/251‑1 (MQ=255)
                                                                                                                                                                                        gATAAACTCAATCCCATCTCCGTCGGGCAGGCCGAGATCTAGAATAATCAAAGCTTGCTTACAGGTTGCTGCTTTCAGCAAGCCGCGTTGCAGCGTTTCGGCCCCAAAGACGCGCATGCCGTCGCCCGCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCATCAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGaaa                                                          <  1:88955/251‑1 (MQ=255)
                                                                                                                                                                                                       aTCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTc                                           >  1:366887/1‑251 (MQ=255)
                                                                                                                                                                                                         cACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTcc                                          <  1:211028/250‑1 (MQ=255)
                                                                                                                                                                                                            cATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGtt                                                                                            <  2:156094/197‑1 (MQ=255)
                                                                                                                                                                                                            cATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGtt                                                                                            >  1:156094/1‑197 (MQ=255)
                                                                                                                                                                                                                         gagaTCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGt                         <  1:524952/251‑1 (MQ=255)
                                                                                                                                                                                                                                    taatCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGcccc              >  2:88031/1‑251 (MQ=255)
                                                                                                                                                                                                                                      atCAAGTCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCg            >  2:370137/1‑251 (MQ=255)
                                                                                                                                                                                                                                      atCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCg            <  1:277497/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aaaTCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCTTCAAAGACACGCATCCCGTCGCCATCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGAGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAGCGCGGTAATAGTGCCCCCgtgt         >  2:511085/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ggCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCGTGTACagc    <  2:452468/250‑3 (MQ=255)
                                                                                                                                                                                                                                               ggCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCGTGTACATcc   >  1:318669/1‑251 (MQ=255)
                                                                                                                                                                                                                                                gCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGTCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCGTGTACATCCa  <  1:127446/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
CGGGCGCGGTGGTGGCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCGTGTACATCCA  >  NC_000913/721358‑721848

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: