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	92.7%	109.4 / NA	42	S141N (AGC→AAC)	emrA	multidrug efflux pump membrane fusion protein EmrA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (22/16);  minor base C (0/3);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 9.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAATCAAATTCAAATTATGTCTCAGGTGTCTGGCAGCGTGACGAAAGTCTGGGCCGATAACACCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAGCGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATG  >  NC_000913/2811602‑2812092                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     ggAATCAAATTCAAATTATGTCTCAGGTGGCTGACAGCGTGACGAAAGTCTGGGCCGATAACACCGATTGTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTGTGAAAAAGCCAAAACTGCACTGTCTTCCAGCGTTCGCCAAACCCACCAGCGGATGATTAACAGCAATCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGAc                                                                                                                                                                                                                                                  <  2:178980/251‑1 (MQ=255)                 tatttccaagttgccTGGCAGCTTGACGAAAGTCTGGGCCGATAACACCGTTTTTGTAAAAGAAGGCGACGTGCTGGTCCCTCTCGACCCGACAGATGCTCGCAAGGCTTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTACCAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGtgtg                                                                                                                                                                                                                                   <  1:269181/236‑1 (MQ=255)                     ccccaggtgccgggaggCGTGAGGAAGGTCTGGGCCGATAACCCGGATTTGGTAAAAGAAGGCGACGTGCTGGCCACTCCCGACCCGACAGATGCTCGCCAGGCTTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGAGGATTAACAGCAAGCAGTTGCAGGCGAATATGGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCt                                                                                                                                                                                                                              <  2:186826/235‑1 (MQ=255)                            tgtCTGGCATCGGGACGAAAGTCTGGGCCGATAACCCCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTGGCAGGCGAATATTGAGGTGCAGAAAACCGCCCCCGCGAAAGCACAAAACGACTCCAAccgcc                                                                                                                                                                                                                                        <  2:77258/234‑1 (MQ=255)                            tgtCTGGCAGCGTGACGAAAGTCTGGGCCGATAACACCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAAccgcc                                                                                                                                                                                                                                        >  1:77258/1‑234 (MQ=255)                             gtCTGGCAGCGTGACGAAAGTCTGGGCCGATAACACCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATg                                                                                                                                                                                                                      <  1:114322/251‑1 (MQ=255)                                       gTGACGAAAGTCTGGGCCGATAACACCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGAt                                                                                                                                                                                                            >  2:5537/1‑251 (MQ=255)                                          aggAAAGTTTGGGCCGATAACACCGATTTTGTAAAAGAAGGCGACGTGCGGGTCCCTCTCGCCCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTGCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCTAAAGCACAAAACGACTTCATCCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTgg                                                                                                                                                                                                         <  2:357604/249‑1 (MQ=255)                                           cGAAAGTCTGGGCCGATAACACCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGt                                                                                                                                                                                                        >  1:246468/1‑251 (MQ=255)                                             aagggcggggcagaaaaacccGATTTGGAAAAAGAAGGCGCCGTGCGGGACACTCGCGACCGGAAAGGTGCTCGCAAGGCGTGTGAAAAAGCCAAAACTGCACGCGCTGCCAGCGTTGGCCAAACAAACAAGCTGCGGATAAACAGCAAGCAGGGGCAGGCAAATATTGGGTTGAAGAAAATCGCCCTCGCAAAACCACAAACCGACTACAACCGCCGTGGGCCGCTGGGCAATGCCAACCTGATTGGTcg                                                                                                                                                                                                      <  2:285990/232‑1 (MQ=255)                                                               cccGATTTTGTAAAAGAAGGCGACGTGCGGGTCCCTCTCGCCCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCAAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTGGCAGGCGAATATTGAGTGGCAGAAAATCGCCCTCGCGAAAGCAAAAAGCGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGc                                                                                                                                                                                    <  2:110697/250‑1 (MQ=255)                                                                                    gACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAACACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAACACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGc                                                                                                                                                               >  2:277347/1‑251 (MQ=255)                                                                                                cctctCGACCCGACAGATGCTCGCCAGCCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAACCCCCCCAGCTGATGATTAACAGCAAGCAGTGGCAGGCGAATATTGAGGTGCAGAAAACCGCCCGCGCAAAACCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACt                                                                                                                                                   <  2:194446/250‑1 (MQ=255)                                                                                                  ccccGACCCGACAGAGGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACGGGCTTCCAGCGTTCGCCAAACCCACCAGCGGATGATTAACAGAAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCAAAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTgg                                                                                                                                                 <  1:54556/248‑1 (MQ=255)                                                                                                   ctcGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGCTGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGAAGAAAATCGCCCTAGCGAAAGCACAAAACGACTACAACCGCCGGGTGCCGCTGGGCAATGCCAACCTGATTGGACGCGAAGAGCTGCACCACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCACCTgg                                                                                                                                                 >  2:310072/1‑250 (MQ=255)                                                                                                   ctcGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACGGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTgg                                                                                                                                                 >  2:164018/1‑250 (MQ=255)                                                                                                             cAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGAtt                                                                                                                                      >  2:176198/1‑251 (MQ=255)                                                                                                                        cAGGCGTTGGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTCCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGCCTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAacaatacaa                                                                                                                           <  1:5537/251‑1 (MQ=255)                                                                                                                                tGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATc                                                                                                                   >  1:127090/1‑251 (MQ=255)                                                                                                                                  aaaaaGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGCGGAGCAGAAAATCACCCTCGCGAAAGCACAAAACGACTACAACCGACGAGAGCCGCAGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACAAGCCCGCGACGCCGTCACCAGTGCCCAGGCGCACCTGGACGTCGCGATTCAACAATACAGTGCCAATc                                                                                                                   >  1:281771/1‑249 (MQ=255)                                                                                                                                   aaaaGCCAAATCTGCACTGGCTTCCAGCGTTCGCCAACCCCACCAGCGGATGATTAACAGCAGGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCGGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAgg                                                                                                                <  1:216991/251‑1 (MQ=255)                                                                                                                                    aaaGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAAGGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGCCGCGCTTCAACAATACAATGCCAATCAGGc                                                                                                               >  2:308094/1‑251 (MQ=255)                                                                                                                                        ccaaaacgccccTGGCTTCCAGGGGACGCCAACCCCACCAGCTGATGATTAGCAGCAAGCAGTTGCAGGCGAAAATGGAGGTAGAGAAAATCGCCCCCGCGAAAGCACAAACCGTCCCCAACCGCCGCGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCTCCAGTGTCCAGGCGCAACTGGACTTCGCGATTCAACAATACAATGCCAATCAGGCgatg                                                                                                           <  2:307371/240‑1 (MQ=255)                                                                                                                                               tgaaCTGGCTTCGAGCGTTAGCCAAAACAACCAGCTGATGATTAACAGCAAGGCGCTGCAGGAGAATATTGAGGAGGAGAAAATCGCCCTAGCGAAAGAAGAACACGACTACAACAGCCGTGTGCCGCTGGGCAATGCCCACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGCCGCAGTCACCACTGCCCAGGCGAATCGGGACGTCGCGATTCAACAATACAATGACAATCAGGCGCTGATTATgg                                                                                                    >  2:136832/4‑251 (MQ=255)                                                                                                                                                               gTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGa                                                                                       >  2:271609/1‑248 (MQ=255)                                                                                                                                                                                          aaCAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGTCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCAc                                                         <  1:58251/251‑1 (MQ=255)                                                                                                                                                                                             agcaagcaGTTGCAGGCGAAAATGGAGGTGAAGAAAACCGCCCTCGCGAAAGCACAAAACGACTACAGCCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACGGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGa                                                      <  2:181923/251‑1 (MQ=255)                                                                                                                                                                                                aagcaGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTAAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCCCCGAAGt                                                   >  1:276482/1‑251 (MQ=255)                                                                                                                                                                                                 agcaGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAga                                                   >  2:22661/1‑249 (MQ=255)                                                                                                                                                                                                        gaaggcaaaaagtGAGGTGCAGAAAATCGCCACCGCAAAAGCAAAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACg                                           <  2:205336/239‑1 (MQ=255)                                                                                                                                                                                                         cAGGCGACTATTGAGGTGCAGAAAATCGCCCTCGCGAAAGAACAAAACGAATAAAACCGCAGTGTACCGCTGGGAAATGCCAACATGCTTGGCCGCGAAGAGATGCCACACGCCCGCGACGCAGTCACCACAGCCCAGGCGCAACTGGCCGTCGCGATTCAACAATACGATGACAATCCGGCGATGCTTCTGGGGACTAAACTGGAAGCTAAGCCAGCCGAGCAACAGGATGCGACCGAAGTACGTAACGc                                          >  2:344030/1‑251 (MQ=255)                                                                                                                                                                                                          aGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGcc                                         >  1:281161/1‑251 (MQ=255)                                                                                                                                                                                                             cGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCCCCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCgatg                                                                                                           >  2:104454/1‑182 (MQ=255)                                                                                                                                                                                                             cGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCgatg                                                                                                           <  1:104454/182‑1 (MQ=255)                                                                                                                                                                                                                atattagggcgcagacactCGCCCTGGCGAAAGCACAAAACGACTCCCCCCGCCGTGCGCCACAAGGCAATGCCAACCTGATTGATCCCGAAGAGATGCAACACGCCCGCGACGCCGCCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCCACAATACAATGCCAATCAGGCGATGATTCTGGGGCCTAAACTGGAAGATCAGCCAGCCGTGCCACAGGCTGCCACCGACGTACGTAACGCctggccg                                   >  2:22516/19‑249 (MQ=255)                                                                                                                                                                                                                  atggAGGTGCAGAAACTCGCCCTCGCAAAACCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACATGCTGCCACCGAAGTACGTAACGCCtggctggc                                 <  1:133397/248‑1 (MQ=255)                                                                                                                                                                                                                     gAGGTGCAGAAAACCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCt                              <  1:22661/251‑1 (MQ=255)                                                                                                                                                                                                                       ggTGCAGAAAATCGCACTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGCCGCTgg                            >  1:275826/1‑251 (MQ=255)                                                                                                                                                                                                                                aaTCGCCCTCGCGAAAGCACAAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGAACGTAACGCCTGGCTGGCGCTGGAGCGTACTc                   >  2:83573/1‑251 (MQ=255)                                                                                                                                                                                                                                      ccTCGAGAAAGCACAAAACGACTACAACCGCCGAGGGCCGCTGGGCAATGCCAACATGATTGGTAGCGAAGAGCTGCAACACGCACGAGAAGCCGACACCAGTGCCCAGGCGCAAATGGACGTCGCGATACACAAATACAATGCCAATAAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGTCAGCCGTGCAACAGGCTGCCACCGAAGTAAGTAACGCCTGGCTGGCGCTGGAGCGTACTCGaatt              >  1:273242/1‑250 (MQ=255)                                                                                                                                                                                                                                          gcgAAAGCACAAACCGCCTACACCCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTTGCGCTGGAGCGTACTCGTATTATCAg         <  2:241585/251‑1 (MQ=255)                                                                                                                                                                                                                                                  acaAAACGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATg  >  1:54987/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     GGAATCAAATTCAAATTATGTCTCAGGTGTCTGGCAGCGTGACGAAAGTCTGGGCCGATAACACCGATTTTGTAAAAGAAGGCGACGTGCTGGTCACTCTCGACCCGACAGATGCTCGCCAGGCGTTTGAAAAAGCCAAAACTGCACTGGCTTCCAGCGTTCGCCAAACCCACCAGCTGATGATTAACAGCAAGCAGTTGCAGGCGAATATTGAGGTGCAGAAAATCGCCCTCGCGAAAGCACAAAGCGACTACAACCGCCGTGTGCCGCTGGGCAATGCCAACCTGATTGGTCGCGAAGAGCTGCAACACGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATG  >  NC_000913/2811602‑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: ‑