Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,130,161	G→A	L19F (CTT→TTT)	yghO ←	putative DNA‑binding transcriptional regulator YghO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,130,161	0	G	A	100.0%	149.2 / NA	45	L19F (CTT→TTT)	yghO	putative DNA‑binding transcriptional regulator YghO
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/19);  total (26/19)

CTTTGCCGTAACGCTCGCGGTGCAAGGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCG  >  NC_000913/3129918‑3130368                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                cTTTGCCGTAACGCTCGCGGTGCAAGGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTg                                                                                                                                                                                                           >  2:487551/1‑250 (MQ=255)    tGCCGTAACGCTCGCGGTGCAAGGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGata                                                                                                                                                                                                        <  1:349288/250‑1 (MQ=255)            cgctcgctGTGCAAGGTATCTAGTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTTCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATACCTGATATGAGATCAt                                                                                                                                                                                               <  2:51497/251‑1 (MQ=255)            cGCTCGCGGTGCAAGGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCAt                                                                                                                                                                                               >  1:27450/1‑251 (MQ=255)                     tGCAAGGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCAt                                                                                                                                                                                      >  1:72087/1‑251 (MQ=255)                       cAAGGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTGATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCt                                                                                                                                                                                    <  2:72087/251‑1 (MQ=255)                        aaGGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTg                                                                                                                                                                                   >  1:211114/1‑251 (MQ=255)                          ggTATCTATTTGCGCGGTAATCCGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGAGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTggg                                                                                                                                                                                 >  2:423314/1‑250 (MQ=255)                                 aTTTGCGCGGTAATACGCCCATCTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGcc                                                                                                                                                                              <  1:45108/247‑1 (MQ=255)                                 aTTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATTAATAAATGGTGTATCCAGTTTGGATCGTCGGGATAAAGTGATTACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGAGCATGTTTGTCATCTGGAGCCAt                                                                                                                                                                            <  2:305358/249‑1 (MQ=255)                                         ggTAATACGCCCAACTATCTGCCCCGCTTTTTT‑GCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGtt                                                                                                                                                                 >  1:532173/1‑251 (MQ=255)                                            aaTACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGAGCATGTTTGTCAGCTGGAGCCATAGAACAGGGTtca                                                                                                                                                               <  2:193950/251‑1 (MQ=255)                                              tACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtcatc                                                                                                                                                             >  1:26439/1‑251 (MQ=255)                                                         aTCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGGTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCAGCATGAGtcat                                                                                                                                                    <  2:530990/249‑1 (MQ=255)                                                              ccccGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACtt                                                                                                                                              >  1:565542/1‑250 (MQ=255)                                                                             cacccacGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTg                                                                                                                              <  1:58233/251‑1 (MQ=255)                                                                              acccacGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTg                                                                                                                              >  1:497806/1‑250 (MQ=255)                                                                                ccacGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAAt                                                                                                                           <  1:320155/251‑1 (MQ=255)                                                                                   cGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTc                                                                                                                         >  1:152272/1‑250 (MQ=255)                                                                                            gatgatATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAgc                                                                                                               >  1:403175/1‑251 (MQ=255)                                                                                            gatgatATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAgc                                                                                                               <  2:532173/251‑1 (MQ=255)                                                                                                    ggTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTc                                                                                                         >  2:316963/1‑249 (MQ=255)                                                                                                    ggTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCa                                                                                                        <  1:573037/250‑1 (MQ=255)                                                                                                     gtccgtccCCGGATTTTTCGCAGACAACTGCTCGATGCGTTCAATGAATCACGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGCCGGACAAGCGATAAATCACTTAAGGACATTTTTATTAAGGAAGGTCCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTc                                                                                                         >  1:392046/1‑248 (MQ=21)                                                                                                       ccgtccCCGGATTTTTCGCAGACAAATGATCGTTGCGTTCAATGAATCAAGAGGGTATCCAGTTTGCATCGTCGGGATAAAGAGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGATGCGACCAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGACTTCAGCCGTAAGCGCCGTCAGAcc                                                                                                    >  1:394286/1‑251 (MQ=25)                                                                                                        cgtctcaggATTATTAGCAGACAACTGCTCGTTGCGTTCAATGAATAACGGGGGTATCCAGTTTGGATCGTCGGGACAATGTGCTGACGGAACAGCGGTAAATGCCTGAAGGCCATTTTTATTAAGGAAGGTGCGAACAAATCCCTAATATGAGTTCATGTTTGTCATCTGGAGCCATAGAACAGGCATCATCATGAGTGATCAACTTACATTCGCCGACGGAGAATTCAGCAGTAAGCAACGTCAGAcc                                                                                                    >  2:512252/8‑250 (MQ=18)                                                                                                        cgtccCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCa                                                                                                   <  2:459655/251‑1 (MQ=255)                                                                                                          tccCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGCATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCACTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGtt                                                                                                                                                                 >  1:36497/1‑187 (MQ=255)                                                                                                          tccCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCACTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGtt                                                                                                                                                                 <  2:36497/187‑1 (MQ=255)                                                                                                                      tCGCAGACAAATGCGCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATTTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACGTACCTTCTCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGTAAAGAGAtttt                                                                                      <  2:295311/250‑1 (MQ=255)                                                                                                                              aaaTGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCg                                                                             >  2:339793/1‑251 (MQ=255)                                                                                                                                     cGTTGCGATCAATGAATAAAGGGGGTATACAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATgg                                                                        <  2:127528/249‑1 (MQ=255)                                                                                                                                                aTGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGACATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGCTTTTCTTGTCCCGCATGGAGCAGATTCTGcc                                                           >  2:423683/1‑251 (MQ=255)                                                                                                                                                   aaTAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATg                                                        >  1:19188/1‑251 (MQ=255)                                                                                                                                                             gTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCAATGGCAAAACAtg                                              >  1:21558/1‑251 (MQ=255)                                                                                                                                                                 ccAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACAtggtg                                           <  2:69349/250‑1 (MQ=255)                                                                                                                                                                     tttGGATCGTCGGGATAAAGTGATGACGGAAAAGCGCTCAATGCCTTACGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCGAAACATGGTGGAAGt                                      >  2:493987/1‑251 (MQ=255)                                                                                                                                                                                   aTAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAACAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGtttt                        >  1:189724/1‑251 (MQ=255)                                                                                                                                                                                    tAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGtttt                        >  2:577436/1‑250 (MQ=255)                                                                                                                                                                                      aaGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTAca                     <  2:152272/251‑2 (MQ=255)                                                                                                                                                                                      aaGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATGAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGAAAAACATGGGGGAAGTCATCGAGCCGTTTTAcc                     <  2:497806/251‑1 (MQ=255)                                                                                                                                                                                           atgaCGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCa                  >  2:124189/1‑249 (MQ=255)                                                                                                                                                                                                  aaaaGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTa         <  1:487551/251‑1 (MQ=255)                                                                                                                                                                                                       cGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGc    >  1:150582/1‑251 (MQ=255)                                                                                                                                                                                                         aTAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAATCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCg  >  2:154307/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                CTTTGCCGTAACGCTCGCGGTGCAAGGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCG  >  NC_000913/3129918‑3130368

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