Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,334,353	0	G	C	74.4%	76.0 / 26.1	40	T413S (ACC→AGC)	eptA	phosphoethanolamine transferase EptA
Reads supporting (aligned to +/- strand):  ref base G (7/3);  new base C (17/12);  total (25/15)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTTGATAATCCTCCGACAGCCACAGCAGCATCGGCACCTGTTTTTGGCTATCCGGGGCGATGGCATAAGGCAGACCGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGGTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTTGCAGGTTATTGATGTACTCTTCAAGCCCGTGGAACAGCACTTCGTCATAGCATTCGCCGTTGATGCACTGATCAGGTAGATTCAGCGCGGTGACGTTCTGGTGAGGCACGCGGTCGCAGGCACCTTT  >  NC_000913/4334110‑4334603                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                           ttttGATAATCCTCAGACAGCCACAGCAGCATTGGCACCTGTTTTTGGCTATCCGGGGCGATGGCATAAGGCAGACCGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGt                                                                                                                                                                                                                                                     >  1:335441/1‑251 (MQ=255)           cctccGACAGCCACAGCAGCATCGGAACCTGTTTTTGGCTATCCGGGGCGATGGCATACGGCAGCCCGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGGTACAGGTCTGGAtctc                                                                                                                                                                                                                                            >  2:193146/1‑250 (MQ=255)              ccGACAGCCACAGCAGCATTGGCACCTGTTTTTGGCTATCCGGGGCGATGGCATAAGGCAGACCGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCAt                                                                                                                                                                                                                                          <  2:210779/249‑1 (MQ=255)                                               ggCTATCCGGGGCGATGGCATAAGGCAGACCGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTgc                                                                                                                                                                                                       >  1:129231/1‑250 (MQ=255)                                                   aTCCGGGGCGATGGCATAAGGCAGACCGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGGTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGAGGCg                                                                                                                                                                                                   >  2:185364/1‑251 (MQ=255)                                                        gggCGATGGCATAAGGCAGACCGTGCCGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAACGATCACCCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTACACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAAATGCGGCGGAt                                                                                                                                                                                                >  2:185524/1‑249 (MQ=255)                                                              tGGCATAAGGCAGAACGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACAAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCCTTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCgggtg                                                                                                                                                                                        >  2:142924/1‑248 (MQ=255)                                                                   tAAGGCAGACCGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAAt                                                                                                                                                                                    >  1:394712/1‑250 (MQ=255)                                                                    aaGGCAGACCGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAAt                                                                                                                                                                                    >  1:136593/1‑249 (MQ=255)                                                                     aGGCAGACCGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACATGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAgg                                                                                                                                                                                 <  2:394712/251‑1 (MQ=255)                                                                                   gatagatGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGATGGTAAATTTATCCTGATGTTATTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGGCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGTTCGCGGTTGTAATAGGTAGGACCGAGGcgg                                                                                                                                                                   >  2:361002/1‑249 (MQ=255)                                                                                       gatGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGGTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCg                                                                                                                                                                >  1:298131/1‑250 (MQ=255)                                                                                                ttCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGgtgtgt                                                                                                                                                       <  2:136593/250‑1 (MQ=255)                                                                                                   aCCTAACGATTCACCCTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTATTTCAGCAGATTACTCGCTTTATCAAAAATATAGTCGCCGTAAACCAGAGTGTTGTCGTAAGTGTTCAACAGATGCTCTTTGGACCAGGCCTGGAGCTCATTGGTGTCGCAGGATTGGGTAAATTTACTGAACTGAGGCGTATATCGGTTAGACTAGGTCGGACCGTGTCTGCAGATGGTGGGTacg                                                                                                                                                    >  2:143607/1‑248 (MQ=255)                                                                                                          gATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCCGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCAc                                                                                                                                            >  1:388712/1‑251 (MQ=255)                                                                                                                    ggTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCAACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCAccc                                                                                                                                  >  1:404942/1‑250 (MQ=255)                                                                                                                        aGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGTCGGATAGCGGTTGTAATAGTTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCa                                                                                                                              <  1:179506/251‑1 (MQ=255)                                                                                                                                   ccAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATg                                                                                                                   <  1:371896/251‑1 (MQ=255)                                                                                                                                     aGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGt                                                                                                                  >  1:357041/1‑250 (MQ=255)                                                                                                                                         tggtggTAAATTTAACCTGATGTTCTTTCAGCAGATTAATCGCTTTATCACCAATATAGTCGACGTAAACCCGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTActct                                                                                                             >  1:304836/1‑251 (MQ=255)                                                                                                                                             ggtAAATTTATCCGGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGGCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTcc                                                                                                          <  2:335441/250‑2 (MQ=255)                                                                                                                                                aaaTTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGcc                                                                                                      >  1:300881/1‑251 (MQ=255)                                                                                                                                                 aaTTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGccc                                                                                                     <  2:357041/251‑1 (MQ=255)                                                                                                                                                    ttATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCt                                                                                                                                                                    <  1:17956/185‑1 (MQ=255)                                                                                                                                                    ttATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTATTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTACATTTCCTGAACTGCGGCGGATTGCGGTTGTAATAGGTCGGACCGTGGCt                                                                                                                                                                    >  2:17956/1‑185 (MQ=255)                                                                                                                                                     tattctgatgttagtgtatcagATTGATCGCTTTATCAACAATATAGCAGACGTGAACCAGCTTGTTTTCGAAAGTTTTCACCTGTTGCTCTTTGCTACAGTTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGta                                                                                                  <  2:14387/231‑2 (MQ=255)                                                                                                                                                                         agATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGATATTGATGTACTCTTCCAGCCCGTGGCACAGTACTTCGTCATAGc                                                                              >  1:93431/1‑250 (MQ=255)                                                                                                                                                                         agATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTGGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCTGACCGTGGCTGCCGATGGTTTGTAAGACAATCACGCAATCACCCTGGAGGTTATTGATGTACTCTTCCATCCCGTGGAACATTACTTCGTCATAGCa                                                                             <  2:162046/251‑1 (MQ=255)                                                                                                                                                                          gATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGGTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTTGCAGGTTATTGATGTACTCTTCAAGCCCGTGGAACAGCACTTCGTCATAGCAt                                                                            <  2:298131/251‑1 (MQ=255)                                                                                                                                                                             tAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTACGTGTTCACCAGTTGCTCTTTGGTACCGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTTGCCGGTTATTGCTGTACTCTTAAAGCCCGTGGAACAGCACTTCGTCATAGCATTct                                                                         >  2:208256/1‑250 (MQ=255)                                                                                                                                                                              aaTCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGTGGAACAGTACTTCGTCATAGCATTCGc                                                                        >  1:146699/1‑251 (MQ=255)                                                                                                                                                                               aTCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGTGGAACAGTACTTCGTCATAGCATTCGcc                                                                       <  1:125520/251‑1 (MQ=255)                                                                                                                                                                                       aTCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGGTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTTGCAGGTTATTGATGTACTCTTCAAGCCCGGGGAACAGCACTTCGTCATAGCATTCGCCGTTGATg                                                                <  1:388564/250‑1 (MQ=255)                                                                                                                                                                                                      gACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGGTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTTGCAGGTTATTGATGTACTCTTCAAGCCCGTGGAACAGCACTTCGTCATAGCATTCGCCGTTGATGCACTGATCAGGt                                                    >  1:192709/1‑247 (MQ=255)                                                                                                                                                                                                         gTAAAACAGCGTGCTGTCGTAAGCGTTCACCAGTTCCTCTTTGGTACAGGTCTGGATCTCATTGGTGTCGCCGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGCACTAGGTCGGACCGCGGCTGCCGCTGGTGCGTAAGACAATCACGCCATCACCTTGCAGGTTATTGCTGTACTCTTCAAGCCCGTGGAACAGCACTTCGTCATAGCATTCGCCGTTGATGCACTGATCCGGCAGATTc                                              >  2:382466/1‑250 (MQ=255)                                                                                                                                                                                                            aaCCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGATACAGGTCTGGATCTAATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGTGGAACAGTACTTCGTCATAGCATTCGCCGTTGATGCACTGACCAGGCAGATTCAgcg                                          >  1:182128/1‑251 (MQ=255)                                                                                                                                                                                                             aCCAGCGCGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGCTCTCATTGGTGTCGCAGGTTGGGGTAACTTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGTGGAACAGTACTTCGTCCTAGCATTCGCCGTTGCTGCCCTGACCAGGCAGATTCAgcgc                                         >  2:201703/1‑251 (MQ=255)                                                                                                                                                                                                                          cGTAAGTGTTCACCAGTTGCTCTTTGGTACAGGTCTGGATCTCATTGGTGTTGCAGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGCCCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTGGCATGTTATTGATGTACTCTTCAAGCCCGTGGAACAGCACTTCGTCATAGCATTCGCCGTTGATGCACTGATCAGGTAGATTCAGCGCGGTGACGTTCTgg                            <  2:405631/251‑1 (MQ=255)                                                                                                                                                                                                                          cGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGTGGAACAGTACTTCGTCATAGCATTCGCCGTTGATGCACTGACCAGGCAGATTCAGCGCGGTGACGTTCTgg                            >  2:88803/1‑251 (MQ=255)                                                                                                                                                                                                                                        aGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGGGGAACAGTACTTCGTCATAGCATTCGCCGTTGATGCACTGACCAGGCAGATTCAGCGCGGTGACGTTCTGGTGAGGCACGCGGTc              <  1:185524/251‑1 (MQ=255)                                                                                                                                                                                                                                                    gTACAGGTCTGGATCTCATTGGTGTCGCAGGCTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTTGCAGGTTATTGATGTACTCTTCAAGCCCGTGGAACAGCACTTCGTCATAGCATTCGCCGTTGATGCACTGCTCAGGTAGATTCAGCGCGGTGACGTTCTGGTGAGGCACGCGGTCGCAGGCACCttt  >  2:208251/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                           TTTTGATAATCCTCCGACAGCCACAGCAGCATCGGCACCTGTTTTTGGCTATCCGGGGCGATGGCATAAGGCAGACCGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGGTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTTGCAGGTTATTGATGTACTCTTCAAGCCCGTGGAACAGCACTTCGTCATAGCATTCGCCGTTGATGCACTGATCAGGTAGATTCAGCGCGGTGACGTTCTGGTGAGGCACGCGGTCGCAGGCACCTTT  >  NC_000913/4334110‑4334603

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