Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,123,938 G→A G98G (GGC→GGT glcG ← putative heme‑binding protein GlcG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,123,9380GA97.1% 219.8 / NA 70G98G (GGC→GGTglcGputative heme‑binding protein GlcG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (44/23);  minor base C (0/2);  total (44/26)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCCCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGC  >  NC_000913/3123694‑3124183
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
ggcgcCAGAACGCCGCAGCGTCCAGCCCTGTTCCCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCCCCAATAATTTGCCCATCCACAACATCCGGTACTCCGCCt                                                                                                                                                                                                                                                 <  2:380190/250‑1 (MQ=255)
  cgcCAGAACGCCGCAGCGTCCAGCCCTGTTCCCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTc                                                                                                                                                                                                                                               >  2:114472/1‑251 (MQ=255)
   gcCAGAACGCCGCAGCGTCCAGCCCTGTTCCCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTcc                                                                                                                                                                                                                                              >  1:163757/1‑251 (MQ=255)
   gcCAGAACGCCGCAGCGTCCAGCCCTGTTCCCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTcc                                                                                                                                                                                                                                              <  1:349002/251‑1 (MQ=255)
    ccAGAACGCCGCAGCGTCCAGCCCTGTTCCCGGTAAAACTTCTGCATCCACAAATCGTTGAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGTCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCa                                                                                                                                                                                                                                             <  2:319830/251‑1 (MQ=255)
     cAGAACGCCGCAGCGTCCAGCCCTGTTCCCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGAGACATGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACAAATAATTTGACCATCCACAACAACCGGTACACCGCCTTcc                                                                                                                                                                                                                                              >  2:282515/1‑249 (MQ=255)
       gAACGCCGCAGCGGCCAGCCCTGTTCCCGGTAAAACTTCTTCATCCACGAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTATGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCAGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACACCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACAGCCTTCCAGCg                                                                                                                                                                                                                                          >  2:451331/1‑251 (MQ=255)
            ccgcAGCGTCCAGCCCTGTTCCCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGtt                                                                                                                                                                                                                                     >  1:588081/1‑251 (MQ=255)
             cgcAGCGTCCAGCCCTGTTCCCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGCCGtt                                                                                                                                                                                                                                     >  1:417044/1‑250 (MQ=255)
             cgcAGCGTCCAGCACTGTTCCCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTATGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTCTTTCGCCACCACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAACCCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACAGCCTTCCAGCGACGtt                                                                                                                                                                                                                                     >  2:85960/1‑250 (MQ=255)
                  cGTCCAGCCCTGTTCCCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAAc                                                                                                                                                                                                                               <  1:347167/251‑1 (MQ=255)
                     cccgCCCTGTTTCCGGTAAACCTTCTTCATCCACCAAACGTTCAAAATTCGCGTCAATGCGTAACCGGCTCTGGGTTATGGTTTGACTAATTGTTTATATCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGAGCATCCTGCGCTCCGGTTAAACCAGAATCCCCCACGGCACCAATAATTCGCCCATCAACAACAACCGGTACACCGCCTTCCAGCGACGTGAATAACgac                                                                                                                                                                                                                            >  2:5478/4‑249 (MQ=255)
                       aGCCCTGTTCCCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGgcgc                                                                                                                                                                                                                          <  1:273461/251‑1 (MQ=255)
                             gTTCCCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCAACAATAATTTGCCCATCCACAACAACCGGTACCCCGCCTTCCAGCGACGTTAATAACGGCGCAGTc                                                                                                                                                                                                                      <  1:491690/249‑1 (MQ=255)
                                         aCTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACg                                                                                                                                                                                                        >  1:481497/1‑251 (MQ=255)
                                                   ccACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTc                                                                                                                                                                                              >  1:147191/1‑251 (MQ=255)
                                                        aaaCGTTTAAAATTAGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCAt                                                                                                                                                                                         >  1:490532/1‑251 (MQ=255)
                                                           cGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCAtctc                                                                                                                                                                                      >  1:252284/1‑251 (MQ=255)
                                                                  aaTTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAg                                                                                                                                                                               <  1:366944/251‑1 (MQ=255)
                                                                     tGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGccc                                                                                                                                                                            >  2:103905/1‑251 (MQ=255)
                                                                         gTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAg                                                                                                                                                                        <  1:103905/251‑1 (MQ=255)
                                                                           cAATGCGTAAACGGCTCTGGTTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCTCGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGt                                                                                                                                                                       <  2:541587/250‑1 (MQ=255)
                                                                           cAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGtt                                                                                                                                                                      >  2:483208/1‑251 (MQ=255)
                                                                            aaTGCGTAAACGGCCCTGGGTTATGGTTTGACCCATTGTTTACCTCCTCGTTTTCGCTTATTTCGCCCACACCGCTTCTTCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTCAACCAGAATCGCCCACGGCCCCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGGCGTTTATGCCGTCGCAGTCACGAACGCGGTTCGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGtgt                                                                                                                                                                     <  2:652953/251‑3 (MQ=255)
                                                                                   aaaCGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCacgacg                                                                                                                                                              >  1:622995/1‑251 (MQ=255)
                                                                                         ctctGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGcccc                                                                                                                                                          >  1:584567/1‑249 (MQ=255)
                                                                                          tctGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCAAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGcacc                                                                                                                                                          >  1:274182/1‑245 (MQ=255)
                                                                                             gggTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAgcgc                                                                                                                                                     <  1:288170/250‑1 (MQ=255)
                                                                                             gggTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAgcg                                                                                                                                                      <  2:638481/249‑1 (MQ=255)
                                                                                             gggTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGATGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAACCCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAgcg                                                                                                                                                      <  1:169606/249‑1 (MQ=255)
                                                                                                 tATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGcggcgg                                                                                                                                                >  1:208076/1‑251 (MQ=255)
                                                                                                   tGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGt                                                                                                                                               >  1:581623/1‑250 (MQ=255)
                                                                                                             cATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCttt                                                                                                                                      <  2:236423/249‑1 (MQ=255)
                                                                                                             cATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTtct                                                                                                                                    <  2:15253/251‑1 (MQ=255)
                                                                                                              aTTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTACGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGCTTCACGACGCCCCAGCGCGGCGGTACGCGCTTtctc                                                                                                                                   >  2:380384/1‑251 (MQ=255)
                                                                                                                tGTTTATATCCTCGTTTACGCTTATTTCGCCAACACCGCTGCAGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGCAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCt                                                                                                                                 >  2:409100/1‑251 (MQ=255)
                                                                                                                         cctcGTTTTCGCTTATTTCGCCAACGCCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGgaga                                                                                                                           <  1:411905/248‑1 (MQ=255)
                                                                                                                         cctcGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGatat                                                                                                                        >  1:31586/1‑251 (MQ=255)
                                                                                                                            cGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGAtata                                                                                                                       <  2:31586/249‑1 (MQ=255)
                                                                                                                             gTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGc                                                                                                                    >  1:300026/1‑251 (MQ=255)
                                                                                                                               tttCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTACAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTAAGCGCTTTCTCCTGGGAGATATAAGCCg                                                                                                                  >  2:368875/1‑251 (MQ=255)
                                                                                                                                ttCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGc                                                                                                                 >  1:32957/1‑251 (MQ=255)
                                                                                                                                      tATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGCGATATAAGCCGCAATCgg                                                                                                           >  2:556504/1‑251 (MQ=255)
                                                                                                                                       aTTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGgc                                                                                                          <  1:477930/251‑1 (MQ=255)
                                                                                                                                                acacCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCg                                                                                                 >  1:6488/1‑251 (MQ=255)
                                                                                                                                                 cacCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGt                                                                                                <  1:97098/251‑1 (MQ=255)
                                                                                                                                                   ccGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTcc                                                                                              >  1:106813/1‑251 (MQ=255)
                                                                                                                                                     gctgctGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCAt                                                                                            >  2:78387/1‑251 (MQ=255)
                                                                                                                                                      ctgctgCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATg                                                                                           >  1:395648/1‑251 (MQ=255)
                                                                                                                                                        gctgcCGATTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTCCTCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATGCg                                                                                         >  1:174025/1‑251 (MQ=255)
                                                                                                                                                             cGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCATTCTCCTGGGAGATCTAATCCGCAATCGGCGCGCAATCGTCCATGCGACtt                                                                                     >  2:633952/1‑250 (MQ=255)
                                                                                                                                                               cTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTGCGCGCTTTCTCCTGGGAGATATAAGCCGAAATCGTCGCGCAATCGTCCATGCGACTTaa                                                                                   >  1:282139/1‑250 (MQ=255)
                                                                                                                                                                   cgcgACCTGCGAATCCTGTGCTCCGGTTAAACCAGTAACGCCCACGGCACCAATACTTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATATCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTGCATGCGCCTTAACGccc                                                                              >  2:640340/1‑250 (MQ=255)
                                                                                                                                                                            cgcATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATGCGACTTACCGCCAGCAGATGAc                                                                     >  2:445946/1‑251 (MQ=255)
                                                                                                                                                                              cATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGAccg                                                                   >  1:440119/1‑251 (MQ=255)
                                                                                                                                                                               aTCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGAccgc                                                                  >  2:519011/1‑251 (MQ=255)
                                                                                                                                                                               aTCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGAccgc                                                                  >  2:542826/1‑251 (MQ=255)
                                                                                                                                                                                     tgCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCAt                                                            >  1:459808/1‑251 (MQ=255)
                                                                                                                                                                                               aaaCCAGAAACGCCCACGGCACCAAGAATTTGCCCATCCACAACAACCGGTCCACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAAGCGGCGCGCAATCGGCCATGCGACTTAACGCCAGCAGATGCCCGCCGTCATCGGCAACAGc                                                  <  2:191650/251‑1 (MQ=255)
                                                                                                                                                                                                aaCCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCCTCGGCAACAGcc                                                 >  2:25623/1‑250 (MQ=255)
                                                                                                                                                                                                aaCCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGc                                                  >  1:483699/1‑250 (MQ=255)
                                                                                                                                                                                                     gAAACGCCCACGGCACCAATAATTTACCCATCCACAACAACAGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCCGTACGTCCGCTGTCCCCCACATCTTCATAGCCCTTAGTTTCACGACGCCCAAGCGCGGCGGTACGAGACTCCCCCTGGGAGACATAAGCCGCAATCGGCGCGCAATCGCCCTTGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAgaactcg                                             >  2:470612/1‑244 (MQ=255)
                                                                                                                                                                                                      aaaCGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGa                                            >  2:17441/1‑250 (MQ=255)
                                                                                                                                                                                                             cACGGCACCAATAATTTGCCCAGCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGAcc                                    <  1:110572/251‑1 (MQ=255)
                                                                                                                                                                                                                 gCACCAATAATTGTCCCATCCGCAACAACCGGTACCCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTATTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCACCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACCGACCAGtt                                <  2:440119/251‑1 (MQ=255)
                                                                                                                                                                                                                       agaatTTGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCGTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGAGCAGTTAttttt                          <  2:417044/249‑1 (MQ=255)
                                                                                                                                                                                                                             tGCCCATCCACAACAACCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTgcg                     <  1:300539/250‑1 (MQ=255)
                                                                                                                                                                                                                              ccccATCCAAAACACCCGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGTTTCGTCCTTTGTTCACCATCTCTTCATAGCCCTTAGTTTCCCGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATAGAAGCCGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGcc                   <  2:198353/250‑1 (MQ=255)
                                                                                                                                                                                                                                             ccGGTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCg    <  1:542826/251‑1 (MQ=255)
                                                                                                                                                                                                                                                gTACACCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACCGACCAGTTATTTTTCTGCGCCTCTTCCTTACCTGCGGc  >  2:225292/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
TGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCCCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGCGCAATCGTCCATGCGACTTAACGCCAGCAGATGACCGCCGTCATCGGCAACAGCAATGGAAACAGACCAGTTATTTTTCTGCGCCTCTTCCTGACCTGCGGC  >  NC_000913/3123694‑3124183

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