Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,303,755	C→T	V119I (GTT→ATT)	mdtN ←	putative multidrug efflux pump membrane fusion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,303,755	0	C	T	98.5%	189.0 / NA	67	V119I (GTT→ATT)	mdtN	putative multidrug efflux pump membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (36/29);  minor base G (0/1);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAATATCCGCTTCGACCGCCGCACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTT  >  NC_000913/4303511‑4303996                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  gcaaaacccGCTGCGACCGCCGCACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTc                                                                                                                                                                                                                                             <  2:174112/244‑1 (MQ=255)      atCCGCTTCGACCGCCGCACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTACGATCCGGTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGACAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACAGTGGCATTAATAGAGTCGGCAc                                                                                                                                                                                                                                        >  1:359955/1‑251 (MQ=255)      atCCGCTTCGACCGCCGCACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCAc                                                                                                                                                                                                                                        >  1:557267/1‑251 (MQ=255)             tcgcccgccgCACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAAc                                                                                                                                                                                                                                   <  2:192378/245‑1 (MQ=255)               gACCGCCGCACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTACCAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTTCCGGGGCATTAATTGAGTCGGCACCCAACTGtt                                                                                                                                                                                                                               >  2:192075/1‑251 (MQ=255)                 ccgccgCACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTgc                                                                                                                                                                                                                             >  1:521364/1‑251 (MQ=255)                    cagcACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTTTGGCCTGTTTCGCGGCGGCACGGGGCTTTTCTACCGTGGCATTAATCGATTCGGCAccctcctggtgcg                                                                                                                                                                                                                            >  2:440544/3‑238 (MQ=255)                     cgcACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATCCGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGACGTTCTTGCACGGTCAACCTCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATATGTGGCCTGTTTCGCGGCGGACCGGGCTTTTTCTACCGTGGCATTAATCGTGTCGGCCCCAAACTGTTGCGCGt                                                                                                                                                                                                                         >  2:589266/1‑251 (MQ=255)                         cGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGTCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAac                                                                                                                                                                                                                     >  2:456113/1‑251 (MQ=255)                          gCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAacg                                                                                                                                                                                                                    >  1:7675/1‑251 (MQ=255)                          gCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATGAATCGAGTCGGCACCAAACTGTTGCGCGTCAacg                                                                                                                                                                                                                    <  2:115865/251‑1 (MQ=255)                           cTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGGGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCAAGGGCTTTTTATACCGTGGCATTAATCGAGGCGGCACCAAACTGTTGCGCGTCAacgc                                                                                                                                                                                                                   >  2:323395/1‑251 (MQ=255)                              ggCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAacgctac                                                                                                                                                                                                                <  1:250223/251‑1 (MQ=255)                              ggCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGCTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAacgctac                                                                                                                                                                                                                >  1:204576/1‑251 (MQ=255)                                cAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAacgctacgc                                                                                                                                                                                                              >  2:190058/1‑251 (MQ=255)                                cAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAAGACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGGCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAacgctacgc                                                                                                                                                                                                              <  2:119319/251‑1 (MQ=255)                                          aTCCACGCCTCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCa                                                                                                                                                                                                    <  1:589266/251‑1 (MQ=255)                                                         ggCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCCCGGGCTTTTTCTACCGTGGCATTAATCGCGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGGGCCTAATTTGCTTATccc                                                                                                                                                                                     >  2:254350/1‑250 (MQ=255)                                                         ggCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATcc                                                                                                                                                                                      <  1:254350/250‑1 (MQ=255)                                                            gctggcgcctgcatgagcgcgTAAAAATACGCAATTAAAATCCGATTCTGCGACGCACTCCGCCGTCCTTGCACGGTCACCATCTTCCGCTGAGCCAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTCATGTATCCCTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCCAGGCATTACTCGAGTCGGCACCAAACAGTTGCGCGTAAACGCTCCGCTGGGTGGACCTAATTTGCTTATCCAgcg                                                                                                                                                                                  >  2:369084/21‑251 (MQ=255)                                                                     tGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCgaggg                                                                                                                                                                         >  2:165770/1‑251 (MQ=255)                                                                       aCTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCgagggag                                                                                                                                                                       <  1:191755/251‑1 (MQ=255)                                                                           cgcCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCtt                                                                                                                                                                   >  1:232025/1‑251 (MQ=255)                                                                                 tAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAACCTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGCGGCTTCAGCtt                                                                                                                                                             >  2:522231/1‑251 (MQ=255)                                                                                  aaCAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTGCGCGGCGGCACGGGCTTTTTCGTCCGTGGCATTAATCGAGTCGGCTCCAAACTGTTGCGCGTCAGCGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCttt                                                                                                                                                            <  2:272453/251‑1 (MQ=255)                                                                                                  gATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCCCGGTCAACATCTTCCGCTGAGACCAACCCTTCTTTCTGTAATAGCTCGGTGCGGCGTAATGTAACTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTAACGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGATTATCCCGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCg                                                                                                                                            >  2:239277/1‑251 (MQ=255)                                                                                                  gATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCg                                                                                                                                            >  1:94031/1‑251 (MQ=255)                                                                                                               ggcgCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTCTCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCCACGCTACGCTGGGGGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCGGATTGGCTTCGTACGGGCGCgggg                                                                                                                               >  2:577473/1‑250 (MQ=255)                                                                                                               ggcgCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGt                                                                                                                               >  1:43325/1‑251 (MQ=255)                                                                                                                    gcTGCGCCGGTCGTGCACGGTCAACCTCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGTGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGACGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATg                                                                                                                          <  2:415920/251‑1 (MQ=255)                                                                                                                      tGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCg                                                                                                                        >  1:371705/1‑251 (MQ=255)                                                                                                                      tGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGGCGATGCg                                                                                                                        >  2:287610/1‑251 (MQ=255)                                                                                                                              tCTTGCCCGGTCAACATCTTCCGCTGAGACAAACCCTTCTTGCAGTAATTGCGCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCTAGTGAGGATTCAGCTTTCGCCAGATTGGCTTTGTACTGGCGCTGGTCTATGCGCAACAGc                                                                                                                 <  2:352100/250‑1 (MQ=255)                                                                                                                                ttGCACGGTCAACATATTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTATACCGTGGCATTAATCGAGTCGGAACCAACCTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGAGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGACGATGCGGAACAGCaa                                                                                                               >  1:561293/1‑250 (MQ=255)                                                                                                                                 tGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGCTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAAt                                                                                                             >  1:469153/1‑251 (MQ=255)                                                                                                                                  gCACGGGCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATc                                                                                                            <  2:251323/251‑1 (MQ=255)                                                                                                                                   cACGGTAAACATCTTCCGCTGCGAAAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCCAGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGGGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGGCGATGCGGAACAGCAAATCg                                                                                                           >  2:561159/1‑251 (MQ=255)                                                                                                                                           aCATTTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAAGGTATCTGTGGCCTGTTTCGCTGCCGCACGGGCTTTTCCTACCGTGGCATTAATCGAGTCGGCGCCAAACTGTTGCGCGTCAACGCTGCGCGGTGTGAGCATAATTTGCTTATCCAGCGCCGCGATGGATGCTTCAGCTGTCGCCAGATTGGTTTCGTACGGGCTCGTGTCGATGCCTAGCAGTAAATCGCCGTGttt                                                                                                   <  2:277533/251‑1 (MQ=255)                                                                                                                                                  ccGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGcc                                                                                            >  1:518946/1‑251 (MQ=255)                                                                                                                                                   cGCTGAGACAAATCCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGTCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCt                                                                                           <  1:24404/251‑1 (MQ=255)                                                                                                                                                    gCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTg                                                                                          >  1:105751/1‑251 (MQ=255)                                                                                                                                                    gCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTGTTCTACCGTGGCATTAAGCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTg                                                                                          <  1:527277/251‑1 (MQ=255)                                                                                                                                                        agaCAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTg                                                                                      <  2:94031/251‑1 (MQ=255)                                                                                                                                                            aaacccttcttcaattaagtgcgcGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTGCTTCCGTGGTATTAATCGAGTCTGCCCCAAGCTGTTGCGCGTCAACGCTACGCTGGGTGAGTATAATTTTCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGACgg                                                                                  <  2:462034/228‑1 (MQ=255)                                                                                                                                                            aaaaCCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCgg                                                                                  <  1:165770/251‑1 (MQ=255)                                                                                                                                                                      tgcAGTAATGGCTCGGTGCGGCTTAATGTAGCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGcc                                                                          <  2:317713/247‑1 (MQ=255)                                                                                                                                                                       tCAGTAATGGCTCCGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCGTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTGCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGTCTGGTTGTCGGTGACCGCCAg                                                                        <  2:564758/250‑1 (MQ=255)                                                                                                                                                                            aaTGGCTGGGTGGGGCTTAATGTACCTGTGGCCTGTTTCGCGGCGGCACGGGTTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTTCTTATCCAGCGCCGCGAGGGAGTCTTCAGCTTTCGCCAGATTGGCTTCGTTCGGGCGCGGGTCGATGCGTAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTa                                                                   <  1:521890/250‑1 (MQ=255)                                                                                                                                                                               ggCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAAt                                                               <  1:550494/251‑1 (MQ=255)                                                                                                                                                                                 cccGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCTCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGc                                                             <  1:71109/249‑1 (MQ=255)                                                                                                                                                                                  ccGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCg                                                            <  2:140204/250‑1 (MQ=255)                                                                                                                                                                                   cGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCg                                                            >  1:140204/1‑250 (MQ=255)                                                                                                                                                                                    ggTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGc                                                          >  1:29889/1‑251 (MQ=255)                                                                                                                                                                                     gTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCTTGTTTGTCGGTGACCGCCAGTTCTACAATGCGGcc                                                         <  1:199871/251‑1 (MQ=255)                                                                                                                                                                                           cgTAATGTATCTGTGGCCTGCTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGAc                                                   >  1:407717/1‑251 (MQ=255)                                                                                                                                                                                              aaTGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTc                                                >  1:295112/1‑251 (MQ=255)                                                                                                                                                                                              aaTGTATCTGTGGCCTGTGTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGACCAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACtt                                                 <  2:371705/250‑1 (MQ=255)                                                                                                                                                                                                tgcaTCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCAAACGTAGCATTAATCGAGTCGGCAACAAACTGTTGCGCGTAAACGCTACGCTGGGTGAGCCTAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGACTTCGTACGGGCGCGGGTCGATGCGGAACAGCACATCGCCCTGTTCGACTGACTGGTTGTCGGCGTCCGACAGTTCTACAATGCGTCCGCTGCCTTCCg                                              >  2:44533/4‑251 (MQ=255)                                                                                                                                                                                                     cTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTAcc                                         >  2:164837/1‑251 (MQ=255)                                                                                                                                                                                                                  cgcgGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTg                                                                                                                                                                                             >  1:108782/1‑90 (MQ=255)                                                                                                                                                                                                                  cgcgGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGGCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTg                                                                                                                                                                                             <  2:108782/90‑1 (MQ=255)                                                                                                                                                                                                                    ggtcgggatgggCTTTTCCTACCTGGGCATGAAAGGAGGCGGCACCAGATGGTGGCGCGTAAACGCGACTCTGGGTGATCATAATTGGTTGATCCAGCGCCGCGAGGGAGGCTTAAGCTTTCGCCAGAGTGGCTTGGTACGGGCGCGGGTCGATGGGGAACAGCATGTCGCCCGGTTTGACTGCCTGTTGGGCGGTGTCCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATct                          <  2:580674/242‑1 (MQ=255)                                                                                                                                                                                                                     ggcggcACGGGCCTTTTCTACCGTGGAATGAATGGGGTCGGCACCAAACGGTGGCGCGTCACCGGTACGCGGGGTGAGCATAATTGTCTGATCCCGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCTCGGGTCGATGCGGAACAGCAAATCGCCCGGTTTGACGGCCGGGTGGTCGTTGACCGCCAGTGCTTCAATGCGTCCGCTGACTTCCGGTGCCACATCAATTGGATctg                         <  2:241041/251‑1 (MQ=255)                                                                                                                                                                                                                             gggCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTACCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGAcgcg                 <  1:107358/251‑1 (MQ=255)                                                                                                                                                                                                                                     tctcccGTGGCATTAATCGAGTCGTCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGGCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGt          <  1:334937/246‑1 (MQ=255)                                                                                                                                                                                                                                      cTACCGTGGCCTTAATCGAGTCGGCAACAAACTGTTGCGCGTCAACGCTACGCTGGGTGCGCATAATTTGCTTATCCCGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGCTGCGGACCAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTCAGCGTcc        >  2:237255/1‑250 (MQ=255)                                                                                                                                                                                                                                            tGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGtt  >  1:385275/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  GCAATATCCGCTTCGACCGCCGCACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTT  >  NC_000913/4303511‑4303996

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