Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,678,043	0	C	T	43.5%	68.2 / 27.7	62	E107K (GAA→AAA)	yphE	putative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (14/21);  new base T (27/0);  total (41/21)
Fisher's exact test for biased strand distribution p-value = 1.99e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGG  >  NC_000913/2677895‑2678192                                                                                                                                                      |                                                                                                                                                      tgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTcc                                                                                                                                                      >  1:3371/1‑150 (MQ=255)  gctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCACCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGc                                                                                                                                                    <  2:41142/151‑1 (MQ=255)   ctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGGCGA‑GGGGTTTCGcc                                                                                                                                                   >  1:858835/1‑151 (MQ=255)     tttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGcc                                                                                                                                                   <  1:356564/149‑1 (MQ=255)         cgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCTACCGAGGCAAAGGGTTTTCcccc                                                                                                                                                   >  2:1091954/1‑146 (MQ=255)          gcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGAGGGAGA‑GGGTTTTCGCCCCTGt                                                                                                                                              >  2:559680/1‑149 (MQ=255)           ccGGGCTTAGCGTTGAACCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTc                                                                                                                                             <  2:438791/149‑1 (MQ=255)           ccGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCAGA‑GGGTTTTCGCCACTGt                                                                                                                                              >  1:39344/1‑148 (MQ=255)                cTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTTCCCCACCGTCAGCCCtt                                                                                                                                      >  1:446557/1‑151 (MQ=255)                   aGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACCGAGGGAGA‑GGGTTTTCGCCCCCGTCAGCCCCTcc                                                                                                                                    >  2:1074986/1‑150 (MQ=255)                       ttGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGA‑GGGTTTCCCCCACTGTCAGCCCCTccacc                                                                                                                                 >  1:1174454/1‑149 (MQ=255)                               aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTTCGCCCCCGTCAGCCCTTCCACCCGACtt                                                                                                                           >  1:954859/1‑147 (MQ=255)                               aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTTCCCCACTGTCAGCCCCTCCCCCACACTTaa                                                                                                                         >  1:32728/1‑149 (MQ=255)                                gttgttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGAGGCCCCGGCCGCGGGAGG‑GGTTTTTCGCCCCCGTTGGCCCCTCCCCCAACCTTAACt                                                                                                                       >  2:1137076/1‑150 (MQ=255)                                   gttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCt                                                                                                                    <  1:651978/150‑1 (MQ=255)                                   gttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCt                                                                                                                    <  2:593215/150‑1 (MQ=255)                                    ttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTCCGCCACTGTCAGCCCCTCCACCAGACTTTACTcc                                                                                                                     >  1:184098/1‑148 (MQ=255)                                     tCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTTCGCCACTGTCAGCCCTTCCCCCAGAATTAACt                                                                                                                       >  1:974775/1‑145 (MQ=255)                                          aCTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAAcc                                                                                                            <  1:864937/151‑1 (MQ=255)                                            tAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCTGACCGAGGGAGA‑GGTTTTTCCCCACCGTCAGCCCCTCCACCAGACTTAAATCCCGATaaaaccc                                                                                                          >  1:372130/1‑147 (MQ=255)                                            tAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCCAGGCAGA‑GGGTTTTCGCCCCTGTCAGCCCCTCCCCCAGACTTAACTCCTGATAAAAcgc                                                                                                          >  1:917155/1‑151 (MQ=255)                                              aCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACcgcgc                                                                                                        <  1:164552/151‑1 (MQ=255)                                                   aaCGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGAGGCAGA‑GGGTTTTCGCCCCCGTCAGCCCCTCCCCCCCACTTAACTCCTGATAAACCCCGCGCAcc                                                                                                   >  2:252672/1‑151 (MQ=255)                                                      gCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGCGGCAGA‑GGGTTTTCGCCCCCGTCAGCCCTTTCACCCGACTTAACTCCTCATAAACCCCGCGCAccccc                                                                                                >  2:50879/1‑151 (MQ=255)                                                         cAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCCACCGCGGGAGA‑GGGTTTCCCCCACTGTCAGCCCCTCCCCCCGACTTAACTCCCGATAAAACCCGCGCcccccc                                                                                                >  2:208225/1‑148 (MQ=255)                                                                   aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGAACGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTTACTCCTGATAAAACCCGCGCACCCCCCGTTCAGCGGc                                                                                     >  1:735081/1‑149 (MQ=255)                                                                   aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCCGA‑GGGTTTTCGCCCCCGTCAGCCCCTCCACCCGACTTAAATCCTGATAAAACCCGCGCACCCCCCGATCAGCGGCAc                                                                                   >  2:736897/1‑151 (MQ=255)                                                                   aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCAGA‑GGGTTTTCGCCCCCGTCAGCCCCTCCACCCAACTTAACTCCTGGTAAAACCCGCGCAaccccc                                                                                               >  2:973837/1‑138 (MQ=255)                                                                    aGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACg                                                                                  <  2:181241/151‑1 (MQ=255)                                                                          gTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGAGGGAGA‑GGGTTTCCGCCACTGTCAGCCCTTCCCCCCGACTTAACTCCTGATAAACCGCGCGCACCCCCCCTTCAGCGGCACGGCCAg                                                                             >  1:364784/1‑150 (MQ=255)                                                                          gTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGt                                                                            <  1:1189193/151‑1 (MQ=255)                                                                           ttGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGt                                                                            <  1:756470/150‑1 (MQ=255)                                                                            tGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGAGGCGGA‑GGGTTTTCGCCACCGTCAGCCCTTCCCCCCGACTTTACTCCTGATAAAACCCGCGCCCCCCCAGCTCCACGGCACGGCggg                                                                             >  2:812246/1‑146 (MQ=255)                                                                              ggCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCCACCGAGGCAGA‑GGGTTTTCGCCACCGTCAGCCCCTCCCCCAGACTTAACTCCTGATAAAACGCGCGCCCCCCCCGCTCAGCGGCCCGGCGAGGCAGc                                                                        >  2:1074659/1‑151 (MQ=255)                                                                              ggCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGc                                                                        <  1:660085/151‑1 (MQ=255)                                                                                 aTCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCCACCCGGGGGGA‑GGGTTTTCGCCCCTGTCCGCCCCTCCCCCAGACTTTACTCCTGATAAAACGCGCCCCCCCCCCGTTCAGCGGCACCGCGGGTCcgcg                                                                       >  2:25567/1‑149 (MQ=255)                                                                                  tCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTa                                                                     <  1:728712/150‑1 (MQ=255)                                                                                    cTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGCGGCAGA‑GGGTTTTCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGCTCAGCGGCCCGGCGAGGCAGCGTAg                                                                    >  2:378515/1‑149 (MQ=255)                                                                                       cgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAAATCCTGATAAACCGCGCGCACCCCCACTTCAACGGCACCGCGAGTCAGCGTAGCTTc                                                                >  1:1029855/1‑150 (MQ=255)                                                                                         ccATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGAGGGCGA‑GGGTTTTCGCCACTGTCAGCCCTTCCCCCCGACTTAACTCCCGATAAAACCCCCGCACCCCCCGTTCAGCGGCCCGGCGGGGCCGCGGAGCtttg                                                               >  2:1025455/1‑147 (MQ=255)                                                                                         ccATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTc                                                             <  2:801062/151‑1 (MQ=255)                                                                                         ccATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTc                                                             <  2:777238/151‑1 (MQ=255)                                                                                         ccATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTc                                                             <  1:220791/151‑1 (MQ=255)                                                                                              tGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCTGAACGAGGCAGA‑GGGTTTTCCCCACTGTCAGCCCCTCCCCCCCACTTTAATCCTGATAAAACGCCCGCACCCCCCATTCAGCGGCACGGCGAGTCAGCGGAGCTTTGTCTcc                                                          >  1:49986/1‑149 (MQ=255)                                                                                                 aGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTcc                                                      <  2:111821/150‑1 (MQ=255)                                                                                                 aGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGGGGGAGA‑GGGTTTTCCCCCCTGTCAGCCCCTCCACCCGACTTAACCCCTTATAAAACCCGCGCCCCCCCCGTTCCGCGGCACGGCGGGGCAGCGTAGACTCGTCACCCTccc                                                     >  2:1130562/1‑150 (MQ=255)                                                                                                   gTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGAGGGCGA‑GGGTTTTCGCCCCTGTCACCCCCTCCACCCGACTTAACTCCCGATAAACCCCGCGCACCCCCCGGTCCGCGGCACGGCGCGTCAGCGTAGCTTTGTCCCCCTcc                                                      >  2:418982/1‑148 (MQ=255)                                                                                                     aatcaatcATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTc                                                  <  1:556716/150‑1 (MQ=255)                                                                                                     aatcaatcATGCCGTTGCGGCGGGGCCACTGAACGAGGGGGA‑GGGTTTTCGCCCCTGTCAGACCTTCCACCACACCTTACCCCTGCTAAACAACGCGCACCCCCCGTTCAGCGGGACGGCGCGGCAGCGGAGCTTCGTCACCTTCCCGTc                                                  >  2:793515/1‑150 (MQ=255)                                                                                                         aatcaTGCCGTTGCGGCGGGGCCCCCGACCCCGGCAGA‑GGGTTTTCGCCCCCGTCCGCCCCTCCACCCAACCTAACTCCTGGTAAAACGCGCCCACCCCCCCTTCAACGGGACGGCGGGGTAGCGGAACTTCGTCACCCTCCCGTc                                                  >  2:926319/1‑146 (MQ=255)                                                                                                            caTGCCGTTGCGGCGGGGCCACTGACCGCGGCAGA‑GGGTTTTCGCCCCTGTCAGCCCCTCCACCCGACTTAACTCCTGATAAACCCCCCGCACCCCCCGTTCAGCGGCCCGGCGAGTCAGCGTAGCTTCGTCACCCTCCACTCGCGgctc                                           >  2:1059724/1‑150 (MQ=255)                                                                                                              tGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTTCGCCCCTGTCAGCCCCTCCACCAGACCTAACTCCCGATAAAACCCGCGCACCCCCCGTTCCGCGGCACGGCGGGGCAGCGGAGCTTCGTCACCTTCCCGTCGCGTCTCGcc                                        >  2:206116/1‑151 (MQ=255)                                                                                                                 cGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAAt                                     <  1:1076052/151‑1 (MQ=255)                                                                                                                     gcggcgGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCa                                  <  1:559696/150‑1 (MQ=255)                                                                                                                     gcggcgGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCa                                  <  2:1108079/150‑1 (MQ=255)                                                                                                                      cggcggGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATcc                                   >  1:298072/1‑148 (MQ=255)                                                                                                                        gcggGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGata                              <  2:578661/151‑1 (MQ=255)                                                                                                                        gcggGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGat                               >  2:957055/1‑150 (MQ=255)                                                                                                                         cggGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGatat                             >  2:770101/1‑151 (MQ=255)                                                                                                                               cACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGc                       >  1:1179352/1‑151 (MQ=255)                                                                                                                                    aCCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATcc                  >  1:61806/1‑151 (MQ=255)                                                                                                                                       gAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCggg               >  2:287131/1‑151 (MQ=255)                                                                                                                                                     tCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCgg  >  1:267858/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGG  >  NC_000913/2677895‑2678192

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