Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,174,610	0	A	C	44.3%	97.9 / 14.2	62	W186G (TGG→GGG)	ycfT	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (25/9);  new base C (0/27);  total (25/37)
Fisher's exact test for biased strand distribution p-value = 7.09e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCAT  >  NC_000913/1174462‑1174757                                                                                                                                                     |                                                                                                                                                    cAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccacc                                                                                                                                                   >  2:513351/1‑151 (MQ=255)  aGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccacc                                                                                                                                                   >  1:291267/1‑150 (MQ=255)   ggCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCcaccac                                                                                                                                                 <  2:859797/151‑1 (MQ=255)       ccgaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGc                                                                                                                                              <  1:179626/150‑1 (MQ=255)          accgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCCCACTGTTCATTCCCCACCACGGCGTg                                                                                                                                           <  1:543577/150‑1 (MQ=255)          accgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTgg                                                                                                                                          <  1:179077/151‑1 (MQ=255)           ccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTggg                                                                                                                                         <  1:172703/151‑1 (MQ=255)                aGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGa                                                                                                                                    <  1:858789/151‑1 (MQ=255)                 gAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGa                                                                                                                                    >  2:439779/1‑150 (MQ=255)                    cGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAAtt                                                                                                                                >  1:94129/1‑151 (MQ=255)                         aGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCg                                                                                                                           >  1:332919/1‑151 (MQ=255)                            aaaGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGcc                                                                                                                         >  2:941663/1‑150 (MQ=255)                                aagCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAACCCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCCCCACGGGGTGGGAACGAAATTAACCGCCacac                                                                                                                     <  2:17783/150‑1 (MQ=255)                                   ccATCGGCAGATGGCGGCGAAACGGCACCTTTTTAACACAGGTCATTATTGTTGCGCAACACAAGCCGCAAGGGTGGAAAAAAGGAAAATTGCGGATCCCACTTTTCATCCCCCCCCCCGGCGGGGAAACAAAATAAGCCGCCACACTca                                                                                                                  <  2:143509/150‑1 (MQ=255)                                       cagcagATGGCGGCGCAACGGCACCTCTTAAACACAGGTCATTATTTTTGCCCAAAACCATGCGCAAAGGCTGTAATAAGGAAAATTGCGGATCACCCTGTTCATTCCCCCCCCGGGCGGGGGACCGAAATTAACCGCCACACTcag                                                                                                                 <  2:701270/147‑1 (MQ=255)                                            gATGGCGGCGCAACGGCACCTCTTAAACACAGGTCATTATTGTTGCGCAAAACCATGCGCCAAGGCTGTAAAAAGGCAAATGGGGGATCACACTGTTCATCCCCCCCCCCGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTac                                                                                                           <  1:293834/148‑1 (MQ=255)                                             aTGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTacaaac                                                                                                       >  2:1003272/1‑151 (MQ=255)                                              tGGCGGCGCAACGGCACCTCTTAAACACAGTTCATTATTGTTGCGCAAAACCAGGCGCCAGGGCTGTAATAAGGAAAATGGCGGATCACACTTTTCATCCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTacaaac                                                                                                       <  1:881336/150‑1 (MQ=255)                                               ggcggcGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTacaaaca                                                                                                      >  1:1069189/1‑150 (MQ=255)                                               ggcggcGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCAAAACCATGCGCCAAGGCTTTAATAAGGCAAATTGCGGATCACCCTGTTCATTCCCCCCCACGGCGGGGGAACGAAATTAACCGCCACACTCAGCAGTacaa                                                                                                         <  1:155016/147‑1 (MQ=255)                                               ggcggcGAAACGGCCCCTCTTTAACACAGGTCATTTTTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGAAAATTGCGGATCACCCTTTTCATCCCCCCCCCCGGCGGGGGAACGAAATTAACCGCCACACTCAGCAGTacaaaca                                                                                                      <  2:240034/150‑1 (MQ=255)                                                 cgtccccgcGGCACCTCTTAACAAAAGGTCTTTATTTTTGCGCCACCCTGTGGGCCAGGGCTATAAAAAGGGAAAGTGCGGACCCCCTTGTTCATCCCCCCCCACGGCGTGGAAAGAAATTTAACCGCCACACTCAGCAGTACAAGCAAgg                                                                                                   <  2:631193/143‑1 (MQ=255)                                                  ggcgCAACGGCACCTCTTAAACACAGGTCATTATTGTTGCGCAAACCCATGCGCCAGGGCTGTAAAAAGGAAAATGGCGGATCCCACTGTTCATCCCCCCCCACGGGGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGc                                                                                                  <  2:1021502/151‑1 (MQ=255)                                                  ggcgAAACGGCACCTCTTAAACACAGGTCATTATTGTTGCGCCAACCCATGCGCCAAGGCTGAAATAAGGAAAATTGCGGATCACACTGTTCATTCCCCCCCACGGCGTGGGAACGAAATAAACCGCCACACTCAGCAGTACAAACAAg                                                                                                    <  1:568463/149‑1 (MQ=255)                                                       aaggcctccTCTTTAACACAGGTCATTATTGTTGCCCAAAACCACCCGCCAAGCCTGTAAAAAGAAAAATTGCGGATCCCACTGTTCATCCCCCTCCAGGGGGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAAt                                                                                              <  1:9057/143‑1 (MQ=255)                                                       aaCGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAAt                                                                                              >  1:630200/1‑150 (MQ=255)                                                          ggCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAGGCGCCAAGGCTGTAATAAGGAAAATTGCGGATCACCCTTTTCATTCCCCACCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTg                                                                                          <  1:514503/151‑1 (MQ=255)                                                                 ctTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAgcgc                                                                                   >  1:699376/1‑151 (MQ=255)                                                                  tttAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                                                  >  1:587918/1‑151 (MQ=255)                                                                   ttAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATGGCGGATCACACTTTTCATTCCCCACCAGGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAgcg                                                                                    <  1:567790/148‑1 (MQ=255)                                                                   ttAACACAGGTATTTATTGTTGCGCCAAACCATGCGCCAAGTCTGTAAAAAGGAAAATGGGGAATCACCCTGTTCATCCCCCCCCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                                                  <  2:773664/150‑1 (MQ=255)                                                                     aaCACAGTTCTTTTTTGTGGCCCAAACCCTTGCGCAATGGCTGTAATAAGGAAAATGGCGCATCCCTCTTTTCGTCCCCCCCCGGGGCGGGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAg                                                                                <  2:623131/150‑1 (MQ=255)                                                                      acacAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGAc                                                                              >  2:57588/1‑151 (MQ=255)                                                                        acaGTTCATTTTTGTTGCGCAAAACCATGCGCCAAGGTTGAAAAAAGCCAAATGGGGGATCACACTGTTCATCCCCCCCCCCGGGGGGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACgg                                                                            <  2:548307/151‑1 (MQ=255)                                                                          agttCATTATTGTTGCGCAAAACCAGGCGCCAGGTCTGAAATAAGGCAATTTGCGGATCACACTGTTCATCCCCCCCCCCGGCGTGGAAAGAAAATGAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGc                                                                           <  2:583032/147‑1 (MQ=255)                                                                           ggTCATTTTTGTTGCCCAAAACCAGGCGCCAAGGCTGTAATAAGGAAAATTGCGGATCACACTGTTCATTCCCCCCCACGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCt                                                                          <  2:242078/150‑1 (MQ=255)                                                                                 tattGTGGCCCCAAACCATGCGCAAAGGCTTTAAAAAGGCAAATTGCGGATCACACTGTTCATCCCCCCCCCCGGCGTGGAACGGAAATTAGCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAt                                                                    <  1:273758/150‑1 (MQ=255)                                                                                   ttgttgCGCAAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCTTTCCCCCCCCCGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAttt                                                                  <  1:490039/150‑1 (MQ=255)                                                                                   ttgttcCGCCAAACCATCCGCCAGGGCTATAATAAGTAAAATTGCGGATCACACTTTTCATCCCCCCCCCCGGCTGGGAAACAAATTTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAttt                                                                  <  1:340403/150‑1 (MQ=255)                                                                                         cgcCAAACAATGCGCAAGGGCTGTAAAAAGGCAAATTGCGGATCCCACTGTTCATCCCCCCCCCCGGCGTGGAAACGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACAt                                                             <  2:951143/149‑1 (MQ=255)                                                                                           caaacccATGCGCCAAGGCTGTAATAAGGAAAATTGCGGATCACACTGTTCATCCCCCCCCACGGGGGGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACAt                                                             <  2:505842/142‑1 (MQ=255)                                                                                             aaaCCATGCGCAAAGGCTGAAATAAGGCAAATTGGGGATCACACTTTTCATCCCCCCCCACGGCGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACg                                                       <  2:849768/151‑1 (MQ=255)                                                                                              aaCCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACg                                                       >  1:609025/1‑150 (MQ=255)                                                                                                ccATGCGCCAAGGCTGAAATAAGGCAAATTGCGGACCACCCTTTTCATCCCCCACCACGGCGGGGAAACGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGaa                                                     <  2:415741/150‑1 (MQ=255)                                                                                                 cATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAAt                                                   >  2:792851/1‑151 (MQ=255)                                                                                                  aTGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAAt                                                   >  2:396510/1‑150 (MQ=255)                                                                                                  aTGCGCCAAGGCTGTAATAAGACAAATTGCGGATCACCCTGTTCATTCCCCCCCCCGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATa                                                  <  1:1007099/151‑1 (MQ=255)                                                                                                  aTGCGCCAAGGCTGAAAAAAGGCAAAGTGCGGATCACCCTGTTCATCCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAAt                                                   <  2:134895/150‑1 (MQ=255)                                                                                                   tGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAAt                                                   >  1:800795/1‑149 (MQ=255)                                                                                                       ccAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACa                                              >  1:1022428/1‑150 (MQ=255)                                                                                                         aaGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGa                                                <  2:1067820/146‑1 (MQ=255)                                                                                                                taataaGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTTAGAACGAAATTAACCGCCACACTCAGCAGAACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCAt                                    >  1:74980/1‑151 (MQ=255)                                                                                                                 aataaGGCAAATGGCGGATCACACTGTTCTTTCCCCCCCCCGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCAt                                    <  1:309962/150‑1 (MQ=255)                                                                                                                  ataaGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATAc                                  >  1:730617/1‑151 (MQ=255)                                                                                                                    aaGGAAAATTGCGGATCACTCTATTCATCCCCCCCCGCGGCGGGGGAAGGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCAt                                    <  1:450560/147‑1 (MQ=255)                                                                                                                        cAAATTGCGGATCACACTTTTCTTCCCCCCCCACGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATa                             <  1:699129/150‑1 (MQ=255)                                                                                                                        cAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAc                            >  2:861871/1‑151 (MQ=255)                                                                                                                         aaaTTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAcc                           >  1:871049/1‑151 (MQ=255)                                                                                                                              gCGGATCACCCTGTTCATTCCCCACCACGGCGTGGGAAGGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACaa                       <  2:335777/150‑1 (MQ=255)                                                                                                                                  aTCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTgg                  >  2:447430/1‑151 (MQ=255)                                                                                                                                  aTCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTgg                  >  1:555217/1‑151 (MQ=255)                                                                                                                                   tCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTgg                  >  2:353614/1‑150 (MQ=255)                                                                                                                                       acTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTg             >  2:79776/1‑151 (MQ=255)                                                                                                                                        cTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTg             >  1:84233/1‑150 (MQ=255)                                                                                                                                                 ccccccccCGGCGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGatc    <  1:823581/150‑1 (MQ=255)                                                                                                                                                  cccACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGatcat  >  1:115793/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCAT  >  NC_000913/1174462‑1174757

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