Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,261	A→T	I362L (ATA→TTA)	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,261	0	A	T	100.0%	129.0 / NA	41	I362L (ATA→TTA)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (22/19);  total (22/19)

GTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGAC  >  NC_000913/4639015‑4639506                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        gTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAAt                                                                                                                                                                                                                                                     <  1:93008/251‑1 (MQ=255)         cacgagcccAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGta                                                                                                                                                                                                                                             <  2:638004/244‑1 (MQ=255)                 cAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTACACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTACTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCGTTCTGATTGGAGCGTTAATGAACGGAATTTATTTg                                                                                                                                                                                                                                     >  1:379025/1‑251 (MQ=255)                         cAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGt                                                                                                                                                                                                                             >  1:12310/1‑251 (MQ=255)                                     acggACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAACCGCTAACCGAGCGACGTTTAAACCCAATGTAATATTTGCTGGTGGGGCTTTCATTGGGTATGTTTTATTTGATCTTGCTGGCGCTTTCTGAACATATCGGTTTTTCCGAGGCAGGGATAATCGACAGTCTTATTGGGGCGTTAATGCACGGGatttatttacagggggtttttgaaagcgg                                                                                                                                                                                                                 >  2:262603/4‑230 (MQ=255)                                       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)                                                                                                                                                                  ggTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGtattta                                                                                                                                                                                                                                         >  1:508298/1‑102 (MQ=255)                                                                                                                                                                  ggTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGtattta                                                                                                                                                                                                                                         <  2:508298/102‑1 (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)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        GTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGC‑AACA‑GCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGAC  >  NC_000913/4639015‑4639506

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