Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,038,228	G→A	8.8%	A389A (GCC→GCT)	PP_0897 ←	class 1 fumarate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,038,228	0	G	A	8.8%	188.1 / 7.0	68	A389A (GCC→GCT)	PP_0897	class 1 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base A (4/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

CAACGGCCATCAGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCGGTCAGCATCTTG  >  NC_002947/1038080‑1038375                                                                                                                                                     |                                                                                                                                                      cAACGGCCATCAGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCg                                                                                                                                                       >  4:10891/1‑149 (MQ=255)  aaCGGCCATCAGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCgg                                                                                                                                                      >  3:147126/1‑149 (MQ=255)          tCAGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTAGGCCCAGCCGGcccg                                                                                                                                              >  7:345298/1‑148 (MQ=255)            aGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TAc                                                                                                                                          <  3:10891/149‑1 (MQ=255)              gTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTAGGCCCAGCCGGCCC‑‑GACCa                                                                                                                                        <  8:345298/149‑1 (MQ=255)                 cacGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACtt                                                                                                                                     >  1:120666/1‑149 (MQ=255)                 cacGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATAAGGGTAGCAGTTGTCGGGCCGGCTGGGCC‑‑TAACACt                                                                                                                                      >  6:91758/1‑148 (MQ=255)                             gtCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTcacccacc                                                                                                                          <  3:185704/148‑1 (MQ=255)                             gtCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATATGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCg                                                                                                                         >  2:283383/1‑149 (MQ=255)                                       gCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCgg                                                                                                               >  2:368604/1‑149 (MQ=255)                                       gCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCgg                                                                                                               >  6:205932/1‑149 (MQ=255)                                        ccTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCggg                                                                                                              >  1:90764/1‑149 (MQ=255)                                             aTGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGAc                                                                                                         <  5:361445/149‑1 (MQ=255)                                             aTGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGAc                                                                                                         <  7:300336/149‑1 (MQ=255)                                              tGGCGGTCGGGCCACTCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGTCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGAc                                                                                                         <  8:331768/148‑1 (MQ=255)                                              tGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACg                                                                                                        <  1:368604/149‑1 (MQ=255)                                                 cGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACgtag                                                                                                     <  4:173509/149‑1 (MQ=255)                                                  ggTCGGGCCACGCTCGGACTTGCCAATCTTGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACgtag                                                                                                     >  1:207705/1‑148 (MQ=255)                                                          cACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAg                                                                                            <  2:131443/149‑1 (MQ=255)                                                          cACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAg                                                                                            >  4:99757/1‑149 (MQ=255)                                                             gCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCgg                                                                                         >  1:148387/1‑149 (MQ=255)                                                                  gACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCttt                                                                                    >  8:151888/1‑149 (MQ=255)                                                                  gACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCttt                                                                                    >  6:273731/1‑149 (MQ=255)                                                                     ttGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAg                                                                                 >  6:133015/1‑149 (MQ=255)                                                                        ccAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTc                                                                              >  2:314421/1‑149 (MQ=255)                                                                        ccAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTc                                                                              >  4:314657/1‑149 (MQ=255)                                                                          aaTCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGtagtag                                                                                                  <  1:29821/127‑1 (MQ=255)                                                                          aaTCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGtagtag                                                                                                  >  2:29821/1‑127 (MQ=255)                                                                           aTCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTAc                                                                           <  1:283383/149‑1 (MQ=255)                                                                              aTGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCg                                                                         <  3:314657/148‑1 (MQ=255)                                                                                  ccAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAAc                                                                    >  8:231774/1‑149 (MQ=255)                                                                                   cagcagGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAAc                                                                    >  5:8970/1‑148 (MQ=255)                                                                                    agcagGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAActc                                                                  <  7:151888/149‑1 (MQ=255)                                                                                    agcagGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAActc                                                                  <  5:273731/149‑1 (MQ=255)                                                                                       agGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTc                                                               <  3:99757/149‑1 (MQ=255)                                                                                       agGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTc                                                               <  2:120666/149‑1 (MQ=255)                                                                                        ggCCAGTCTGCTCAAGGAGCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCg                                                              <  6:251291/149‑1 (MQ=255)                                                                                         gCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGc                                                                                           >  4:183064/1‑119 (MQ=255)                                                                                         gCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGc                                                                                           <  3:183064/119‑1 (MQ=255)                                                                                              tCTGCTCAAGGATCTGGCGAGTGAACTTGTCGATACGGGTAGCAGTTTTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCTGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACgg                                                        <  2:90764/149‑1 (MQ=255)                                                                                              tCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACgg                                                        <  3:78594/149‑1 (MQ=255)                                                                                                  cTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTc                                                               >  7:3711/1‑138 (MQ=255)                                                                                                  cTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTc                                                               <  8:3711/138‑1 (MQ=255)                                                                                                  cTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGtt                                                      <  5:91758/147‑1 (MQ=255)                                                                                                  cTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGa                                                    <  4:356072/149‑1 (MQ=255)                                                                                                  cTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGa                                                    <  3:126511/149‑1 (MQ=255)                                                                                                   tCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCGATACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGtt                                                      >  1:290109/1‑148 (MQ=255)                                                                                                   tCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAg                                                   >  1:20502/1‑149 (MQ=255)                                                                                                   tCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAg                                                   >  5:338952/1‑149 (MQ=255)                                                                                                   tCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAg                                                   >  4:170094/1‑149 (MQ=255)                                                                                                           cTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCa                                            <  1:314421/148‑1 (MQ=255)                                                                                                              gCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCgg                                                                                         <  6:216364/100‑1 (MQ=255)                                                                                                              gCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCgg                                                                                         >  5:216364/1‑100 (MQ=255)                                                                                                              gCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTc                                                                                   <  6:282949/106‑1 (MQ=255)                                                                                                              gCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCaa                                           <  5:205932/146‑1 (MQ=255)                                                                                                              gCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCa                                        >  6:354366/1‑149 (MQ=255)                                                                                                                gAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCttt                                                                                    <  2:114493/103‑1 (MQ=255)                                                                                                                gAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCAGGCCGACGTAGTAGATGAAGCGGCCttt                                                                                    >  1:114493/1‑103 (MQ=255)                                                                                                                   tGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATgcgc                                   >  8:57842/1‑149 (MQ=255)                                                                                                                      aCTTGTCCATACGGGTAGCAGTGGGCGGGCCGGCGGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTtg                                <  5:302890/149‑1 (MQ=255)                                                                                                                                 cGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATc                     <  6:338952/149‑1 (MQ=255)                                                                                                                                     acgaaGTGGTAGGCCCAGCCGGCCC‑‑GACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTgcgc                           >  6:99245/5‑139 (MQ=255)                                                                                                                                     acgaaGTGGTAGGCCCAGCCGGCCC‑‑GACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTgcgc                           <  5:99245/135‑1 (MQ=255)                                                                                                                                     acgaaGTGGTAGGCCCAGCCGGCCC‑‑GACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACAATGCGATTGTGCTCGGCGTCACgg                 >  6:256416/5‑149 (MQ=255)                                                                                                                                     acgaaGTGGTAGGCCCAGCCGGCCC‑‑GACACCTTAGTCACCCACAGGGTCGACCGGGCCGACGTAGTAGATGTCGCGGGCTTTCAGGTCTACCGGGCAATCTTCGCAACGGTTGAGAATTTCACAAATACGAATGTGCGAGTCATCACg                  >  6:256418/5‑148 (MQ=255)                                                                                                                                        cAGTGGTCGGGCCGGCTGGGCC‑‑TACCCATTCGTAACACACCGTGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCg              >  7:37928/1‑149 (MQ=255)                                                                                                                                        cAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCg              <  2:20502/149‑1 (MQ=255)                                                                                                                                           tGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCGGTc           <  4:147126/149‑1 (MQ=255)                                                                                                                                                 ggCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCGGTCAGCATc     <  3:170094/149‑1 (MQ=255)                                                                                                                                                    cGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCGGTCAGCATCTTg  <  5:354366/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CAACGGCCATCAGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCGGTCAGCATCTTG  >  NC_002947/1038080‑1038375

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