Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	350,615	G→A	5.4%	G10S (GGC→AGC)	hisH →	imidazole glycerol phosphate synthase subunit HisH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	350,615	0	G	A	5.4%	225.8 / 2.0	74	G10S (GGC→AGC)	hisH	imidazole glycerol phosphate synthase subunit HisH
Reads supporting (aligned to +/- strand):  ref base G (35/35);  new base A (2/2);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACAACACCCACCACCAGATCGAAACCGTATT‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGACCGTGTGGTGTTCCCGGGCGTGGGCGCGATTCGCGACTGCATGGCCGAGA  >  NC_002947/350467‑350762                                                                                                                                                        |                                                                                                                                                    ccACAACACCCACCACCAGATCGAAACCGTATT‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGCGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATg                                                                                                                                                     >  4:352626/1‑149 (MQ=255) ccACAACACCCACCACCAGATCGAAACCGTATT‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATg                                                                                                                                                     <  7:35196/149‑1 (MQ=255)  cacaACACCCACCACCAGATCGAAACCGTATT‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATgg                                                                                                                                                    >  5:423599/1‑149 (MQ=255)     aacaCCCACCACCAGATCGAAACCGTATT‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCAt                                                                                                                                                 <  3:158453/149‑1 (MQ=255)     aacaCCCACCACCAGATCGAAACCGTATT‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCAt                                                                                                                                                 >  3:158904/1‑149 (MQ=255)       cacccaccACCAGATCGAAACCGTATT‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATgg                                                                                                                                               <  2:314393/149‑1 (MQ=255)        acccaccACCAGATCGAAACCGTATT‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATgtg                                                                                                                                              >  2:235938/1‑147 (MQ=255)          ccaccaccaGATCGAAACCGTATTACAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATggg                                                                                                                                              >  7:257248/1‑148 (MQ=255)           caccaccaGATCGAAACCGTATT‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTaa                                                                                                                                           >  3:194796/1‑149 (MQ=255)            accaccaGATCGAAACCGTATT‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGATTTCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTaa                                                                                                                                           <  1:121576/149‑1 (MQ=255)            accaccaGATCGAAACCGTATT‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAAc                                                                                                                                          >  5:335228/1‑149 (MQ=255)              caccaGATCGAAACCGTATT‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCt                                                                                                                                        <  3:75292/149‑1 (MQ=255)              caccaGATCGAAACCGTATT‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCt                                                                                                                                        >  4:324798/1‑149 (MQ=255)                         aaCCGTATT‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTgg                                                                                                                             >  5:90891/1‑149 (MQ=255)                            cGTATT‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCAGTAATCG‑ACTATAGCATGGGTAACCTGCACTCGGTGGcc                                                                                                                           >  4:118275/1‑148 (MQ=255)                            cGTATT‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCAGTAATCG‑ACTATAGCATGGGTAACCTGCACTCGGTGGcc                                                                                                                           >  4:123989/1‑148 (MQ=255)                                 t‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGgc                                                                                                                      <  3:324798/148‑1 (MQ=255)                                   cAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGgcgc                                                                                                                    >  1:165731/1‑149 (MQ=255)                                    aaGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGgcgc                                                                                                                    <  4:96483/148‑1 (MQ=255)                                    aaGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCt                                                                                                                   >  2:246487/1‑149 (MQ=255)                                    aaGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCt                                                                                                                   >  3:401560/1‑149 (MQ=255)                                    aaGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCt                                                                                                                   >  7:356241/1‑149 (MQ=255)                                     aGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTGGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCGCACTATGGCATGGGTAACCTGCACTCGGTGGCCAAg                                                                                                                        <  1:60489/144‑1 (MQ=255)                                     aGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAg                                                                                                                        >  2:60489/1‑143 (MQ=255)                                      ggCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTg                                                                                                                 <  1:246487/149‑1 (MQ=255)                                       gCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCtt                                                                                                                  >  5:20095/1‑147 (MQ=255)                                       gCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGa                                                                                                                <  6:423599/149‑1 (MQ=255)                                         tttCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTg                                                                                                                 <  6:159661/146‑1 (MQ=255)                                         tttCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTg                                                                                                                 >  5:159661/1‑146 (MQ=255)                                           tCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGc                                                                                                              <  2:165731/147‑1 (MQ=255)                                           tCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGc                                                                                                              <  6:229668/147‑1 (MQ=255)                                               ccGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAg                                                                                                        >  4:222842/1‑149 (MQ=255)                                                cgcgCTCTGCGCATGGCCATCAACCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTACACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGTCATGGGTAACCTCCACTCGGTGGCAAATGCGCAAGAGCACGAAgg                                                                                                       >  6:292107/1‑149 (MQ=255)                                                          gcATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTg                                                                                                                              <  3:205392/116‑1 (MQ=255)                                                          gcATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTg                                                                                                                              >  4:205392/1‑116 (MQ=255)                                                          gcATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTa                                                                                               >  8:131295/1‑147 (MQ=255)                                                               gCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGc                                                                                         <  4:401560/148‑1 (MQ=255)                                                                ccATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTg                                                                                       <  4:194796/149‑1 (MQ=255)                                                                       ccTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGa                                                                                                                >  4:200295/1‑117 (MQ=255)                                                                       ccTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGc                                                                                                              >  8:152749/1‑119 (MQ=255)                                                                       ccTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGc                                                                                                              <  7:152749/119‑1 (MQ=255)                                                                       ccTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTa                                                                                                         <  3:260871/124‑1 (MQ=255)                                                                       ccTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTa                                                                                                         >  4:260871/1‑124 (MQ=255)                                                                       ccTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGt                                                                                                >  3:403783/1‑133 (MQ=255)                                                                       ccTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGt                                                                                                <  4:403783/133‑1 (MQ=255)                                                                       ccTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTAc                                                                                  >  5:369252/1‑147 (MQ=255)                                                                           gacgaGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCa                                                                                <  7:298440/145‑1 (MQ=255)                                                                           gacgaGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCa                                                                                >  8:298440/1‑145 (MQ=255)                                                                              gaGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGc                                                                         <  6:90891/149‑1 (MQ=255)                                                                                gCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCt                                                                       >  1:206329/1‑149 (MQ=255)                                                                                 cGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTc                                                                      <  6:419553/149‑1 (MQ=255)                                                                                 cGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTc                                                                      <  2:206329/149‑1 (MQ=255)                                                                                  gTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCa                                                                     <  6:20095/149‑1 (MQ=255)                                                                                      ggCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGa                                                                 >  3:174676/1‑149 (MQ=255)                                                                                       gCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGa                                                                 <  1:354731/148‑1 (MQ=255)                                                                                             cAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGa                                                          <  8:356241/149‑1 (MQ=255)                                                                                             cAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGa                                                          <  1:378469/149‑1 (MQ=255)                                                                                              aGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAg                                                         >  2:385143/1‑149 (MQ=255)                                                                                                     cTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGAc                                                   >  1:179384/1‑149 (MQ=255)                                                                                                     cTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGAc                                                   >  1:185500/1‑149 (MQ=255)                                                                                                          accaAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGACCgtgt                                              >  4:154995/1‑149 (MQ=255)                                                                                                           ccaAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGACCGtgtg                                             <  8:257248/149‑1 (MQ=255)                                                                                                           ccaAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGACCGtgtg                                             <  3:352626/149‑1 (MQ=255)                                                                                                           ccaAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGACCGtgtg                                             <  7:131295/149‑1 (MQ=255)                                                                                                                                   cGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTgg                                                                                      <  7:242569/84‑1 (MQ=255)                                                                                                                                   cGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTgg                                                                                      >  8:242569/1‑84 (MQ=255)                                                                                                                                          cGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGACCGTGTGGTGTTCCCGGGCGTGGGCGCGATTCGCGACt              >  3:119917/1‑149 (MQ=255)                                                                                                                                          cGTAATCC‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGACCGTGTGGTGTTCCCGGGCGTGGGCGCGATTCGCGACt              >  4:6620/1‑149 (MQ=255)                                                                                                                                            tAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGACCGTGTGGTGTTCCCGGGCGTGGGCGCGATTCGCGACTGc            <  4:158904/149‑1 (MQ=255)                                                                                                                                             aaTCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGACCGTGTGGTGTTCCCGGGCGTGGGCGCGATTc                    >  2:327394/1‑140 (MQ=255)                                                                                                                                             aaTCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGACCGTGTGGTGTTCCCGGGCGTGGGCGCGATTc                    <  1:327394/140‑1 (MQ=255)                                                                                                                                             aaTCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGACCGTGTGGTGTTCCCGGGCGTGGGCGCGATTCGCGACTGCa           >  1:39898/1‑149 (MQ=255)                                                                                                                                                 g‑ACTATAGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGACCGTGTGGTGTTCCCGGGCGTGGGCGCGATTCGCGACTGCATgg        <  3:118275/148‑1 (MQ=255)                                                                                                                                                 g‑ACTATAGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGACCGTGTGGTGTTCCCGGGCGTGGGCGCGATTCGCGACTGCATgg        <  3:123989/148‑1 (MQ=255)                                                                                                                                                   aCTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGACCGTGTGGTGTTCCCGGGCGTGGGCGCGATTCGCGACTGCATGGcc      <  6:369252/149‑1 (MQ=255)                                                                                                                                                       tGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGACCGTGTGGTGTTCCCGGGCGTGGGCGCGATTCGCGACTGCATGGCCgaga  <  3:154995/149‑1 (MQ=255)                                                                                                                                                       tGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGACCGTGTGGTGTTCCCGGGCGTGGGCGCGATTCGCGACTGCATGGCCgaga  <  3:6620/149‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                    CCACAACACCCACCACCAGATCGAAACCGTATT‑CAAGGCTTTCGGCCGCGCTCTGCGCATGGCCATCACCCTCGACGAGCGTATGGCCGGGCAGA‑TGCCTTCCACCAAAGGGTGCCTGTAAATGCAAACGGTAGCCGTAATCG‑ACTATGGCATGGGTAACCTGCACTCGGTGGCCAAGGCGCTTGAGCACGTAGGTGCCGGTAAAGTGCTGGTTACCAGCGATGCCTCAGTGATCCGCGAGGCCGACCGTGTGGTGTTCCCGGGCGTGGGCGCGATTCGCGACTGCATGGCCGAGA  >  NC_002947/350467‑350762

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