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	97.7%	131.7 / NA	43	V119I (GTT→ATT)	mdtN	putative multidrug efflux pump membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/25);  minor base G (0/1);  total (17/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCTTCGACCGCCGCACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATT  >  NC_000913/4303518‑4303992                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              ccgctcggACCGCCGCTCGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTTGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAATCCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGTCGTTTTCGACCGGGGCATTAATCGAGTCGGCAcc                                                                                                                                                                                                                                   <  2:53739/243‑1 (MQ=255)  cGCTTCGACCGCCGCACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCa                                                                                                                                                                                                                                  >  1:319033/1‑250 (MQ=255)  cGCTTCGACCGCCGCACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCa                                                                                                                                                                                                                                  >  1:152805/1‑250 (MQ=255)        gACCGCCGCACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGtt                                                                                                                                                                                                                           >  1:153348/1‑251 (MQ=255)         cccgccgCACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTg                                                                                                                                                                                                                          <  2:134712/250‑1 (MQ=255)                       ggCAACTAATGCACCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAacgctac                                                                                                                                                                                                            <  2:138050/251‑1 (MQ=255)                                ggCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGa                                                                                                                                                                                                   <  1:185978/250‑1 (MQ=255)                                         gccgcTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAAttt                                                                                                                                                                                          <  2:91588/251‑1 (MQ=255)                                          ccgcTGACGGCGCTGGCGGCTGACTGCTCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTAAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCATGTTTCGCGGCGGCACGGGACTTTTCTCCCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGGCAAAGCTACGCTGGGTGAGCATAATTTg                                                                                                                                                                                         >  2:140830/1‑251 (MQ=255)                                           cgcTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGc                                                                                                                                                                                        <  1:72152/251‑1 (MQ=255)                                               gACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCtt                                                                                                                                                                                      <  1:12561/249‑1 (MQ=255)                                                     gcTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAgcg                                                                                                                                                                              <  1:120019/251‑1 (MQ=255)                                                         gcggcTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGcgccgc                                                                                                                                                                          <  2:162839/251‑1 (MQ=255)                                                                     gcCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTc                                                                                                                                                              >  1:7344/1‑251 (MQ=255)                                                                      ccTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGATGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGCTGCGCGTCAACTCTACGCTGGGTGAGCATAATTTGCTTATCCAGCCCCGCAAGGGAGGCTTc                                                                                                                                                              >  1:207017/1‑250 (MQ=255)                                                                           aaCAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCtt                                                                                                                                                         <  1:9071/250‑1 (MQ=255)                                                                                    gCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGCGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCCGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAgttt                                                                                                                                               >  2:240627/1‑248 (MQ=255)                                                                                            aTCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCAACAAACTGTTTCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTCTCCATCGCCGCGCGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGt                                                                                                                                       >  2:12920/1‑251 (MQ=255)                                                                                                   tCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCCACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGgcg                                                                                                                                >  1:11215/1‑251 (MQ=255)                                                                                                   tCTGCGCCTCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAATCCTTCTTTCAGTAATAGCTCGGTGCGGCGTAATGTATCTATGGCCTGTTTCGCGGCGGAACGGGCTTTTTCTACCGTGGCATTAATAGCGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGTGCATAATTTGATTATCCAGCGCCGCGAGGGAGGCGTCAGATTTCGCCCGAGTTGCTTCGTTCGGgcg                                                                                                                                >  2:192928/1‑251 (MQ=255)                                                                                                         gcgcgcTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTc                                                                                                                          <  2:329978/251‑1 (MQ=255)                                                                                                                  gcCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGCTGCGGa                                                                                                                  >  1:4338/1‑250 (MQ=255)                                                                                                                    cGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACa                                                                                                               <  1:16371/251‑1 (MQ=255)                                                                                                                          ttcactgTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAAGGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCATCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCa                                                                                                            <  1:246553/242‑1 (MQ=255)                                                                                                                          tGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCaaa                                                                                                          <  1:223537/250‑1 (MQ=255)                                                                                                                          tGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGGGCGGCGTAATGTATCTGTGGCCTGGTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCaaa                                                                                                          <  1:34918/250‑1 (MQ=255)                                                                                                                          tGCACGGTCAACATCTTCCGCTGAGACAAAACCGTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGGCGGCACCAAACTGTTGCGCGTCAACGCTTCGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCGTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCaaa                                                                                                          <  2:197610/250‑1 (MQ=255)                                                                                                                              cGGTCAACATCTGCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTTTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGcc                                                                                                     <  1:298281/251‑1 (MQ=255)                                                                                                                                gTCAACATCTTCCGCTGAGACAAAAACTTCTTTCAGTAATGGCTCGGTGGGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCt                                                                                                   >  1:22278/1‑251 (MQ=255)                                                                                                                                              cTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTgg                                                                                     >  1:349702/1‑251 (MQ=255)                                                                                                                                               tGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGAGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGCGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGACCAGCAAATAGCCCTGTTTGACTGCCTGGt                                                                                    >  2:164633/1‑251 (MQ=255)                                                                                                                                                    cAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTGTGACTGCCTGGTTGTCg                                                                               <  1:111894/251‑1 (MQ=255)                                                                                                                                                                  aGTAATGGCTCGGTGCGGCGTAATGTATCTGGGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTc                                                                 <  1:55875/251‑1 (MQ=255)                                                                                                                                                                        ggCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAAt                                                           >  1:182632/1‑251 (MQ=255)                                                                                                                                                                               tGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAAAGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCg                                                    >  1:228718/1‑251 (MQ=255)                                                                                                                                                                                 cgtcgTAATGTATTTGTGGCCTGTTTCCCGGCGGTACGGGCTTTTTCTACCGTGGCATTAATCGATCCGGCTCCAATCTGTTGCGCGTCAACTCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCt                                                  <  1:174439/248‑1 (MQ=255)                                                                                                                                                                                          gTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTCCCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCgg                                         <  1:240627/251‑1 (MQ=255)                                                                                                                                                                                                           cgcgGCGGCACGGGCTTTTCCCACCGTGGCATTAATCGAGTCGGCACCAAACTTTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTGGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTAt                        <  1:42817/251‑1 (MQ=255)                                                                                                                                                                                                               gcggcACGGGCTTTTTCTCCCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACCCTACTCTGGGTGAGCATAATTTGCTTAGCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTTTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATctgc                    <  1:96594/251‑1 (MQ=255)                                                                                                                                                                                                                  gcACGGGCTTTTTCTACCGGGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTAGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGa                 <  1:317650/251‑1 (MQ=255)                                                                                                                                                                                                                                 aacGTGGCATGAATCTAGTCGGCACCAACCGGTTGCGCGTCAACGCTACGCGGGGTGGGCATAATGTGCTAAGCCAGCGCCGCGAGGGAGTCTCCAGCGTTCGCCAGATTGGCTTCGTGCGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTGGGACTGCCTTGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTGCCGGTGCCACATCAATGGTATCTGCTGACTCGTAAGCGTCa    <  2:59305/247‑1 (MQ=255)                                                                                                                                                                                                                                   cGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCAtt  >  2:86185/1‑249 (MQ=255)                                                                                                                                                                                                                                        cATTAAGCGAGTCGGCACCAATCTGTTGCGCGTCAACGCTACGCTGGGTTAGCATAAGTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGTTCGATGCGGAACAGCAACTCACCCTGTTTGACTGCCTGGGTGTCGGTGACCGCCAGTTCTACAATGc                                                         <  1:342172/189‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              CCGCTTCGACCGCCGCACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATT  >  NC_000913/4303518‑4303992

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