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%	195.2 / NA	64	P241S (CCA→TCA)	cycA	D‑serine/alanine/glycine/:H(+)symporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/30);  total (34/30)

AGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCT  >  NC_000913/4430339‑4430828                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     ggtgggAGTTCTGGTGTGCGATGATCAAAATCGGTGCCATTGTGTCGCTGATTGTCGGCGGCCTGGTCAGGGTGGCGAGGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGTGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCagag                                                                                                                                                                                                                                                 <  2:57577/247‑1 (MQ=255)           tgggcttgagatgatCAAAATCGTCGCCATTGTGTCGTTGATTGTCTTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCTCATTTGTGGAATGACTGCGGCTGGTTCCCGAAAGGTTTAAGTGTCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTATGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTg                                                                                                                                                                                                                                        <  2:181047/241‑1 (MQ=255)             gggttGCGATTATCAAAAGCGCCGCCCTTGTGTCGCTGATTTTCGTCGGCCGTGTCATGTTGTCGATGCACTTTCAGTCACCGGCTGGTGTGGAAGCGTCATGCGCGCATTTGTGGAATGACTGCGGCTGGTTCCCTAAAGGTTTAAGTGGCTTCTTTGCCGGAGTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGcc                                                                                                                                                                                                                                      <  2:392397/247‑1 (MQ=255)             ggTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGcc                                                                                                                                                                                                                                      >  1:461917/1‑250 (MQ=255)                       gatCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCAcgcgcgctt                                                                                                                                                                                                                            >  1:416782/1‑247 (MQ=255)                        atCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGGGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGCTTGAGCTGGTAGGTACAAAAGCTGCGGAAACCAACGATTCAGAGAAATAACTGCCAcgcgctctta                                                                                                                                                                                                                           >  2:130926/1‑245 (MQ=255)                          cAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGCTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACACCAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGCTTAAcc                                                                                                                                                                                                                        >  2:50358/1‑250 (MQ=255)                           aaaaTCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGGAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGTTTCAGAGAAATCACTGCCACGCGCGATTAACTc                                                                                                                                                                                                                       >  2:296368/1‑251 (MQ=255)                            aaaTCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCCTGGTGGCGATGCACTTTCAGTCACCGACTGGAGTGGAAGCGTCATTCGCGCATTTGTGGAATGCCGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATCGCGGTTTTCGCTTTCGTGGGGTTTGCGCTGGTAGGTACAACAGCTGCGGAAACCAAAGCTTCAGGGAAATCACTGCCCCGCGCGATTAACTcc                                                                                                                                                                                                                      >  2:284967/1‑251 (MQ=255)                                      cATTGTGGCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCCCTTTCAGTCACCGACTGGTTTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATcc                                                                                                                                                                                                            <  2:54975/251‑1 (MQ=255)                                       aTTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCg                                                                                                                                                                                                           >  2:290777/1‑251 (MQ=255)                                            gtCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCCCTGCCACGCGCGATTAACTCCATTCCGATCCGTATcc                                                                                                                                                                                                      >  2:321232/1‑250 (MQ=255)                                             tCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCAt                                                                                                                                                                                                     >  1:245606/1‑251 (MQ=255)                                                    ttGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTc                                                                                                                                                                                              >  1:366086/1‑251 (MQ=255)                                                     tGTCGTCGGCCTGGTCATGGTGGCGATGCCCTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTc                                                                                                                                                                                              <  1:296368/250‑1 (MQ=255)                                                     tGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCt                                                                                                                                                                                             >  1:205357/1‑251 (MQ=255)                                                      ggcgtcgGGCTGGTCATGGTGGCGATTCCCTTTCAATCACCGACTGGTGGTGAAGCGGTAATCGCGCATTTGTGGAATGACGGGGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGGTTTCGCTTTCGGGGGGATTGAGCTGGTAAGGACAACAGCTTCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTa                                                                                                                                                                                            <  1:388433/249‑1 (MQ=255)                                                        cgtcgGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACg                                                                                                                                                                                          <  2:245606/251‑1 (MQ=255)                                                            ggcatggtaatgatgtcGATGCATCTTCAGTCACCGACTTGTGTGGAAGCGTCATTCGCGCATTTGTTGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTGTCTTGGGGACTGAGCTGGTATGTCCAACAGCTGCGTAAACCAAAGATTCAGTGAAATCACTGCCACGCTCTATTAACTCCATTCCGATCCGTATCATTCTGTTCTACGTCtt                                                                                                                                                                                      <  1:237253/235‑1 (MQ=255)                                                            ggCCTGGGGATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGCAAGGTTTCAGGGGCTTCTTTGCCGGATTCCAGATAGCGGGTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGGAATCACTGCCACGTGCGATTAACTCCATTCCGATCCCTATCATTATGTTCTAAGTCtt                                                                                                                                                                                      >  2:88557/1‑251 (MQ=255)                                                              ccTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATGCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAATGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCTTGGGGAATGAGCTGGTAGGTACACCAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCTATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTcg                                                                                                                                                                                    <  1:284967/251‑1 (MQ=255)                                                               cTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTcgc                                                                                                                                                                                   >  1:383114/1‑251 (MQ=255)                                                                   tCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTg                                                                                                                                                                               >  2:476089/1‑251 (MQ=255)                                                                        gtggCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATtgt                                                                                                                                                                          >  2:195278/1‑251 (MQ=255)                                                                          ggCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGCAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCCTTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATtgtgg                                                                                                                                                                        >  2:240250/1‑250 (MQ=255)                                                                               tGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATAACTGCCACGCGCGATTAACTCCATTCCGGTCCGTATCATTATGTTCTACGTCTTCGCGCTGATtgtggttttg                                                                                                                                                                   >  2:26003/1‑245 (MQ=255)                                                                                                aCTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAg                                                                                                                                                  <  2:247862/251‑1 (MQ=255)                                                                                                aCTGGTGTGGAAGCGTCATTCGCGACTTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGGTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGAGATTATGTCCGTGACGCCGTgg                                                                                                                                                    >  1:341424/1‑249 (MQ=255)                                                                                                       gggAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTa                                                                                                                                           <  1:130926/250‑1 (MQ=255)                                                                                                              gTCATTCGAGCATTTGAGGAATGACGGCGGCTTGTTCACGCACGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGAATGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCCGAGAAATCACTGCCACGCGCGCTTAACTCAATTCCGATCCGTTTCATTATTTTCTACGTCTTAGCGCTGATTGTGTTTATGTCCGTGAATCCGTGTTGTTCTGTAGTccgac                                                                                                                                    >  2:375939/1‑248 (MQ=255)                                                                                                                cATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACCACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCAAGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTCTGTCCGTAACGCCGTGGAGTTCGGTAGTCCCgg                                                                                                                                    >  2:378171/1‑249 (MQ=255)                                                                                                                    cgcgCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGaaaa                                                                                                                              <  2:461917/251‑1 (MQ=255)                                                                                                                     gcgcATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAg                                                                                                                             <  1:240250/251‑1 (MQ=255)                                                                                                                                 aaTGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGa                                                                                                                 <  1:111694/251‑1 (MQ=255)                                                                                                                                  agtaCGGCGTCTGGTTCCCGAAAGTTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTGCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGAGTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGAGTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGaa                                                                                                                <  2:504944/248‑1 (MQ=255)                                                                                                                                      cggcggCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGt                                                                                                            >  1:341308/1‑251 (MQ=255)                                                                                                                                      cggcggCTGGTTCCCGAAAGGTTTAAGTGACTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGCAACCAACGATTCAGAGGAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTCCGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGCTCGGTAGTCCGGGCGAAAAGCCCGTTTGTTGGACTGt                                                                                                            >  2:83153/1‑251 (MQ=255)                                                                                                                                       ggcggcTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGtt                                                                                                           <  1:124281/251‑1 (MQ=255)                                                                                                                                        gcggatgtgttctgctaGGTGTAAGTGGTATCTTTGCCGGATTCCAGATAGCTGGTTTCGACGTCGGGGGGAGTGAGCTGGTAGGTTCAGCAGCTGCGTAAACCAAAGATTCAGAGAAATCACTGCCTCGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCGACGTCATCGCGCTGATTGTGATTATTTTCGTGACGCCGTTGAGTTCTGTAGTCCCGGAGAAAAGCCCGTTTGTTGCACTGTTc                                                                                                          <  2:287102/235‑1 (MQ=255)                                                                                                                                          ggcTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCgt                                                                                                        <  1:211952/251‑1 (MQ=255)                                                                                                                                               gTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGCGCTGGTAGGTACCACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTCTGTTCTACGTCTTCGCGCTGCTTGTGATTATGTCCGTGACGCCGTGGCGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGACCTGTTCGTGTTgg                                                                                                   >  2:223310/1‑251 (MQ=255)                                                                                                                                                   ccGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAgg                                                                                               >  1:167005/1‑251 (MQ=255)                                                                                                                                                    cGAAAGGTTTAAGTGACTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAggg                                                                                              >  2:99681/1‑251 (MQ=255)                                                                                                                                                         ggTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGATCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGctgcc                                                                                         <  1:101188/251‑1 (MQ=255)                                                                                                                                                              aaGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCctgct                                                                                     <  1:440970/250‑1 (MQ=255)                                                                                                                                                                gTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGcc                                                                                  >  1:81049/1‑251 (MQ=255)                                                                                                                                                                         ttGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGt                                                                          <  2:334101/250‑1 (MQ=255)                                                                                                                                                                                aTTCCAGATAGCGGTTGTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACt                                                                  <  2:366086/251‑1 (MQ=255)                                                                                                                                                                                     aGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGGAGTCACTGCCACGCGCGATTAAATCCATTCCGATCAGTATAATTATGTTCTGCTTCTTCGCGCTGATTGTGATTATGTCCGTGGCGACGTGGAGTTCGGGAGTCCCAGAGAAAAGCACGTTTGTTGCACTGTTCGTGTTAGTGGGACTGCCTGATGCCGGAAGCGTGATCAACTTtg                                                               >  1:417547/1‑249 (MQ=255)                                                                                                                                                                                     aGATAGAGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTAGAACAGCTGCGGAAACCAAAGATTCAGAGACATCACTGCCACGCGCGATTAACTCCATTCCGATCAGTATCATTATGTTCTACTTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTCGTCCCGGAGAAAAGCCCGTTTGTTGACCTGTTCGTGTTGGTAGGGCTGCCTGCTGCATCAAGCGTGATCAACTTtgtg                                                             >  1:314563/1‑251 (MQ=255)                                                                                                                                                                                      gATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTgtgg                                                            >  1:234796/1‑251 (MQ=255)                                                                                                                                                                                        tAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTgtggtg                                                          >  1:498865/1‑251 (MQ=255)                                                                                                                                                                                         aGCGTTTTTCGCTGTCGTGGGGATTGAGCTGGTAGGTTCAACAGTTGTGGAAACCAAAGATTCAGAGAAATCACTGTCCCGCGCGATTAACTCCATTCCGCTCCGTATCATTATGTTCTACGTCGTCGCGCTGATTGTGATTATGTCCGGGACGCCGTGGAGTTCGGTAGTTCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTgtggtg                                                          <  1:470633/250‑1 (MQ=255)                                                                                                                                                                                         aGCGGTTTTCGCTTTCGGGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTgtggtg                                                          <  2:383114/250‑1 (MQ=255)                                                                                                                                                                                                ttCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACctc                                                  >  1:486213/1‑251 (MQ=255)                                                                                                                                                                                                      ttCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTgcg                                              <  1:26003/249‑1 (MQ=255)                                                                                                                                                                                                           ggggATTGAGCTGGTAGGTACAACAGCGGCGGAAACCAAAGATTCAGAGAAATCACTGCCTCGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCt                                        <  2:269537/250‑1 (MQ=255)                                                                                                                                                                                                              gATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTACATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCg                                    >  1:35624/1‑251 (MQ=255)                                                                                                                                                                                                                       ggtaggtaCAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAgcg                           <  1:196466/251‑1 (MQ=255)                                                                                                                                                                                                                       ggtaggtaCAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAgcg                           <  1:132160/251‑1 (MQ=255)                                                                                                                                                                                                                              acaacaGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGtct                     <  1:4915/250‑1 (MQ=255)                                                                                                                                                                                                                               caacaGCTGCGGAAACCAAAGATTACGAGAAATCACAGCCACGCGCGGTTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGACAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACCACGGCGtcttc                   >  1:292028/1‑251 (MQ=255)                                                                                                                                                                                                                                                aaaGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCt  >  1:289451/1‑251 (MQ=255)                                                                                                                                                                                                                                                aaaGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTGCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCATCCGTATGCt  <  1:440073/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     AGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCT  >  NC_000913/4430339‑4430828

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: ‑