Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,430,584	C→T	P241S (CCA→TCA)	cycA →	D‑serine/alanine/glycine/:H(+)symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,430,584	0	C	T	100.0%	124.5 / NA	42	P241S (CCA→TCA)	cycA	D‑serine/alanine/glycine/:H(+)symporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/20);  total (22/20)

AGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGT  >  NC_000913/4430345‑4430830                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       aGTTCTGGTTTGCGATGATCAAAAGCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCTGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGaaagca                                                                                                                                                                                                                                             <  2:184301/251‑4 (MQ=255)    tCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGGCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACt                                                                                                                                                                                                                                           <  2:225807/250‑1 (MQ=255)       tgTTTGCTATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGc                                                                                                                                                                                                                                         <  2:206785/248‑1 (MQ=255)               atgatCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGAt                                                                                                                                                                                                                               >  2:100201/1‑251 (MQ=255)                 gatCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTa                                                                                                                                                                                                                             <  1:100201/251‑1 (MQ=255)                 gatCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTa                                                                                                                                                                                                                             <  1:257602/251‑1 (MQ=255)                                 aTTGTGTCGATGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCg                                                                                                                                                                                                             >  1:60218/1‑251 (MQ=255)                                    gtgtCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTAt                                                                                                                                                                                                          >  1:40985/1‑251 (MQ=255)                                             aTTGTCATAGGCCTGGTCATGGAGGCGATGCACTATCAGTCACAGACAGGAATGGAAGCGTCATTCGCGCATTTGTGGAATGACAGCAGCTGGTTCCCGAAAGGTTAAAGTAGCTTCTTTACCGGATTACAGATAGCGGTTTTCGCATTCGTGGGGATTCAGCTGGTAGGTACACCAGCTGCGGAAACCAAAGATTGAGAGAAATCACTGACCCGCGCGATTAACTCCATTTCGCTCCGTATCACTATGtt                                                                                                                                                                                                 >  1:233800/1‑251 (MQ=255)                                                      ggCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGCGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGCTTAAGTGGCTTCTTTGCCGGATTCCAGCTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTCGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCtt                                                                                                                                                                                        >  2:146531/1‑251 (MQ=255)                                                       gCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTc                                                                                                                                                                                       <  1:45064/251‑1 (MQ=255)                                                            gTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCTCGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCGTCGCGCt                                                                                                                                                                                  <  2:41380/251‑1 (MQ=255)                                                            gTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTATTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAAATCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCt                                                                                                                                                                                  >  1:165977/1‑251 (MQ=255)                                                                 ggtggCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGAGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGAtt                                                                                                                                                                              <  1:65914/250‑1 (MQ=255)                                                                              tttCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGt                                                                                                                                                                >  2:233419/1‑251 (MQ=255)                                                                                    tCACCGACTGGTGTGGAAGCGTCATTCGAGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAAACAAAGATTCAGAGAAATCACTGCCACGAGCGCTTAACTCCATTCCGCTCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGcc                                                                                                                                                          >  1:142197/1‑251 (MQ=255)                                                                                    tCACCGACTGCTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGGTCCCGGAGGGTTTAAGTGGCTTCTTTGCCGGATTGCAGATAGCGGTTTTCGCTTTCGTGGGGCTTGAGCTGGTAGGTACCACAGCTGCGGAAACCAAAGATTCAGAGAAATCATTGCCACGCGCGATTAACTCCATTCAGATACGTATCAATATGTTCTACGTCTTCGCGATGATTGGGATTATGTcaggtaagcc                                                                                                                                                          >  2:212615/1‑242 (MQ=255)                                                                                       ccGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCAGTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTg                                                                                                                                                       <  2:93068/251‑1 (MQ=255)                                                                                        cGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATGCCATATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTgg                                                                                                                                                      <  1:180318/251‑1 (MQ=255)                                                                                                           atttgcgcATTTGTGGAACGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTGTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCTCGATTAACTCCATTCCGATCCGGATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGATTTCGGTAGTCCCGgaga                                                                                                                                   <  1:94410/247‑1 (MQ=255)                                                                                                           aTTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCgg                                                                                                                                      >  1:127433/1‑248 (MQ=255)                                                                                                                                         gTTCCAGAATGGTTTCCGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGAGGGGCTTGCGCTGGTAGGTAAAACAGCTGCGGAAACCAAAGATTCAGAGAAATCAATGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGCTTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTACCGGAGCAAAGCCCGTTTGTTGAACTGTTCGTGTTgg                                                                                                     >  2:266364/1‑251 (MQ=255)                                                                                                                                               gAAAGGATTAAGGGGATTCTTTACCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGACACCAAAGATTCAGAGAAATCACTGCAACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGc                                                                                               >  2:166729/1‑251 (MQ=255)                                                                                                                                                    gTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGAGTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGTAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGctgcct                                                                                          <  2:225576/251‑1 (MQ=255)                                                                                                                                                      ttAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGctgcctgc                                                                                        >  1:56321/1‑251 (MQ=255)                                                                                                                                                      ttAAGTGGCTTCTTTGCCGGATTCCAGAGAGAGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGctgcctg                                                                                         <  1:146531/250‑1 (MQ=255)                                                                                                                                                        aaGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCctgct                                                                                       <  2:56321/250‑1 (MQ=255)                                                                                                                                                          gTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGcc                                                                                    <  1:19595/251‑1 (MQ=255)                                                                                                                                                                      ccGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCAGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATc                                                                        >  1:129223/1‑251 (MQ=255)                                                                                                                                                                            tCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACttt                                                                  >  2:248369/1‑251 (MQ=255)                                                                                                                                                                                gATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTgtgg                                                              >  1:15936/1‑251 (MQ=255)                                                                                                                                                                                    gCGGTTTTCGTTTTCGTGGTGGTTGGGCTTGTAGGTACAACTGCTGCGGAAACCAAAGATTCAGAGAAATCGCTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATTTTCTTCGTCTTCGCGCTGATGGTGAGTATGTCCGTGACGCTGTGGAGTTCGGTGGGTTCGGAGAAAAGTCCGTTTGATGAACTGTTCTTGTTGGTGTGGCTGCCGGCGGCCGCGAGCTTGATctaatgtgttgtgct                                                          <  2:263486/251‑15 (MQ=255)                                                                                                                                                                                          ttCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACct                                                     <  1:128591/250‑1 (MQ=255)                                                                                                                                                                                             gCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTgc                                                 >  1:184046/1‑251 (MQ=255)                                                                                                                                                                                                     ggggATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGGTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCtt                                         <  2:40985/251‑1 (MQ=255)                                                                                                                                                                                                             aGCTGGTAGGTACAACCGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCCTTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGACCTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCACCTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAAc                                 >  2:215089/1‑251 (MQ=255)                                                                                                                                                                                                                  gtaggtaCAACAGCTGCGGAAACCAAAGATTAAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATAATTATGTTCTACGTCTTCGCGCTGCTTGTGATTATGTCCGTGACGCCGTGGCGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAgcgg                            >  1:116746/1‑251 (MQ=255)                                                                                                                                                                                                                    aggtaCAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAgcggcg                          >  2:60630/1‑251 (MQ=255)                                                                                                                                                                                                                        acaacaGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGtctt                      >  1:77751/1‑251 (MQ=255)                                                                                                                                                                                                                                        ccAAAGATTCATAGTAATCACTGCCACGCGCGATTAACTCCATTCATATCCGTATCATAATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGGGTTCGGTAATCCCGGAGAAAAGCCCGTTTGTTGACCTGTTCGTGTTGGTAGGGCTGCCTTCTGCCGCAAGCGTGTTCAACTTTGGTGTGCTGACCTCTGCGGCGTCTTCCGCTACCAGCTACGCCTTCTCTACCAGCCTTATg      >  2:45770/1‑251 (MQ=255)                                                                                                                                                                                                                                          aaaGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCt    <  2:158866/251‑1 (MQ=255)                                                                                                                                                                                                                                            aGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCTCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGt  <  1:106151/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       AGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGT  >  NC_000913/4430345‑4430830

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