Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,299	T→C	G374G (GGT→GGC)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,639,299	0	T	C	100.0%	98.5 / NA	43	G374G (GGT→GGC)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/23);  total (20/23)

TGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAA  >  NC_000913/4639053‑4639543                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       ggACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTggg                                                                                                                                                                                                                                                    >  1:598895/2‑250 (MQ=255)     cGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGGGC‑a                                                                                                                                                                                                                                                >  1:629753/1‑250 (MQ=255)             aaaGTACGCTATGCGGCGGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGTTCCCCGCGCAACGTTTACACCCAATGCAAGTTTTGCTGGTTGGGCTTTCATTGGTGATGTTTTATTTGTTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCAtgt                                                                                                                                                                                                                                      <  2:614657/250‑1 (MQ=255)                 tACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTTGCGC‑AACA‑GCATgttgtt                                                                                                                                                                                                                                  >  1:171693/1‑251 (MQ=255)                  aCGCCAGTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGttt                                                                                                                                                                                                                                 <  1:236609/251‑1 (MQ=255)                                  tGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTTTTGGCAGGCTGGCGC‑AACA‑GCATTTTGTTTACACTctggctgtggt                                                                                                                                                                                                                 >  1:176875/1‑241 (MQ=255)                                      cTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCTCCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGTGCTTTCAGTTGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGttgttg                                                                                                                                                                                                               <  2:641437/249‑1 (MQ=255)                                         aCTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGGGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGCTTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGCAAGGCCGGCGC‑G‑CACGCAAGTTGTTTTCCCTCGCGCTGTTGTTGCTggg                                                                                                                                                                                                          >  1:640324/1‑250 (MQ=255)                                            tttATGGCGTGCTTTTTGTTTGAAACGCTGTCCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATTTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATg                                                                                                                                                                                                        <  2:400170/250‑1 (MQ=255)                                                     ttCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGGAAGGCTGGCGC‑ACCA‑GCCTGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGGGTGAtttg                                                                                                                                                                                              >  1:218694/1‑248 (MQ=255)                                                           gTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGGAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACt                                                                                                                                                                                        >  1:344289/1‑251 (MQ=255)                                                                   aaCGCTAACCGCGCAACGTTTACACCCAATGCAATATTTGCTAGTGGAGCTTTCATAGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGTTTAGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTTAAATGCTGGCGG‑AACA‑GCATGTTTTTTTCCCCAGAGCTGTTGTTGATTGTGGGGGTTGTGTGGGGACTGCTCaaac                                                                                                                                                                                >  2:524876/1‑249 (MQ=255)                                                                    aCGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGGTGTGGGGGCTGCTCAActc                                                                                                                                                                               >  1:428324/1‑251 (MQ=255)                                                                    aCGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGAtgtg                                                                                                                                                                                              <  1:140393/236‑1 (MQ=255)                                                                    aCGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGAATTTAAAGGCTGGCGC‑AACA‑GACTGTTGTTTACCCTCGCGCTGTTGTTGCTGGGTGGGGTGAtgtg                                                                                                                                                                                              >  2:140393/1‑236 (MQ=255)                                                                      gCTCACCGCGCAACGTTTCCTCCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTg                                                                                                                                                                             <  2:12310/251‑1 (MQ=255)                                                                              cgcAATGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCCGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACC‑GGATATTGTTTTCCCTCGCGCTGTTGTTGATGGATGGTGTGATGTGGGGCCTGCTCAACTCTGCCGAtcgc                                                                                                                                                                     >  1:111975/1‑248 (MQ=255)                                                                                                      atTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGtgtg                                                                                                                                             >  1:595851/1‑251 (MQ=255)                                                                                                      atTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGGGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTggggaaccgtgtg                                                                                                                                             >  1:358760/1‑241 (MQ=255)                                                                                                                                aTGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCatgat                                                                                                                   <  1:370782/251‑1 (MQ=255)                                                                                                                                  gTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTTCCCTCGCGCTGTTGTTGCTGGGTGGTGTGATGTGGGGGCTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGGACCCGTGTGCTGGTGGGGGCGCTGGCCGGCATGgtgt                                                                                                                 >  1:274055/1‑251 (MQ=255)                                                                                                                                     gtATTTGCTCTTGCTGTCGCTTTCTGAACATATCGGTTTGACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGttt                                                                                                               <  2:595851/249‑1 (MQ=255)                                                                                                                                                              tagtcaatcGGTTTTTCCGTGGTATGGTTGATCGCCATGCTGATTGGGGCGTTAACTATCGGTTTTTATTTGCAGTCGGTATTGAAACGCTG‑TGCCAACA‑GCATGTTGGTTTCCCTCGCTGTGTTGTTGCTGGATGGTGTTATGTGGTGACTGCTCATCTCGGCCGATAGCGTGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGTCATGATGTTTGTGACTCGGAATATCGACTGGTATg                                                                                      <  2:203731/242‑1 (MQ=255)                                                                                                                                                               aaCATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGGTGGTGTGATGTGGGGGCTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGGCTGGTATGCg                                                                                    >  1:127403/1‑251 (MQ=255)                                                                                                                                                                        gTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTg                                                                           <  1:474330/251‑1 (MQ=255)                                                                                                                                                                             aCCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGGTGGGGTGATGTGGGGGCTGCTCAACTCTGCCGCTAGCGCGCTGTTGTTTGGGACCGGTGTGCTGGTGGTGGGGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCggg                                                                      >  1:60846/1‑249 (MQ=255)                                                                                                                                                                             aCCGTGGCATGGATAATCGCCAGTCGGATTGGTGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCg                                                                        <  2:60846/249‑1 (MQ=255)                                                                                                                                                                                gTGGCATGGATAATCGGCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCATGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAAt                                                                   <  2:274055/251‑1 (MQ=255)                                                                                                                                                                                     attgATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGAGTGGTGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGGTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGGATGCGGTTTCACTGCCGAAAATGaaa                                                               >  1:641644/4‑250 (MQ=255)                                                                                                                                                                                     aTGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTTCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGGACCAGTGGGCTGGAGGTGGCGCTGGCCGGCCTGATGTTTGTGACTCGTAATATCGGCTGGGATGCGTTTTTACTGCCGAAAATGaaa                                                               >  1:87648/1‑250 (MQ=255)                                                                                                                                                                                       ggATAATCTCAAGTCTGATTGGGGCGTTAATGAACGGTATGTATTTGCAGGCGGTATTGAAAGGCTGGCGT‑AATA‑GCATGTTGTTTTCCCTCGCGCTGTTGTTGCTGGAGGGTGTGATGTGGGGACTGCTGAAATCTTCCGATAGCGCTCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCGGGCCGGCATGATGTTTGTGACTCGTAATAGCGACTGGTATGCGTTTTCACTGCTTAAAATGAAAGc                                                             <  2:509598/250‑1 (MQ=255)                                                                                                                                                                                           aaTCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTGCACTGCCGAAAATGAAAGCCAGt                                                         <  1:246533/250‑1 (MQ=255)                                                                                                                                                                                              cGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTaaa                                                      <  1:587575/250‑1 (MQ=255)                                                                                                                                                                                                            gggCGTTAATGAACGGTAGTTATTTGCAGGCGGTAGTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGGGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGAc                                       <  1:347379/251‑1 (MQ=255)                                                                                                                                                                                                                  tAATTAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGAtca                                  <  1:495889/250‑3 (MQ=255)                                                                                                                                                                                                                     ggAACGGTATTTATTTGCTGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGTTGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGAtcaatta                              <  2:344289/250‑7 (MQ=255)                                                                                                                                                                                                                             atttatttGCAGGCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCt                       <  1:236991/250‑1 (MQ=255)                                                                                                                                                                                                                                    ttcaggcgttATTTAAAGGCTGGCGC‑AACC‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGTTGGTGTGATGTGGGGACTTCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTTGTGCTGGCCGGCAAGATGTTTTTGACTCGTAATATCGACTGGTATGCGTTTTCTCTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGAGCAATTACGTATCTGGAAATa                <  2:528808/241‑1 (MQ=255)                                                                                                                                                                                                                                       aggCGGTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGt            <  1:291061/249‑1 (MQ=255)                                                                                                                                                                                                                                        ggctgTATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCGGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGAGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGtt           <  2:87648/247‑1 (MQ=255)                                                                                                                                                                                                                                             tATTGAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAAAGGCCGATCAATTACGTATCTGGAAATAAGGTTGaaa       >  2:507914/1‑250 (MQ=255)                                                                                                                                                                                                                                                 gAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCTAAAATGAAAGCCAGTAAAGAAGTTAAAATGGACGATCAATTACGTATCTGTAAATAAGGTTGAAAAATa   >  2:356523/1‑250 (MQ=255)                                                                                                                                                                                                                                                 gAAAGGCTGGCGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCGGTATGGTGTGATGTGGGGACTTCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGGTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACTTATCTGGACATAAGGTTGAAAAATaa  <  2:86822/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       TGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAA  >  NC_000913/4639053‑4639543

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