Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,506,393	T→C	D315G (GAC→GGC)	yhfT ←	uncharacterized protein YhfT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,506,393	0	T	C	100.0%	239.2 / NA	72	D315G (GAC→GGC)	yhfT	uncharacterized protein YhfT
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/36);  total (36/36)

GAGGAAGTAAATGGCAACCGCGATAGAGAATCCGGTATAACCCGCCATCTTAATTGCTGCGAAAATCGAACCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGTCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCGGCGTTACGCCTGCGGAATATGCTTTCTCCAGTGTGAAGATCGACACTTCACTGCCAGCAAAAATCTTCATGCTGGCAACGGCGGCAATCAATGCTCCCACGATGG  >  NC_000913/3506146‑3506639                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       gAGGAAGTAAATGGCAACCGCGATAGAGAATCCGGTATAACCCGCCATCTTAATTGCTGCGAAAATCGAACCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCa                                                                                                                                                                                                                                                     <  2:344281/251‑1 (MQ=255) gAGGAAGTAAATGGCAACCGCGATAGAGAATCCGGTATAACCCGCCATCTTAATTGCTGCGAAAATCGAACCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCa                                                                                                                                                                                                                                                     <  1:401160/251‑1 (MQ=255)    gAAGTAAATGGCAACCGCGATAGAGAATCCGGTATAACCCGCCATCTTAATTGCTGCGAAAATCGAACCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCg                                                                                                                                                                                                                                                  >  1:389654/1‑251 (MQ=255)     aaGTAAATGCCAAGCTCGATAGAGAATCCGGTATAACCCGCCATCTTAATTGCTGCGAAAATCGAACCGACCAGCAGCGCCACTTCCCTCAGCATATTCATGGCGTTACGGCTGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATCACCACTGCGCCTAATACCGCCGCCACCATCGGATTCGGCGAGAGATAGCCCCCCGc                                                                                                                                                                                                                                                 >  2:259999/1‑251 (MQ=255)        tAAATGGCAACCGCGATAGAGAATCCGGTATAACCCGCCATCTTAATTGCTGCGAAAATCGAACCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGCAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCAt                                                                                                                                                                                                                                               >  1:355520/1‑250 (MQ=255)         aaaTGGCAACCGAGATAGAGAAGCCGGTATAACCCGCCATCTTAATTGCTGCGAAAATCGAACCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAATCATCAAGACTTCCGCCGAAATATCCACTGCGCCTAATACCGCTGCAACCATCGGATTCTGCGAGAGATAGCCCACTGCATa                                                                                                                                                                                                                                              <  2:479292/250‑1 (MQ=255)              gCAACCGCGATAGAGAATCCGGTATAACCCGCCATCTTAATTGCTGCGAAAATCGAACCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACaaa                                                                                                                                                                                                                                        <  2:236341/251‑1 (MQ=255)                     cgATAGAGAATCCGGTATAACCCGCCATCTTAATTGCTGCGAAAATCGAACCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTaa                                                                                                                                                                                                                                   >  1:28266/1‑249 (MQ=255)                              aTCCGGTATAACCCGCCATCTTAATGGCTGCGAAAATCGACCCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCAGCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACt                                                                                                                                                                                                                        <  2:644638/251‑1 (MQ=255)                                 cGGTATAACCCGAAATCTTACTTGCTGCGCAAATCGAACCGCCCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCGAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGCGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGc                                                                                                                                                                                                                      >  1:598396/1‑250 (MQ=255)                                  ggTATAACCCGCCATCTTAATTGCTGCGAAAATCGAACCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCAt                                                                                                                                                                                                                    <  1:204191/251‑1 (MQ=255)                                     ataACCCGCCATCTTAATTGCTGCGAAAATCGAACCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAACCAAAGGTAAAGCCCGCAACTGCATac                                                                                                                                                                                                                  >  1:719223/1‑250 (MQ=255)                                       aaCCCTCCATATTAATTGATGAGAAACTAGAACCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCGGAAGCCTTACGCACCGACGGGTAGCGTCCCAGCTATTTGCCGCTCGAACGAAGCAGCAAGACTTCCGCCGCAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCCTACACACAGGAAAAGCCCGCAACTGCATacacc                                                                                                                                                                                                               >  2:702152/1‑250 (MQ=255)                                                 cTTAATTGCTGCGAAAATCGAACCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCGACTGCATACACACCCGTTGCt                                                                                                                                                                                                      >  1:331447/1‑250 (MQ=255)                                                  ttAATTGCTGCGAAAATCGAACCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTa                                                                                                                                                                                                     >  1:236624/1‑250 (MQ=255)                                                      ttGCTGCGAAACTCGAACAGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCCTTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGAAAAGCCCGCAACTGCATACACACCCGTTGCTAAcgcg                                                                                                                                                                                                >  2:416060/1‑251 (MQ=255)                                                          tgcGAAAATCGAACCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGGACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTCATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCAACGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCggtgg                                                                                                                                                                                            >  1:105647/1‑251 (MQ=255)                                                          ggcGAAAAGCGAACCGACCAGCAGCGCCCCTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACTCACCGACGGGTAGCGCCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCGGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAACCAAAGGTAAAGCCCGCAACTGCATACAGACCCGTTGCTAACGCggtgg                                                                                                                                                                                            <  2:127273/250‑1 (MQ=255)                                                           gcGAAAATCGAACCGACCAGCAGCGCCACTTCAATCAGCATATTCATGGCGTTACGGATGTTATCCGAAGCATTACGCACCGACGGGGAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACAGCTGCAACCATCGGATTCGGCGAGAGATAGCCCAACGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCAGTTGCTAACGCggtggg                                                                                                                                                                                           >  1:472052/1‑250 (MQ=255)                                                                 aTCGAACCGACCAGCAGCTCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTCGCGTCCCAGCCATTTGCCGCTCGAACGAAGCAGCAAGACTTCCGCCGAAATACCCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATACACAAAGGTAAAGCCCGCAACTGCATCCACACCCGTTGCTAACGCGGTGGTTGCAAt                                                                                                                                                                                     >  2:239455/1‑251 (MQ=255)                                                                 aTCGAACCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAAt                                                                                                                                                                                     >  1:473555/1‑251 (MQ=255)                                                                    gAACCGCCCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACa                                                                                                                                                                                                                                          <  1:62700/195‑1 (MQ=255)                                                                    gAACCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTCGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACa                                                                                                                                                                                                                                          >  2:62700/1‑195 (MQ=255)                                                                     aaCCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAAc                                                                                                                                                                                 >  2:590839/1‑251 (MQ=255)                                                                         gACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCGACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGca                                                                                                                                                                             <  2:331447/251‑1 (MQ=255)                                                                             agcagcGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACaa                                                                                                                                                                          >  1:549092/1‑250 (MQ=255)                                                                              ggagcGCCACTTCCATCAGCAGAGTCAGGGCGTTACGGATGTTATCCGTCGCATTACGCACCGACTGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCAGCGGATTCGGCGAGAGATAGCTCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTGGCTAACGCGATTGTGGTAATCAACgtctcaaat                                                                                                                                                                        <  2:423836/249‑9 (MQ=255)                                                                                      cctcccATCAGCATATTCATGGCGTGACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAATGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTcc                                                                                                                                                                <  2:93866/247‑1 (MQ=255)                                                                                      aCTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTcc                                                                                                                                                                >  2:21522/1‑251 (MQ=255)                                                                                         tCCATCAGCATATTCATGGCGTTACGGATGTTCTCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAACTAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCTGTTGCTACCACGGTGGTGGCAATCACCGGAACAAAACCCAGTCcgcg                                                                                                                                                             >  2:362768/1‑251 (MQ=255)                                                                                          ccATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCgcgc                                                                                                                                                            >  2:19943/1‑251 (MQ=255)                                                                                            aTCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCAt                                                                                                                                                          >  1:66812/1‑251 (MQ=255)                                                                                              cAGCATATTCATCGCGTTACGGATGTTATCCGACGCATTACGCACCGACAGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCACCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGGGGCAATCAACGGCACAAAACCCAGTCCGCGCATa                                                                                                                                                         >  2:452167/1‑250 (MQ=255)                                                                                                   tatTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGCAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCCACTGCATACACACCCGTTGCTa                                                                                                                                                                                                     >  1:643346/1‑201 (MQ=255)                                                                                                   tatTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATATCCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCTTTGCTa                                                                                                                                                                                                     <  2:643346/201‑1 (MQ=255)                                                                                                   tatTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCt                                                                                                                                                   >  1:86473/1‑251 (MQ=255)                                                                                                    atTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTg                                                                                                                                                  <  2:389654/251‑1 (MQ=255)                                                                                                    atTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTg                                                                                                                                                  <  1:609147/251‑1 (MQ=255)                                                                                                      tCATGGCGTTACGGATGTTACCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGc                                                                                                                                                 <  2:491141/250‑1 (MQ=255)                                                                                                            cgttaggcaTGTTATCCGACGCATTACGCACCGACGGGCAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCCCGACTTACCCAGAAATAAACAATGAGCCTAATACCGCTACAACAATAGGATTCGGAGAGAGATAGCACACCGCATAAAAAAAGGTACAGCCCGCCACTGCATACACACCAGCTGCTAACGCCGCGGCGGCAATCAACGGCACAAACCCCAGACCGCGCATAAAttaagaccgaccc                                                                                                                                          >  2:487484/9‑240 (MQ=255)                                                                                                            cGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAgcc                                                                                                                                          <  1:9833/251‑1 (MQ=255)                                                                                                              ttACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGAAGCAAGACTTCCGCCGACATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGCGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTACCGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAgccgc                                                                                                                                        >  1:693843/1‑251 (MQ=255)                                                                                                                  ggATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGa                                                                                                                                    <  2:236624/251‑1 (MQ=255)                                                                                                                   gATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCATAATACCGATGCAACCATCGGATTCGGCGAGCGATAGCCCCACGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAATTTCTGCCAGAGCCGCCTGAt                                                                                                                                   >  1:40686/1‑251 (MQ=255)                                                                                                                      gTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACCCCCGTTGCTAACGCGGTGGTGGCAATCCACGGCCCAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCCCCTGATTa                                                                                                                                 >  2:413289/1‑250 (MQ=255)                                                                                                                         aTCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATc                                                                                                                              <  2:237229/250‑1 (MQ=255)                                                                                                                         aTCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCa                                                                                                                             <  2:86473/251‑1 (MQ=255)                                                                                                                           ccGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGATCCGCCTGATTAATCAGc                                                                                                                           <  2:486816/251‑1 (MQ=255)                                                                                                                             gACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGt                                                                                                                         <  1:413289/251‑1 (MQ=255)                                                                                                                                 cATTACGCACCGACGGGTAGCGGCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGc                                                                                                                     <  2:530677/251‑1 (MQ=255)                                                                                                                                     aCGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGAtt                                                                                                                 <  2:28266/251‑1 (MQ=255)                                                                                                                                       gCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCATCTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGt                                                                                                               >  1:289631/1‑251 (MQ=255)                                                                                                                                           cGACGGGTAGCGTCCCAGCTATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGTCGAGAGATATCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCg                                                                                                           <  1:702152/251‑1 (MQ=255)                                                                                                                                                 gTAGCGGCCCCGCCATTTGCCGATCGAACGAAGCAGCAACACTGCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGATATAGCCCGCCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACTCGGGGGTGGCAATCAACGGCGCAAGCCCCAGTCCGCGCATAAATTTTGCCAGAGCCGCCTGATTAATCAGCTGTTGCTATTGTTCCgtcggta                                                                                                     <  2:48515/251‑7 (MQ=255)                                                                                                                                                      gTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGGGGCAATCAACGGCACAAACCCCAGCCCGCGCATAAATTCTGCCAGAGCCGCCTGGTTAATCAGCGTTTGCGATTGTTCCGGCGTTACGCCt                                                                                                >  2:406429/1‑251 (MQ=255)                                                                                                                                                             ccATTGGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCTCAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCGGCGTTACGCCTGCGGAat                                                                                         <  1:172609/251‑1 (MQ=255)                                                                                                                                                               atgtGCCGATCGAACGAAGCACAAAGCCTTCCGCCGAAATAACCACTTCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATTCCCACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCGGCGTTACGCCTGCGGAatat                                                                                       <  2:402316/248‑1 (MQ=255)                                                                                                                                                                  tGCCGATCGAACGAAGCAGCAAGACGTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCGGCGTTACGCCTGCGGAATATGc                                                                                     <  2:688378/250‑1 (MQ=255)                                                                                                                                                                       aTCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCCGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCGGCGTTACGCCTGCGGAATATGCTTtctc                                                                               >  1:140632/1‑251 (MQ=255)                                                                                                                                                                         cgaacgaaGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCGGCGTTACGCCTGCGGAATATGCTTTCTCCa                                                                             <  1:203672/251‑1 (MQ=255)                                                                                                                                                                                   agcaAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTACAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCGGCGTTACGCCTGCGGAATATGCTTTCTCCAGTGTGa                                                                       >  1:308708/1‑247 (MQ=255)                                                                                                                                                                                            cccgccgAAATAACCCCTGCGCCTAATACCGCTGCAACCATCCGATTCGGCGAGAGATAGCCCCCCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACCCCCGTTGCTAACGCGGGGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATGCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCGGCGTGACGCCTGCGGAATATTCTTTCTCCAGTGTGAAGATCGACACTTc                                                          <  2:172869/250‑1 (MQ=255)                                                                                                                                                                                             ccgccgAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCGGCGTTACGCCTGCGGAATATGCTTTCTCCAGTGTGAAGATCGACACTTCa                                                         <  2:549092/251‑1 (MQ=255)                                                                                                                                                                                                        aCCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCGGCGTTACGCCTGCGGAATATGCTTTCTCCAGTGTGAAGATCGACACTTCACTGCCAGCaaa                                              <  1:610899/251‑1 (MQ=255)                                                                                                                                                                                                           aCTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCTAGAGCCTCCTGATTAATCAGCGTTTTCGATTGTTCCGGCGTTACGCCGGCGGAATATGCTTTCTCCAGTGTGAAGATCGACACTTCACTGCCAGCAAAAAt                                           <  1:522109/251‑1 (MQ=255)                                                                                                                                                                                                              gcgcCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTGACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCGGCGTTACGCCTGCGGAATATGCTTTCTCCAGTGTGAAGATCGACACTTCACTGCCAGCAAAAATCtt                                        <  1:425808/251‑1 (MQ=255)                                                                                                                                                                                                                           cTGCAACCATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCGGCGTTACGCCTGCGGAATATGCTTTCTCCAGTGTGAAGATCGACACTTCACTGCCAGCAAAAATCTTCATGCTGGCAAcg                           >  2:253868/1‑251 (MQ=255)                                                                                                                                                                                                                               aaCCATCGGATTCGGCGAGAGATAGCCCTCCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGGGTGCGATTGTTCCGGCGTTACGCCTGCGGAATATGCTTTCTCCAGTGTGAAGATCGACACTTCACTGCCAGCAAGAATCTTCATGCTGGCAAcggcgg                       <  2:355520/251‑1 (MQ=255)                                                                                                                                                                                                                                  cATCGGATTCGGCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCGGCGTTACGCCTGCGGAATATGCTTTCTCCAGTGTGAAGATCGACACTTCACTGCCAGCACAAATCTTCATGCTGGCAACGGCGGca                     >  2:56827/1‑250 (MQ=255)                                                                                                                                                                                                                                             gCGAGAGATAGCCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCGGCGTTACGCCTGCGGAATATGCTTTCTCCAGTGTGAAGATCGACACTTCACTGCCAGCAAAAATCTTCATGCTGGCAACGGCGGCAATCAATGCTccc         >  2:141808/1‑251 (MQ=255)                                                                                                                                                                                                                                                  agaTAGCCCACCGCATAAACAAAGGTAAAGACCGCAACTGCATACCCACCCGTTGCTAACGCGGTGGGGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCGGCGTTACGCCTGCGGAATATGCTTTCTCCAGTGTGAAGATCGACACTTCACTGCCAGCAAAAATCTTCATGCTGGCAACGGCGGCAATCAATGCTCCCACGAt    >  1:132763/1‑251 (MQ=255)                                                                                                                                                                                                                                                    aTAGCCCACCGCAGAAACAAAGGTAAAGCCCGCAACGGCAGACACCCCCGTTGCTAACGCGGTGGTGGCAATCAACGTCACAAACCCCAGGCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCGGCGTTACGCCGGCGGAATATGCTTTCTCCAGTGTTAAGATCGACACTGCACTGCCAGCAAAAAGCTGCATGCTGGCAACGGCGGcaagcagggcttccatgatgg  <  2:693843/251‑19 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       GAGGAAGTAAATGGCAACCGCGATAGAGAATCCGGTATAACCCGCCATCTTAATTGCTGCGAAAATCGAACCGACCAGCAGCGCCACTTCCATCAGCATATTCATGGCGTTACGGATGTTATCCGACGCATTACGCACCGACGGGTAGCGTCCCAGCCATTTGCCGATCGAACGAAGCAGCAAGACTTCCGCCGAAATAACCACTGCGCCTAATACCGCTGCAACCATCGGATTCGGCGAGAGATAGTCCACCGCATAAACAAAGGTAAAGCCCGCAACTGCATACACACCCGTTGCTAACGCGGTGGTGGCAATCAACGGCACAAACCCCAGTCCGCGCATAAATTCTGCCAGAGCCGCCTGATTAATCAGCGTTTGCGATTGTTCCGGCGTTACGCCTGCGGAATATGCTTTCTCCAGTGTGAAGATCGACACTTCACTGCCAGCAAAAATCTTCATGCTGGCAACGGCGGCAATCAATGCTCCCACGATGG  >  NC_000913/3506146‑3506639

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