Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,947,481	1	.	A	27.8%	173.4 / 102.7	72	intergenic (+18/‑16)	metZ/metW	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base . (9/43);  new base A (0/20);  total (9/63)
Fisher's exact test for biased strand distribution p-value = 5.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01

CGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGG  >  NC_000913/2947330‑2947603                                                                                                                                                          |                                                                                                                           ggtcttAAAAAAAAGGCTGAATGAAAACCACAGCGGAGTATTCTCCGAATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGt                                                                                                                             <  6:1327288/146‑1 (MQ=14)  ggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCAGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGT‑a                                                                                                                           <  5:1067154/151‑1 (MQ=255)  ggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAGGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGt                                                                                                                             <  5:1382182/151‑1 (MQ=255)  ggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGt                                                                                                                             <  6:1389199/151‑1 (MQ=255)  ggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGt                                                                                                                             <  3:199591/151‑1 (MQ=255)  ggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGt                                                                                                                             <  2:933112/151‑1 (MQ=255)  ggcgTAAAAAACCGGCTGCAGGAAAACGCGAGCGGAGTATAGTGCGCATCCCGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGT‑a                                                                                                                           <  3:1433370/148‑1 (MQ=14)   tcagatctcaaaGGTTGCATGAATACGCGAGCGGTGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGT‑a                                                                                                                           <  3:447016/141‑1 (MQ=17)   gggTAAAAATAAGGTTGCATGAAATAGATATGGGAGTATATTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGT‑aa                                                                                                                          <  1:958241/151‑1 (MQ=14)   gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑a                                                                                                                           <  7:1346386/151‑1 (MQ=255)   gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑a                                                                                                                           <  1:1016392/151‑1 (MQ=255)    ggTAACAAGACGGGTGCAGGAAAAAGCGAGCGCTGTATAGTGCGCATCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGT‑aaa                                                                                                                         <  3:20711/151‑1 (MQ=14)    ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGGGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑aa                                                                                                                          <  7:611442/151‑1 (MQ=255)     gaaaaaacacccgtGCAGGAAAACGCGAGAGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGT‑aaaa                                                                                                                        <  3:987503/138‑2 (MQ=14)     gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑aaa                                                                                                                         <  3:964792/151‑1 (MQ=255)      tataattaaGGTTGCATGAAATCGAGAGTGTTGTATTCTTGGCTTCCTCGGACGGGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGT‑aaaa                                                                                                                        <  1:1565964/143‑2 (MQ=11)      tAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTT‑AT‑GAAGT‑a                                                                                                                           >  7:534724/1‑147 (MQ=255)      tAAAAAAAAGGTTGCATCAAAAAGAGTGCGGAGTATAGTCCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGT‑a                                                                                                                           <  6:532396/147‑1 (MQ=17)      aaaaagacaGGTGGAAGGAAAAAGCCAGCGGAGTATTGTCCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGT‑aaaag                                                                                                                       <  3:1447510/150‑3 (MQ=14)       aaaaaaacGGTTCCAGGAAAACGAGAGCGCTGAATACTCGGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGT‑aaaagc                                                                                                                      <  6:1584421/151‑4 (MQ=11)       aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGc                                                                                                                       <  5:1037858/151‑1 (MQ=255)       aaaaaaaaGGGTGCAGGAAAACGCTAGCGCAGTATACTGCGACTCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGT‑aaaag                                                                                                                       <  2:182114/150‑3 (MQ=14)       aaaaaaaaGGGTGAATCAATACCCGATCGTAGTATACTCCGCATACACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTAAT‑GAAGT‑aaaagc                                                                                                                      <  5:1462133/151‑4 (MQ=11)        ataattaGGTTGGATGAAAAGGCGAGCGGAGTATAGTGCGGATCGTGGGACGCGGGGTGGAGCAGCGTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAAGCAATTAAAATTTGAT‑GATGT‑AAAgga                                                                                                                      <  8:773549/145‑3 (MQ=17)        aaaaaacGGTTGCATGAAAATCCCAGCGTAGTATTTTGGGATTCCCGGGACGGGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGT‑aaaagca                                                                                                                     <  4:1498171/151‑5 (MQ=1)        aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGGCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCa                                                                                                                      <  4:696639/151‑1 (MQ=255)         aaatacggTTGCATGAAAACTCGATCGGAGTATACTCGGATTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGT‑aaaagcag                                                                                                                    <  7:565740/145‑6 (MQ=11)         aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAATTTTAAT‑GAGGT‑AAAGc                                                                                                                       >  5:35971/1‑149 (MQ=255)         aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAg                                                                                                                     <  6:1271823/151‑1 (MQ=255)          ataacggTTGCATGAATAATAGATCGGTGAATTGTGCGCTTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGt                                                                                                                    <  3:206467/146‑1 (MQ=11)          accaagcgtGCATGAAAAATAGAGCGCAGTATAGTGGGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGt                                                                                                                    <  6:201769/143‑1 (MQ=11)          aaacaGGGTGCATGAAATAGCGAGGGGAGTATACTGCGCATACCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGt                                                                                                                    <  7:724085/151‑1 (MQ=11)          aaaaaGGTTGAATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGt                                                                                                                    <  3:951435/151‑1 (MQ=17)          aaaaaGGGTCCAGGAAAAAGCCAGGGGAGTATAGTGCGCTTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGt                                                                                                                    <  4:1512744/151‑1 (MQ=11)           aaacGGTTGAATGAAATATCGAGGGTTGTATTTTGGGAATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGTa                                                                                                                   <  8:1375483/151‑1 (MQ=11)           aaaaGGTTGCATGAAAAGGGGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTa                                                                                                                   <  8:207273/151‑1 (MQ=255)           aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTTAA‑GAAGT‑a                                                                                                                           >  4:449672/1‑142 (MQ=255)            aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTAc                                                                                                                  <  8:1043321/151‑1 (MQ=255)             ccccgtGCACGAATACCCGTGCGTAGTATACTGCGAATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTGATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGTa                                                                                                                   <  2:268008/144‑1 (MQ=1)              aGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACgg                                                                                                                <  2:602869/151‑1 (MQ=255)              aGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACgg                                                                                                                <  7:539604/151‑1 (MQ=255)               gggtGCATGAAAACGCGAGCGGAGTATACTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGGCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGt                                                                                                               <  3:1479771/148‑1 (MQ=255)               ggTTGCATGAAAACGCTAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGTACGGt                                                                                                               <  6:210649/151‑1 (MQ=17)               ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGt                                                                                                               <  6:517002/151‑1 (MQ=255)               ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGt                                                                                                               <  3:786402/151‑1 (MQ=255)               ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAAAATTTCTG‑CAGGA‑AAAGCAGTACGGt                                                                                                               >  6:1164935/1‑151 (MQ=18)               cggtgaaTGAAAAATCGATCGGAGTATACTGCGAATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGTACGGt                                                                                                               <  5:782522/145‑1 (MQ=11)                gTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTg                                                                                                              <  3:552782/151‑1 (MQ=255)                  tGCATGAAAAATCGAGTGGAGTATATTGCGCTTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGTACGGTGa                                                                                                             <  7:1169462/150‑1 (MQ=11)                   gcaggaaaaacctatcGGAGTATAGTGCGCATCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGTACGGTGAcg                                                                                                           <  5:1539119/136‑1 (MQ=11)                   gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCTGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGAcg                                                                                                           <  2:1402503/151‑1 (MQ=255)                   gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGAcg                                                                                                           <  5:457829/151‑1 (MQ=255)                   gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGAcg                                                                                                           <  1:334584/151‑1 (MQ=255)                   gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGAc                                                                                                            <  3:1358766/150‑1 (MQ=255)                   gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGAc                                                                                                            <  1:199667/150‑1 (MQ=255)                   gCATGAAAACCCTATTGTAGTATTCTGGGATTCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGTACGGTGAc                                                                                                            <  8:131490/150‑1 (MQ=11)                    caggAAAACCCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGAGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGTACGGTGAcgc                                                                                                          <  2:497021/148‑1 (MQ=14)                    cATGAAAACGCGAGCGGAGTATAGTGCTCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGAcgc                                                                                                          <  4:1553223/151‑1 (MQ=255)                    cATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGAcgc                                                                                                          <  1:634435/151‑1 (MQ=255)                    aaggAAAACTAGAGCGGTGTATAGTGCGCTTCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGTACGGTGAcgc                                                                                                          <  1:35907/148‑1 (MQ=11)                     aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGAcgcg                                                                                                         <  2:567333/151‑1 (MQ=255)                       gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCggg                                                                                                       <  8:317009/151‑1 (MQ=255)                       gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCggg                                                                                                       <  4:917607/151‑1 (MQ=255)                       gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCggg                                                                                                       <  3:630349/151‑1 (MQ=255)                       caattcgcgAGCGTAGAATAGACCCCCTCCCCGGACGCGGGGTGGAGCAGCCTGGGAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGTACGGTGACGCggg                                                                                                       <  6:1278855/146‑1 (MQ=11)                        aaaaatagAGCGGAGTATAGACCCACTCCGCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGTACGGTGACGCgggg                                                                                                      <  6:87182/144‑1 (MQ=11)                        aaaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCgggg                                                                                                      <  4:746572/151‑1 (MQ=255)                        aaaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGGCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCgggg                                                                                                      <  2:222926/151‑1 (MQ=255)                         aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGGt                                                                                                     <  6:1602393/151‑1 (MQ=255)                         aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGGt                                                                                                     <  6:1448963/151‑1 (MQ=255)                         aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGGCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGGt                                                                                                     <  6:1537633/151‑1 (MQ=255)                            cgcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAAATTTTTAT‑GAAGA‑AAAGCAGTACGGTGAcgagcgcggc                                                                                                   >  5:202606/1‑142 (MQ=17)                            atcgttcGGTGTATTTTGGGCATCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGTACGGTGAcg                                                                                                           <  1:202955/136‑1 (MQ=11)                            agagAGGGGAGTATACTGCCATTCCGCGGACGCGGGGTGGAGGAGCCTGGTCGCTCGTCGGTCTCATAACCCTAAGGTCGTCTGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAACTAAAAGCAGTACGGTGACGCGGGGTGGa                                                                                                  <  6:1593014/148‑1 (MQ=11)                             gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGGTGGag                                                                                                 <  4:131023/151‑1 (MQ=255)                             gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGGTGGag                                                                                                 <  8:1113360/151‑1 (MQ=255)                              cgAGGGGAGTATAGTGCGCATCGAGGGACGCGGGTTGGAGCAGCCTGGTCGGTCGTCGGGATCATAACCCGAAGGTCGTCGGTTCAAATCCGTCCCCCGCAACGAATTAAATTTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGTTGGagc                                                                                                <  8:702429/151‑1 (MQ=17)                              cgAGCGGAGTATAGTGAGAATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AACGCAGTACGGTGACGCGGGGTGGagc                                                                                                <  8:131221/151‑1 (MQ=255)                               gAGCGGAGTTTAGTGCGCATGCACGGTCGCGGGGTGGAGGAGTGTGGTTGCTGGTGGGGCTCATATCCGGAAGGTCGTCGGTGGAAATCCGGCGCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGGTGGagca                                                                                               <  7:1033835/151‑1 (MQ=17)                                aGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGGTGGagcag                                                                                              <  2:905082/151‑1 (MQ=255)                                   gttgTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGTACGGTGAcgc                                                                                                          <  8:614850/133‑1 (MQ=11)                                   ggAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAAAA‑TTGATGGAATA‑AAAGCAGTACGGTGAcgcagagcggaagagccc                                                                                           >  2:615157/1‑136 (MQ=14)                                     aGTATAGTGCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAT‑GAAGTAAAAGCAGTACGGTGAcgcg                                                                                                         <  7:440138/135‑1 (MQ=11)                                     aGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAAATTTTGAG‑GAGGT‑AAAGCAGTACGGTGACGCgggctcggaagcgcacc                                                                                         >  6:437963/1‑137 (MQ=17)                                          aGTGCGCATGCACTGACGCGTGGTGGAGCAGCCTGGTAGCTGGTCGGGCTAATAACCCGAAGGTCGTCGGTTCAAATCCTGCCCCCGCAACGAATTATAATTTGAT‑GAAGT‑AAAGCAGTACGGGGACGCGGGGTGGAGCAGCCTGGTAGCt                                                                                    <  5:1151634/151‑1 (MQ=14)                                          aGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCt                                                                                    <  3:532092/151‑1 (MQ=255)                                           gTGCTCATCCAGGGACGCGGCGGGGAGCTGGCTCGTAGCTCGGCTGTATCATAACCCCAAGGCCGGCGCGGCAAATAAGCCCCCCGCAACCAAATAAAAATTGCT‑GAAGA‑AACCCACTACCCGGACGCGGGGTGGAGCAGCCTGGTAGCtc                                                                                   <  2:482301/151‑1 (MQ=9)                                                   ccACGGACGCGGGGTGGAGCAGGGTGGTAGCTCGTCGGTCTCATAACCCGAAGGTCGTCGGTGCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGc                                                                            <  4:260103/150‑1 (MQ=255)                                                     cctgactaggtctGGAGGAGCCTTTTAGCTCGTCTGTATATTAAGCCTAATGTCGGCTGTTCAAATAATCCCCCCTCAACCAATTAAAAATTGAA‑GAAGA‑AACGCACTACCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCa                                                                         <  3:228636/139‑1 (MQ=2)                                                                       agcCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTa                                                                                       <  7:1455868/119‑1 (MQ=255)                                                                       agcCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTa                                                                                       >  6:1449473/1‑119 (MQ=255)                                                                                   tttctgcataagaagccGAATGGCGTCGGGTGAAATCCGGCCCCCTCAAGCAATAAAACTTAGAA‑GTATT‑AACCCACTACTGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGATGGTCGTCGGTTCAAATCCgg                                            <  6:1144261/135‑1 (MQ=2)                                                                                                         tcgtcgGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGGTg                                                                                                    >  8:1364390/1‑72 (MQ=255)                                                                                                                ttCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAAcc                                                                    >  5:1036241/1‑97 (MQ=255)                                                                                                                            ccccTCAATCAATTATACTTTGAT‑GAATT‑AATGGAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACgg  <  5:84834/151‑1 (MQ=18)                                                                                                                                             aTTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGt      <  6:1313960/130‑1 (MQ=255)                                                                                                                                             aTTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCACCGCAACCAATCAAATTTGATGAAGTAAAAGCAGt      >  7:1320650/1‑130 (MQ=255)                                                                                                                                                          |                                                                                                                           CGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAT‑GAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGG  >  NC_000913/2947330‑2947603

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