Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,791,128	A→T	7.6%	I197N (ATC→AAC) ‡	hemE ←	uroporphyrinogen decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,791,128	0	A	T	7.6%	193.4 / 5.6	66	I197N (ATC→AAC) ‡	hemE	uroporphyrinogen decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (27/34);  new base T (2/3);  total (29/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 = 5.72e-01

GATCACCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAGGGTCCAC  >  NC_002947/5790980‑5791274                                                                                                                                                     |                                                                                                                                                   gATCACCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGa                                                                                                                                                    <  3:254601/149‑1 (MQ=255)   tCACCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATc                                                                                                                                                  >  2:318398/1‑149 (MQ=255)     aCCGGCACCTTGCGCCCCGCGTGCGGGGCGCTCCGGCCGCCGCCGCTCCTGCGCCTGTACGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTg                                                                                                                                                <  1:39304/149‑1 (MQ=255)      ccGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTgg                                                                                                                                               <  8:33332/149‑1 (MQ=255)      ccGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTgg                                                                                                                                               >  2:54582/1‑149 (MQ=255)        ggCACCTTGCGCCCCTCGTGCTCGCGGATCAGACCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGGCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGcc                                                                                                                                             >  1:33579/1‑149 (MQ=255)         gCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCg                                                                                                                                            <  3:241965/149‑1 (MQ=255)            ccTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGtt                                                                                                                                          >  4:177293/1‑148 (MQ=255)            ccTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGATCCTGGCCggcg                                                                                                                                         <  4:171572/149‑4 (MQ=255)             cTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTaa                                                                                                                                        <  7:99625/149‑1 (MQ=255)                gcgcCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAg                                                                                                                                       <  1:318398/147‑1 (MQ=255)                gcgcCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAg                                                                                                                                       <  1:54582/147‑1 (MQ=255)                     ccTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTg                                                                                                                                >  2:63747/1‑149 (MQ=255)                      cTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTgg                                                                                                                               <  1:299195/149‑1 (MQ=255)                         gTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTg                                                                                                                             <  1:63747/148‑1 (MQ=255)                          tGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGAc                                                                                                                           <  3:177293/149‑1 (MQ=255)                             tCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGCGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTggcccagtc                                                                                                                        >  5:277356/1‑142 (MQ=255)                               gcgGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACt                                                                                                                      >  6:83530/1‑149 (MQ=255)                               gcgGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACt                                                                                                                      >  6:115735/1‑149 (MQ=255)                               gcgGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACt                                                                                                                      >  4:59642/1‑149 (MQ=255)                                 ggATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTgg                                                                                                                    <  5:250559/149‑1 (MQ=255)                                      aGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCa                                                                                                                <  3:59642/148‑1 (MQ=255)                                            cTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTc                                                                                                         >  8:237777/1‑149 (MQ=255)                                              gACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTc                                                                                                         >  3:35669/1‑147 (MQ=255)                                                       gcgcATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCATGTGTCGAAGATCTGTACGGGCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCcagcagca                                                                                                 <  4:87284/146‑1 (MQ=255)                                                       gcgcATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCcagcagcag                                                                                                <  7:237777/147‑1 (MQ=255)                                                       gcgcATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGt                                                                                              >  7:141751/1‑149 (MQ=255)                                                         gcATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTg                                                                                             >  6:3217/1‑148 (MQ=255)                                                         gcATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGc                                                                                            <  4:38363/149‑1 (MQ=255)                                                          cATGTAGGCCAGGGACAACTCGTGATCGGCCGCCGCCGAGTGGTTGGCGCCCCAGTTGTCGAAGATCTGTACGGCCTGGGCGCGGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACGGACTGGGCCAGCTTGTCCAGCAGCACGTGCa                                                                                           <  2:33579/149‑1 (MQ=255)                                                                ggCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCtt                                                                                     >  1:46068/1‑149 (MQ=255)                                                                  ccAGGGAGAACTCCTGATTGGCCGGCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTCTACGGCCTGGGCGCCGGCCTGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGc                                                                                   <  1:263716/149‑1 (MQ=255)                                                                   cAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCg                                                                                  <  2:46068/149‑1 (MQ=255)                                                                    aGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCgg                                                                                 <  7:83658/149‑1 (MQ=255)                                                                    aGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGACAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCgg                                                                                 >  8:268123/1‑149 (MQ=255)                                                                      ggAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGt                                                                               >  5:248722/1‑149 (MQ=255)                                                                           aCTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGGCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCg                                                                          <  6:97109/149‑1 (MQ=255)                                                                             tCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGc                                                                                   >  2:296681/1‑138 (MQ=255)                                                                             tCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGc                                                                                   <  1:296681/138‑1 (MQ=255)                                                                             tCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCg                                                                          >  7:29187/1‑147 (MQ=255)                                                                                tGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGc                                                                     >  4:184968/1‑149 (MQ=255)                                                                                    aGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATg                                                                 >  4:73908/1‑149 (MQ=255)                                                                                          ccgccgACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTg                                                           <  8:141751/149‑1 (MQ=255)                                                                                             ccgACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTc                                                        <  2:208974/149‑1 (MQ=255)                                                                                               gACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCtt                                                      >  2:286974/1‑149 (MQ=255)                                                                                                 cAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGc                                                    <  6:119327/149‑1 (MQ=255)                                                                                                     ttGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGaa                                                <  6:248722/149‑1 (MQ=255)                                                                                                       gccgccCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAg                                               <  3:309608/148‑1 (MQ=255)                                                                                                            cccAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTc                                         >  6:28807/1‑149 (MQ=255)                                                                                                              cAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTcga                                       <  5:3217/149‑1 (MQ=255)                                                                                                               aGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTcgac                                      <  8:29187/149‑1 (MQ=255)                                                                                                                    tCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAgc                                  <  5:83530/148‑1 (MQ=255)                                                                                                                            cTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGa                         >  1:217389/1‑149 (MQ=255)                                                                                                                               tACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGAc                        <  3:184968/147‑1 (MQ=255)                                                                                                                                   gCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgta                  <  2:217389/149‑1 (MQ=255)                                                                                                                                    ccTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgtag                 <  7:268123/149‑1 (MQ=255)                                                                                                                                    ccTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgtag                 >  5:332668/1‑149 (MQ=255)                                                                                                                                        ggCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTg                                                     <  7:148990/109‑1 (MQ=255)                                                                                                                                        ggCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTg                                                     >  8:148990/1‑109 (MQ=255)                                                                                                                                        cttaacGGCCAGATCCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTg                                                           <  6:319908/98‑1 (MQ=255)                                                                                                                                        cttaacGGCCAGATCCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTg                                                           >  5:319908/6‑103 (MQ=255)                                                                                                                                         ttaacGGCCAGATCCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGc            >  6:99042/5‑149 (MQ=255)                                                                                                                                           aacGGCCAGATCCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGcc           <  5:99042/146‑1 (MQ=255)                                                                                                                                             cGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAgg        <  4:35669/149‑1 (MQ=255)                                                                                                                                              ggCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCATGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAggg       >  8:88913/1‑149 (MQ=255)                                                                                                                                                ccAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAGGGTc     >  4:43058/1‑149 (MQ=255)                                                                                                                                                   ggATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAGGGTCCAc  <  3:43058/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GATCACCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAGGGTCCAC  >  NC_002947/5790980‑5791274

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