Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,818	G→C	G214A (GGC→GCC)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,818	0	G	C	100.0%	139.3 / NA	47	G214A (GGC→GCC)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (17/30);  total (17/30)

CGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGG  >  NC_000913/4638573‑4639060                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   cgcgCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCCTATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGGATATGGCCCTGAATTTTAGCGGTACCGGCGATCTTTCTTTGGGGCCTgccgggg                                                                                                                                                                                                                                               >  2:155259/1‑250 (MQ=255)         atatCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGc                                                                                                                                                                                                                                       >  1:236351/1‑251 (MQ=255)           atcatcTTAGGCAAGCCCTTTATTGTGAGTAGCGTCGGAGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGTTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGa                                                                                                                                                                                                                                     <  2:170958/251‑1 (MQ=255)                cTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATg                                                                                                                                                                                                                                 >  1:160885/1‑250 (MQ=255)                       aaGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTCTTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGCCCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATTTCCCTTTACCCGAAGGGGCCTGGCGGCAGCCGAAACTGACGCTGAATGTGGCCCTGAATTTAAGCGGTACCGGGGGTCTTTCTTTGGTGCCTGCCGGGCGGAATAGCGAAATGACCTTAAc                                                                                                                                                                                                                         >  2:611817/1‑251 (MQ=255)                              ttATTGTGATTAGCGGCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCaa                                                                                                                                                                                                                   <  1:354458/250‑1 (MQ=255)                                aTTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTg                                                                                                                                                                                                                <  1:507680/251‑1 (MQ=255)                                  tgtgATTAGCGTCGTGGAGGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCGGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCATCTgg                                                                                                                                                                                                               <  2:409751/250‑1 (MQ=255)                                         aGCGTCGTGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGGTACCGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGTGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCGGAATTTAAGCGGTACCGGCGATCTTTCGGTGGTGCCTGCCGGTCGTAATAGCGAAATGACCTTAACCAGCAGCTTGCCGCATcc                                                                                                                                                                                                       <  2:582033/251‑1 (MQ=255)                                            gTCGGGGATGCGCGTGGTATTGGTGTGGTGAACGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGCATATGGCCCTGAATTTAAGCGGTACCGGCGCTCTTTCTGTGGTGCATGCCGGGCGTAATCGCGAAATGACCTTAACCAGCAACTGTCCGcctcac                                                                                                                                                                                                      >  1:329079/1‑244 (MQ=255)                                              cGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTTCCTGCCGGGCGGAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGtt                                                                                                                                                                                                  >  2:28416/1‑251 (MQ=255)                                                 ggATGCGCGTGGTATTTGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTGAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGTACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGTGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGttttt                                                                                                                                                                                               <  2:587707/251‑1 (MQ=255)                                                                                                 ggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTc                                                                                                                                               <  2:196988/251‑1 (MQ=255)                                                                                                        aCCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGgcgc                                                                                                                                         >  1:470130/1‑250 (MQ=255)                                                                                                              gAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGAGCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTGTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCa                                                                                                                                  <  1:572509/251‑1 (MQ=255)                                                                                                                aCCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCaaa                                                                                                                                <  2:337832/251‑1 (MQ=255)                                                                                                                 cccGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCTCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAag                                                                                                                               <  1:611817/251‑1 (MQ=255)                                                                                                                   cGGCACCGGGGTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTGTACCCGAATGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAACGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCAGCATCCCAGGTGGTTAGGTGATTTTTTGCCAGCCAAACTGAAAGGTAGCGAGTCAGGTTGTTAGGCGCAATGGCAAagca                                                                                                                             <  2:229981/251‑1 (MQ=255)                                                                                                                    ggCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTGTCAGGCGCAATGGCAAagcag                                                                                                                            <  1:213149/251‑1 (MQ=255)                                                                                                                           ggTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGttt                                                                                                                     <  2:535000/251‑1 (MQ=255)                                                                                                                           ggTTAGAGCAAGGCGGGCAGGGCGTGCATAGCCCTTTACCCGAAGGGGACGGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTTGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTGTTTAGGTGATGTTTTGCCAGCCAAACGGGAAGTTATCGAGTCAGGTGGGCAGGCGCAATGGCAAAGCAACTGGttt                                                                                                                     <  2:362060/251‑1 (MQ=255)                                                                                                                            gTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTg                                                                                                                    <  1:370931/251‑1 (MQ=255)                                                                                                                                gagCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGGTATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCta                                                                                                                 >  2:98825/1‑250 (MQ=255)                                                                                                                                    aaGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTaataat                                                                                                            >  1:274633/1‑251 (MQ=255)                                                                                                                                         gggCAGGGCGTGCATATCCCTTTACCCGCAGGGGACTGGCGGAAGCAGAACCTGAAGCTGCATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCCGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCgg                                                                                                       >  2:74660/1‑251 (MQ=255)                                                                                                                                                              ttACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGc                                                                                    <  1:174587/249‑1 (MQ=255)                                                                                                                                                                     aaGGGGACTGGCGGAAGCGGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATAcc                                                                           <  2:130359/251‑1 (MQ=255)                                                                                                                                                                     aaGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATAcc                                                                           <  1:180379/251‑1 (MQ=255)                                                                                                                                                                      aGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCg                                                                          <  2:470130/251‑1 (MQ=255)                                                                                                                                                                           aCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTgg                                                                     >  1:6935/1‑251 (MQ=255)                                                                                                                                                                            cTGGCGGAAGCAGAACCTGAAGCTGAATCTGGCCCTGAATTTACGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGCAATGACCTTACCCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTggg                                                                    >  2:320696/1‑251 (MQ=255)                                                                                                                                                                                 ggAAGCCGAACCTGAAACTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTACTAGCGCAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTACTAATCTCGGGGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAAc                                                               >  2:435282/1‑251 (MQ=255)                                                                                                                                                                                   aaGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACtt                                                             <  1:93339/251‑1 (MQ=255)                                                                                                                                                                                        gAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTccccg                                                        <  1:28416/251‑3 (MQ=255)                                                                                                                                                                                          aCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTT‑CAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCTGCTGGGAAAACTTCCCCGCGt                                                     >  1:350631/1‑251 (MQ=255)                                                                                                                                                                                               aaGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAg                                                 <  1:142449/251‑1 (MQ=255)                                                                                                                                                                                                                 aaTTTAAGCGATACCGGCGATCTTACTGGGGGGCCTGCCGGGAGTACTAGCGAAATGACCTTAACCAGCAATTGGACGCATCCCAGTTTTTTAGGCGATTTTTTGCCAGCCCAACGGGAAGTTAGCGAGTCCGGTTTTCACGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGAGTGCGTTTTGCTTCAGGCAATGCTACCGGCTGGGCAAACTTCCCAGCGTTTCGCGTCCCAGtggcggcgcc                               >  2:553631/1‑242 (MQ=255)                                                                                                                                                                                                                  aTTTAAGCGGTACCGGTGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTGTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTGCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCa                              <  2:431291/251‑1 (MQ=255)                                                                                                                                                                                                                      aaGCTGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGTCCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTc                                                                                                                                                       <  2:463908/126‑1 (MQ=255)                                                                                                                                                                                                                      aaGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCCTCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTc                                                                                                                                                       >  1:463908/1‑126 (MQ=255)                                                                                                                                                                                                                      aaGCGGTACCGGCGATCTTTCTGGGGGGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACGGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGcc                           <  2:274633/250‑1 (MQ=255)                                                                                                                                                                                                                          ggTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACgg                                                                                                                                                                     <  2:164361/108‑1 (MQ=255)                                                                                                                                                                                                                          ggTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACgg                                                                                                                                                                     >  1:164361/1‑108 (MQ=255)                                                                                                                                                                                                                                  cGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCCGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAACCTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAAt              >  2:593921/1‑251 (MQ=255)                                                                                                                                                                                                                                      atgtCTGTGGTGCCTGCCGTGCGTAATAGCGAAATGACCGTAACCAGCAACTGGACGCATCCCATTTTTTTAGGTGATGTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGTCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGACACCGGCTGGGAAAACTTCCCCGCGTTTAGGTTCGCAGTGACGACGCCAGCCGATCAAGACCAATTaa           <  2:450811/247‑1 (MQ=255)                                                                                                                                                                                                                                        ttCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTaa           <  1:593921/248‑1 (MQ=255)                                                                                                                                                                                                                                            gtggtgCCTGCCGGGCGTAATAGCGAAAAGGCCTTACCCAGCAACTGGCCGCATCCCAGTTTTTTAGTGGATTTTTTGCCGGCCAAGCGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCg    <  2:135426/251‑1 (MQ=255)                                                                                                                                                                                                                                              ggtgCCTGCCGGGCGTAATAGCGAAAGGACCTTAACCAGCAACGGGCCGCATCCCAGTTTTGTAGGTGATTTTTTGCCAGCCAAACGGTAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAAAACCAATTAACGGACCggg  <  1:435282/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   CGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGG  >  NC_000913/4638573‑4639060

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