Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,591,288	T→A	6.2%	I120N (ATC→AAC)	liuC →	methylglutaconyl‑CoA hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,591,288	0	T	A	6.2%	194.8 / 2.5	64	I120N (ATC→AAC)	liuC	methylglutaconyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base A (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

AGCAGTCGGCGCAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTACGCCCTCACCGCCG  >  NC_002947/4591147‑4591434                                                                                                                                                |                                                                                                                                                    agcagTCGGCGCAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAgctgc                                                                                                                                              >  4:293875/1‑149 (MQ=255)          cgcAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATg                                                                                                                                     <  5:51482/149‑1 (MQ=255)           gcAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATgg                                                                                                                                    >  1:272842/1‑149 (MQ=255)             aaCTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCTACATgtc                                                                                                                                   >  6:303757/1‑146 (MQ=255)                     ttCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGNATCGCCTCAAGGNCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCANCCTTTGGCGGCGCGCTGGGNCTGAACAGCTGCTGCGACATGGCCATCGgcgc                                                                                                                          >  5:221547/1‑149 (MQ=255)                     ttCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGgcgc                                                                                                                          >  6:305226/1‑149 (MQ=255)                     ttCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGgcgc                                                                                                                          >  5:194985/1‑149 (MQ=255)                     ttCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGAACAGCTGCTGCGACATGGCCATCGgcgc                                                                                                                          >  5:216452/1‑149 (MQ=255)                     ttCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGA‑CCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGcc                                                                                                                         >  6:305506/1‑149 (MQ=255)                     ttCAACACCAACCCGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGgcgc                                                                                                                          >  3:85214/1‑149 (MQ=255)                           acCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGa                                                                                                                    <  1:134307/149‑1 (MQ=255)                           acCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGa                                                                                                                    >  5:56042/1‑149 (MQ=255)                              aaCCTGGATGATGCCCACGAACGGGGCGAGCTGTTGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGc                                                                                                                 <  5:582/149‑1 (MQ=255)                               aCCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGcc                                                                                                                <  3:205237/149‑1 (MQ=255)                                ccTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGcc                                                                                                                >  4:205237/1‑148 (MQ=255)                                  tGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑aa                                                                                                            >  5:268731/1‑149 (MQ=255)                                    gatgatGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑aa                                                                                                            <  4:41793/147‑1 (MQ=255)                                    gatgatGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACt                                                                                                          >  2:193856/1‑149 (MQ=255)                                    gatgatGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACt                                                                                                          >  8:137198/1‑149 (MQ=255)                                       gatgCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACtgtg                                                                                                       >  5:75249/1‑149 (MQ=255)                                        atgCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGc                                                                                                      <  7:137198/149‑1 (MQ=255)                                         tgCCCACGATNTGGGCTNGCTGATGTACGCCCTGCATCNCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGAACAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGcc                                                                                                     <  6:221547/149‑1 (MQ=255)                                         tgCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGAACAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGcc                                                                                                     <  6:216452/149‑1 (MQ=255)                                            ccACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGt                                                                                                  >  4:345311/1‑149 (MQ=255)                                             cACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGCTGAAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGATGCGACATGGCCATCGGCGCCGAAGCCGCCC‑AACTGTGCCTg                                                                                                   >  1:107073/1‑149 (MQ=255)                                             cACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGTCTCCTAAGACGCCC‑AACtgtagatcg                                                                                                  >  8:164070/1‑142 (MQ=255)                                             cACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACtg                                                                                                         <  7:164070/141‑1 (MQ=255)                                                 aaCTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGaa                                                                                             <  3:293875/149‑1 (MQ=255)                                                 aaCTGGGCGAGCTGATGTACGCCCTGCATCACCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGaa                                                                                             <  2:181314/149‑1 (MQ=255)                                                  ccggggCGAGCTGATGTCCGCCCTGCATCGCCTCAAGGCCCCGACCCTCGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAg                                                                                            <  8:166682/146‑1 (MQ=255)                                                      ggCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTAc                                                                                         <  5:305226/148‑1 (MQ=255)                                                       gCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGc                                                                                       <  6:56042/149‑1 (MQ=255)                                                              gATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGcc                                                                                >  1:25095/1‑149 (MQ=255)                                                                 gTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGTGCCGAAGACGCCC‑AACTGTGTCTGTCGGAAGTACGCATCGGCCTgg                                                                             >  2:910/1‑149 (MQ=255)                                                                  tACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTTGc                                                                            <  6:268731/149‑1 (MQ=255)                                                                       ccTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCcgg                                                                       >  5:165080/1‑149 (MQ=255)                                                                        cTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGt                                                                                           <  3:226429/128‑1 (MQ=255)                                                                        cTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGt                                                                                           >  4:226429/1‑128 (MQ=255)                                                                           cATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCcggcgg                                                                    <  1:193856/148‑1 (MQ=255)                                                                           cATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGt                                                                   <  2:25095/149‑1 (MQ=255)                                                                             tCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGtca                                                                 >  7:6269/1‑149 (MQ=255)                                                                                 cTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAg                                                             >  5:305130/1‑149 (MQ=255)                                                                                  tCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGAATTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGc                                                            >  1:73065/1‑149 (MQ=255)                                                                                     aGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGc                                                                                                      >  4:274603/1‑104 (MQ=255)                                                                                     aGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGc                                                                                                      <  3:274603/104‑1 (MQ=255)                                                                                        ccccGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTc                                                      <  4:85214/149‑1 (MQ=255)                                                                                         cccGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGa                                                                                              >  2:159157/1‑108 (MQ=255)                                                                                         cccGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGa                                                                                              <  1:159157/108‑1 (MQ=255)                                                                                         cccGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCg                                                                                   >  8:332041/1‑119 (MQ=255)                                                                                         cccGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCg                                                                                   <  7:332041/119‑1 (MQ=255)                                                                                           cGACCCTGGCCGTGGTGC‑‑AAGGCTCCGCCTTTGGCGGCGGGCTGGGACTGATCAGCTGCTGCGACATGTCCATCTGCGCCGAAGAGGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTg                                                   <  8:168582/149‑1 (MQ=255)                                                                                              cccTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTc                                                <  6:75249/149‑1 (MQ=255)                                                                                                cTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCaa                                              <  5:305506/149‑1 (MQ=255)                                                                                                cTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCaa                                              <  6:194985/149‑1 (MQ=255)                                                                                                 tGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAg                                             >  5:174294/1‑149 (MQ=255)                                                                                                      gtggtgC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCAt                                        <  3:345311/149‑1 (MQ=255)                                                                                                         gtgC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTg                                      <  6:305130/148‑1 (MQ=255)                                                                                                                                       tGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTAc               >  5:245651/1‑143 (MQ=255)                                                                                                                                       tGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTAc               <  6:245651/143‑1 (MQ=255)                                                                                                                                       tGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTACGCCCt          >  7:150648/1‑148 (MQ=255)                                                                                                                                           ccTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCCAAACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTACGCCCt          >  1:260365/1‑145 (MQ=255)                                                                                                                                           ccTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTTGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTACGCCCt          <  2:260365/144‑1 (MQ=255)                                                                                                                                              gATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTACGCCCTCAccgccg  >  6:38716/1‑149 (MQ=255)                                                                                                                                               aTCAGCTGCTGCGACATGGCCATCGGTGCCGAAGACGCCC‑AACTGTGTCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTACGCCCTCAccgcc   <  1:910/147‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                    AGCAGTCGGCGCAACTGGACTTCAACACCAACCTGGATGATGCCCACGAACTGGGCGAGCTGATGTACGCCCTGCATCGCCTCAAGGCCCCGACCCTGGCCGTGGTGC‑‑AAGGTGCAGCCTTTGGCGGCGCGCTGGGCCTGATCAGCTGCTGCGACATGGCCATCGGCGCCGAAGACGCCC‑AACTGTGCCTGTCGGAAGTACGCATCGGCCTGGCCCCGGCGGTCATCAGCCCGTTCGTGGTCAAGGCCATTGGCGAACGCGCCGCGCGCCGCTACGCCCTCACCGCCG  >  NC_002947/4591147‑4591434

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: ‑