Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,811,848	G→A	S141N (AGC→AAC)	emrA →	multidrug efflux pump membrane fusion protein EmrA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,811,848	0	G	A	97.6%	108.7 / NA	41	S141N (AGC→AAC)	emrA	multidrug efflux pump membrane fusion protein EmrA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/15);  minor base C (0/1);  total (25/16)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATCAAATTCAAATTATGTCTCAGGTGTCTGGCAGCGTGACGAAAGTCTGGGCCGATAACACCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAGCGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATG  >  NC_000913/2811603‑2812092                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     gAATCAAATTCAAATTATGTCTCAGGTGTCTGGCAGCGTGACGAAAGTCTGGGCCGATAACACCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACt                                                                                                                                                                                                                                                 <  1:329951/251‑1 (MQ=255)     caaattcaaattATGTCTCAGGTTTCTGGCAGCGTGACGAAAGTCTGGGCCGATAACACCGATTTTGTAAAAGAATGCGGCTGGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAATCCAAAACGTCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCTGAAAATCGCCCTCGCGAAAGCACAAAACGACTACa                                                                                                                                                                                                                                              <  1:42273/250‑1 (MQ=255)        attcaaattATGTCTCAGGTGTCTGGCAGCGTGACGAAAGTCTGGGCCGATAACACCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAAccc                                                                                                                                                                                                                                          >  2:104782/1‑250 (MQ=255)                tATGTCTCAGGTGTCTGGCAGCGTGACGAAAGTCTGGGCCGATAACACCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGGGTGc                                                                                                                                                                                                                                  >  2:10638/1‑251 (MQ=255)                  tGTCTCAGGTGTCTGGCAGCTTGACGAAAGTCTGGGCCGATAACCCCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGcc                                                                                                                                                                                                                                 <  1:322675/250‑1 (MQ=255)                         ggTGTCTGGCAGCGTGACGAAAGTCTGGGCCGATAACACCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGc                                                                                                                                                                                                                          >  1:19740/1‑250 (MQ=255)                              cgggaaGCGTGTCGAAAGTCTGGGCCGATAACCCCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCAAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGAACAAAACGACTACACCCGCCGTGTGCCGCTGGGCAATGc                                                                                                                                                                                                                     <  1:93458/245‑1 (MQ=255)                                   aGCGTGACGAAAGACTGGGCCGATAACACCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCCGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCCGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAACACAAAAAGACTACAACAGCCGTGTGCCGCTGGGCAATGCCAAcc                                                                                                                                                                                                                >  2:297403/1‑250 (MQ=255)                                          cGAAAGTCTGGGCCGATAACACCGATTTTGTAAAAGAAAGCGCCGTGCTGGCCACTCTCGACCCGACAGCTGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAAAAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGt                                                                                                                                                                                                        >  1:284886/1‑251 (MQ=255)                                            aaaGTCTGGGCCGATAACACCTATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCGGATGATTAACAGCAAGCAGTTGCAGGCGAATATGGAGGTGCAGAAAACCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTcg                                                                                                                                                                                                      <  1:289818/251‑1 (MQ=255)                                             aaGTCTGGGCCGATAACACCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTAGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTcgc                                                                                                                                                                                                     >  1:448995/1‑251 (MQ=255)                                                    ggCCGATAACACCGATTTTGTAAAAGAAGGCGACGTACTGGTCACTCTAGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCAACAGATGATGATTAACAGCAAGCAGTTGCAGGCGCATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCACCCTGATTGGTAGCGacaggc                                                                                                                                                                                              >  2:51194/1‑246 (MQ=255)                                                           aaCACCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAACCTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTGGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAaca                                                                                                                                                                                       <  1:325885/251‑1 (MQ=255)                                                                             gaagGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCAc                                                                                                                                                                     >  1:3893/1‑251 (MQ=255)                                                                                       ggctggtaccgctcacccctaCAGATGCGCGCAAGGAGTTTGAAAAGGCAAAAACTGCCCGTGCTTCCGGCGGTCGCCAAACCCACCAGGTGTGGATTACCACAAAGCAGTGGAAGGGAAATATTGAGGGGTAGAAACTCGCCCCCGGGAAAGCACAAAAGGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACTAGTGCCCAg                                                                                                                                                           <  2:353362/231‑1 (MQ=255)                                                                                        gCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAACAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGCAAATCGCCCTCGCGAAAGCACAACACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGCTTGGTCGCGAAGAGCTGCAACACGCCCGCGCCGCCGTCACCAGTGCCCAgg                                                                                                                                                          >  2:176769/1‑251 (MQ=255)                                                                                                     gACCAGACAGATACTCCGCAGGCGTTTAAAAAAGCCAAACCTGCACTGGTTTTACGGGTTCGCCAAACCCACCAGCTGAGTATAAACAGCAAGCAGTTGCAGGCGAATATTGAGGGGCAGAAAAGCGCCCGCGAGAAAGCACAAAACGACTACACCCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCTCCAGTGCCCAGGCGCAACTGGACGt                                                                                                                                             <  2:404036/251‑1 (MQ=255)                                                                                                         cGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTcgcg                                                                                                                                         <  1:10638/251‑1 (MQ=255)                                                                                                         cGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTcgcg                                                                                                                                         >  1:355185/1‑251 (MQ=255)                                                                                                                  cTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGCCGTCGCGATTccaca                                                                                                                                 >  2:330701/1‑250 (MQ=255)                                                                                                                              ttGAAAAAGCCAAAATTGTACTGGATTCCAGCCATCGCCAAGCCCACGAGCTGCTGATTAACAGCAAGCAGATGCAGGCGCATGTAGACGTGCAGAAAATCGCCCTCGCGATAGCACCAAAAGACGAAAACCGCCGTGTGCCTCTGCGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCGCCAGGGCCAAGACGCAACTGGGCGCCGCGGGTCAACAATAGACTGCCagt                                                                                                                    >  1:149509/1‑249 (MQ=255)                                                                                                                                gAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATc                                                                                                                   <  1:444370/250‑1 (MQ=255)                                                                                                                                        cAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCgatg                                                                                                           >  1:172820/1‑250 (MQ=255)                                                                                                                                          aaaCTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCAGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCCGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGAtt                                                                                                        >  2:385951/1‑251 (MQ=255)                                                                                                                                               gCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAACCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTg                                                                                                     <  1:298834/249‑1 (MQ=255)                                                                                                                                                                    cAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATcagcc                                                                              >  2:324655/1‑251 (MQ=255)                                                                                                                                                                     aaaCCCCCCAGCTGAGTATGAACAGCAAGCAGTTGCAGGCGAATATGGAGGTGCAGAAAATCGCCCCCGCGAAACCACAAACCGACGACATCCGCCGTGTGCCGCGGGGCAATGCCAACCTGATTGGGCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAAGGCCAATCAGGCGATGATTCTGGGGACTAAGCTGGAAGATcagcca                                                                             <  2:375990/251‑1 (MQ=255)                                                                                                                                                                         ccaccaGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGACAGCACAAAACGACTACAACCGCCGTGTGCAGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCACCTATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCccgccg                                                                         >  2:193598/1‑251 (MQ=255)                                                                                                                                                                                    tgatTAACAGCAAGCAGATCCAGGCGAATATTGAGGTGCAGAAAATCGACCTCGCGAAAGCACAAAACGACTACAACAGCAGTCTGCCGCTGGGCACTGCCAACCTGATTGGTAGCGAAGAGCTGCAACACGCCAGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAAATGGAAGATCAGCCAGCCGTGCAACAGGCt                                                              >  2:440039/1‑251 (MQ=255)                                                                                                                                                                                                 gcaGTTGCAGGCGAATATTGAGGAGCAGAAAATCGCCCTAGCGAAAGCACAAAAAGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGGGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTTCAACAGGCAGCCAACGAAGTAc                                                 >  2:245877/1‑251 (MQ=255)                                                                                                                                                                                                     ttGCAGGCGCATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACACCCGCCGTGTGCCGCTGGGCAATACCAACCTGATTGGTCGCGAAGATCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAAGACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGCAGATCCGCCACCCGTGCAACAAGCTGCCACCGAAGTACgga                                              >  2:418431/1‑248 (MQ=255)                                                                                                                                                                                                       tcgggcacatttggAGGTGAAAAAAAGCGCCCCCGCGAAGGCACAAAACGTCAAAAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGGGCAACAGGCTGCCACCGAAGTACGTAAc                                            <  2:19740/237‑1 (MQ=255)                                                                                                                                                                                                                 atTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATg                                                                                                                         >  2:18911/1‑163 (MQ=255)                                                                                                                                                                                                                 atTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATg                                                                                                                         <  1:18911/163‑1 (MQ=255)                                                                                                                                                                                                                               aaTCGCCCTCGCGAAAACACAAAACGACTACACCCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAAGCAGGCGATGATTCTGGGGACTAAACTGGAAGAGCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGTAGCGTACTc                   <  2:206001/251‑1 (MQ=255)                                                                                                                                                                                                                                 tCGCACTCGCGAAAGCAAAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGACAGCCGTGCAACAGGCTGCACCCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGt                 >  2:331818/1‑251 (MQ=255)                                                                                                                                                                                                                                  cGCCCTCGAGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGACAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGAACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGt                 >  1:1151/1‑250 (MQ=255)                                                                                                                                                                                                                                      cTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAAATGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCCCCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGtattat            >  1:272061/1‑251 (MQ=255)                                                                                                                                                                                                                                          cgAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGt        <  1:253845/251‑1 (MQ=255)                                                                                                                                                                                                                                              aGCACAAAACGACGACAACCGCCGTGAGCCGATGGGCACAGACAACCTGATTGGTCACGAAGAGCTGCAACACGCCCGCGGCGCCGTCACCAGTGCCCAGGCGAAACTGGCCGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCTAAGTACGTAACGCCTGGCTTGCGCTGGAGCGAACTCGTATTATCAGTCCg     >  2:87374/1‑250 (MQ=255)                                                                                                                                                                                                                                                cacaAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCgctg  >  2:74874/1‑248 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     GAATCAAATTCAAATTATGTCTCAGGTGTCTGGCAGCGTGACGAAAGTCTGGGCCGATAACACCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAGCGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATG  >  NC_000913/2811603‑2812092

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