Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,625,248	C→T	D118D (GAC→GAT)	serB →	phosphoserine phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,625,248	0	C	T	100.0%	160.2 / NA	52	D118D (GAC→GAT)	serB	phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/26);  total (26/26)

GTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTG  >  NC_000913/4625004‑4625487                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                gTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCGGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCAcc                                                                                                                                                                                                                                           <  1:441244/251‑1 (MQ=255)                     ggTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTAt                                                                                                                                                                                                                       >  1:352433/1‑251 (MQ=255)                      gTCGTTGGCTTGATAAGCATCGTCTGACACAATACTATAGCATACTGGGTGCGGCTATGGTGATTGTTGCCGCTTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTCCACGCACCGCACGGGCTACGTGCCTGGCCCACGAAGCGCAACTGTGTGTCGCCCCGCTGGGGAAAAGCTCGCACCTACGCACGCCGTGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTAtt                                                                                                                                                                                                                      <  2:7914/251‑1 (MQ=255)                           gggCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGAATGAATGTATtgatga                                                                                                                                                                                                                 <  1:374444/251‑1 (MQ=255)                            ggCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGaa                                                                                                                                                                                                                <  2:3722/251‑1 (MQ=255)                                  ataaacagcgtgagaaccAATACAAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGAGCGTGGAAGATTATCAGGTGATTCGTATGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGAAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGATGGGGCTGGAAATGGATTCCACCGCAATGCAGATTGAATGGATTGAGGAAAtgtac                                                                                                                                                                                                          >  2:319566/17‑247 (MQ=255)                                      acagcGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGCTTCCACCGCCATCCAGATTGCATGTATTGATGAAATTGCCAAAc                                                                                                                                                                                                      >  2:397743/5‑251 (MQ=255)                                      acagcGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGGGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCCCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGACTGTATTGATGAAATTGCCAAAc                                                                                                                                                                                                      >  1:236390/5‑251 (MQ=255)                                           gTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGcc                                                                                                                                                                                                 >  1:346208/1‑251 (MQ=255)                                                      tACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGGTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAAAGGGCgg                                                                                                                                                                                      >  2:446502/1‑250 (MQ=255)                                                          aaaacaaaaggggTGAGGCGATGTTGATTGGTGCCGCCTGGTGCGTGGAAGATTATCAGGTTATTCGGCTGGCAGGTTCACTCACCGCCTGGGCTACGCGCCGGGCCCATTAAGCGCAACTGTATGGCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGACGTGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGAtg                                                                                                                                                                                  <  2:262884/241‑1 (MQ=255)                                                             gCAAACTGGGTTCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTCCGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGAtggtg                                                                                                                                                                               <  2:290144/251‑1 (MQ=255)                                                               aaaCTGAGTGCGGCGATGGTGAATGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCCCCTGCGCACGCCGGGGTTGCTGGTGCTGGACATGGATTCCACAGCCATCAAGCCTGCATGTATTGATGAAATTGCCCAACTGGCCGGAACGGGCGAGAtggtggg                                                                                                                                                                             >  2:344464/1‑250 (MQ=255)                                                                            cgATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGCTTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGGGATGGTGGCGGAAGTAACCGa                                                                                                                                                                 >  2:254499/1‑250 (MQ=255)                                                                             aaTGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCAGCTTGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGGTTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGGAGTAACCGaa                                                                                                                                                                >  2:269072/2‑250 (MQ=255)                                                                                ggTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACggg                                                                                                                                                            >  1:4363/1‑251 (MQ=255)                                                                                   gATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCgg                                                                                                                                                         >  2:162568/1‑250 (MQ=255)                                                                                    aTTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGAt                                                                                                                                                        >  1:247665/1‑251 (MQ=255)                                                                                    aTTGTTGCCGCCTGGTGCGTGGAAGATTAGCAGGTGTCTCGTCTGGCAGGTTCACTCACCGCACGGTCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGTACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGAt                                                                                                                                                        <  1:433473/251‑1 (MQ=255)                                                                                          gccgccgtgTGCATGGAAGATTATCAGGTGATTTGTCCGGCAGGTTCCCTCTCCGCCCGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGcgcg                                                                                                                                                   <  1:263503/242‑1 (MQ=255)                                                                                                 ggggCGTGGAAGATTATCAGGTGATGCGTCGGGCAGGTTCGATCACCGCACGGGCTACGCGCCTGGTCCATGAAGAGCAACTGGATGTCGTCCCGCGGGGGAAAAGCCCGCAACTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGAGGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTc                                                                                                                                           <  2:250018/248‑1 (MQ=255)                                                                                                 ggTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTgca                                                                                                                                                                                                                            >  2:109289/1‑168 (MQ=255)                                                                                                 ggTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGaa                                                                                                                                                                                                                            <  1:109289/170‑1 (MQ=255)                                                                                                             aTTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGCGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGcc                                                                                                                               >  2:380378/1‑251 (MQ=255)                                                                                                                        aTTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTCGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCCGCCTgcgc                                                                                                                      >  2:391503/1‑249 (MQ=255)                                                                                                                                       tcactcacCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCt                                                                                                     >  1:390732/1‑251 (MQ=255)                                                                                                                                       accctcacCGCCCGGGGCACGCGCCTGGCCCATGAACCGCAATTGGATCTTCCCCCGTTGGGGAAAACCCCGCCCCGGTCCACGCCGGGTCTGCTGGTGAGTGACATGGATTCCGCCGCCATCCAGATGGAATGGATTGAAGAAATGGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGATCGGGGGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAgacgcttgtcgccgc                                                                                                      <  2:221253/247‑15 (MQ=255)                                                                                                                                                  aCGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATtgatga                                                                                                                                                                                                                 >  1:128071/1‑132 (MQ=255)                                                                                                                                                  aCGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATtgatga                                                                                                                                                                                                                 <  2:128071/132‑1 (MQ=255)                                                                                                                                                   cGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTg                                                                                          >  1:198040/1‑250 (MQ=255)                                                                                                                                                                        gAAGCGCAACTGGCTTTCGCCCCCCTGGGGAAAATCCCGCGCCGGGGCCCGCCGGGTTTGCTGGGGAGGGACGTGGATTCCGCCGCCATCCAAATTGAAGGTTTTGAGGAAATTGCCAAACTGGCCGGAAGGGGCGAGATGGGGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGTTTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGt                                                                    <  2:352433/251‑1 (MQ=255)                                                                                                                                                                            cgcACCTGGAGGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGt                                                                <  1:397743/251‑1 (MQ=255)                                                                                                                                                                                  tGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAAt                                                          >  2:439629/1‑251 (MQ=255)                                                                                                                                                                                     aTGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTg                                                       <  1:382108/251‑1 (MQ=255)                                                                                                                                                                                       gcccccccccgggggAAAATCCCTCCCCTGCGCCCGCCGGTTTTGCTGGTGATGGCAAGGGATTCCCCCGCCATCCAGATTGAATGTATGGATGAAATTGCCAACCGGGCCGGAACGGGCGAGATGGGGGCGGAGGTAACCGAACGGGCGATGCGCGGCGACCGCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGcc                                                     <  2:377296/240‑1 (MQ=255)                                                                                                                                                                                           ccccGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTg                                                 >  1:141371/1‑251 (MQ=255)                                                                                                                                                                                              cGCTGGGGAAAACCCCGCACCGGCGCACGCCGGGTTTGTTGGTGATGGACTGGGATTCCACCGCCATCCAGATTGAATGGATGGATGAAATGGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATg                                              <  2:343223/251‑1 (MQ=255)                                                                                                                                                                                                  ggggAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCa                                           <  1:254499/250‑1 (MQ=255)                                                                                                                                                                                                       aaaTCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTa                                     <  1:162568/251‑1 (MQ=255)                                                                                                                                                                                                                aCCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTg                            >  1:133707/1‑251 (MQ=255)                                                                                                                                                                                                                  cTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGAtt                                                                                                                                       >  2:222144/1‑142 (MQ=255)                                                                                                                                                                                                                  cTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGAtt                                                                                                                                       <  1:222144/142‑1 (MQ=255)                                                                                                                                                                                                                     cgcACGCTGGGTTTCCGGGTGATGGACATGGATGCCCCCGCCAGCCAGATGGAATGTATTGATGAATTGGCCAACCAGGCGGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTCCCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCt                       <  2:8361/251‑1 (MQ=255)                                                                                                                                                                                                                            cGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCt                                       <  1:27724/228‑1 (MQ=255)                                                                                                                                                                                                                            cGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCt                                       >  2:27724/1‑228 (MQ=255)                                                                                                                                                                                                                            cGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATGGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTg                <  1:94804/251‑1 (MQ=255)                                                                                                                                                                                                                             gggTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTgg               >  1:378423/1‑251 (MQ=255)                                                                                                                                                                                                                             gggTTTGCTGGAGATGGACATGGATTCCACCGCCATCCAGATTGAATTTCTTGCTGAAATTGACAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGGAAGGCGCTGACGCCAATATTCTGCAACACGGGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTgg               >  2:311292/1‑251 (MQ=255)                                                                                                                                                                                                                               gTTTGCTGGTGAGGGACAGGGCTTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGGGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGaa             <  2:167846/251‑1 (MQ=255)                                                                                                                                                                                                                                      ggTGAGGGACATGGATTCCCCCGCCATCCAGATTAAATGGATTGATGAAATGGCCAAACTGGCCGGAAGGGGCGAGATGGGGGCGGAAGTACCCGAACGGGCGATGCGCGGCAACCTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACgctgg      <  2:240607/251‑1 (MQ=255)                                                                                                                                                                                                                                         gATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACgctgggc    <  1:250808/250‑1 (MQ=255)                                                                                                                                                                                                                                          aTGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACgctgggctg  >  2:165128/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                GTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTG  >  NC_000913/4625004‑4625487

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