Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,791,128	A→T	6.0%	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	6.0%	208.7 / 2.5	67	I197N (ATC→AAC) ‡	hemE	uroporphyrinogen decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (30/33);  new base T (2/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

GATCACCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAGGGTC  >  NC_002947/5790980‑5791271                                                                                                                                                     |                                                                                                                                                gATCACCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGa                                                                                                                                                 >  1:149786/1‑149 (MQ=255)     aCCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTg                                                                                                                                             <  7:26880/149‑1 (MQ=255)       cGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGc                                                                                                                                           <  4:312069/149‑1 (MQ=255)         gCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGATCCTGGCCg                                                                                                                                         <  8:192855/149‑1 (MQ=255)                gcgcCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAgg                                                                                                                                   <  7:104657/148‑1 (MQ=255)                 cgcCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGt                                                                                                                                  <  7:348617/148‑1 (MQ=255)                     ccTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTg                                                                                                                             <  2:181039/149‑1 (MQ=255)                         gTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGTCCAGGGATAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCTCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTg                                                                                                                          <  8:39622/148‑1 (MQ=255)                         gTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGt                                                                                                                           <  3:58232/147‑1 (MQ=255)                         gTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCAGCCAGGATCTGGCCGTTAAGGTAGCTGGt                                                                                                                           <  5:67177/147‑1 (MQ=255)                          tGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGAc                                                                                                                        <  3:117323/149‑1 (MQ=255)                             tCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCAGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCAGGCCGTTTAGGTAGCGGGTGACCGa                                                                                                                     >  8:276075/1‑149 (MQ=255)                              cgcgGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGAc                                                                                                                    >  7:348295/1‑149 (MQ=255)                              cgcgGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGATCCTGGCCgg                                                                                                                                        <  2:23661/129‑2 (MQ=255)                              cgcgGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGATCCTGGCCgg                                                                                                                                        >  1:23661/1‑128 (MQ=255)                               gcgGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGa                                                                                                                     <  8:167045/147‑1 (MQ=255)                               gcgGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGAc                                                                                                                    <  5:146847/148‑1 (MQ=255)                                      aGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCGTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAg                                                                                                            <  6:303028/149‑1 (MQ=255)                                          cGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTg                                                                                                        <  8:184306/149‑1 (MQ=255)                                          cGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTg                                                                                                        >  8:142591/1‑149 (MQ=255)                                           gCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGt                                                                                                       >  8:167795/1‑149 (MQ=255)                                            cTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTgg                                                                                                                 >  4:129223/1‑138 (MQ=255)                                            cTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTgg                                                                                                                 <  3:129223/138‑1 (MQ=255)                                             tGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTgg                                                                                                                            >  7:101635/1‑126 (MQ=255)                                             tGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTgg                                                                                                                            <  8:101635/126‑1 (MQ=255)                                                 gATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCcag                                                                                                   >  6:182402/1‑147 (MQ=255)                                                       gcgcATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCcagcagcag                                                                                             <  8:324665/147‑1 (MQ=255)                                                       gcgcATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAgg                                                                                            <  8:348295/148‑1 (MQ=255)                                                        cgcATGGAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGTATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTg                                                                                          <  8:120442/149‑1 (MQ=255)                                                             gTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCca                                                                                                    <  5:139962/134‑1 (MQ=255)                                                             gTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCca                                                                                                    >  6:139962/1‑134 (MQ=255)                                                                  ccAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGc                                                                                >  6:144446/1‑149 (MQ=255)                                                                    aGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCgg                                                                              >  6:384678/1‑149 (MQ=255)                                                                      ggAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCtt                                                                                  <  8:198693/143‑1 (MQ=255)                                                                      ggAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCtt                                                                                  >  7:198693/1‑143 (MQ=255)                                                                       gagaACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGtt                                                                           <  2:176282/149‑1 (MQ=255)                                                                             tCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGAGCGACTGGGCCAGCTTGTGCTGGTGGTGGTGGTGGGGTTGGGGGTTGTGGtt                                                                     >  2:149816/1‑148 (MQ=255)                                                                             tCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCggg                                                                             <  3:258587/141‑1 (MQ=255)                                                                             tCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCggg                                                                             >  4:258587/1‑141 (MQ=255)                                                                               cTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAgg                                                                   >  4:177549/1‑149 (MQ=255)                                                                               cTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCAGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGATAAGGTAGCTGGTGACCGACTCGGCCAGCTTGTCCAGCAGCCGGTGCACGGCTTGCGGGTTGTCGTAgg                                                                   >  7:62898/1‑149 (MQ=255)                                                                                tGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGc                                                                  >  3:123879/1‑149 (MQ=255)                                                                                     ggCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATgg                                                             <  7:167795/149‑1 (MQ=255)                                                                                       ccgccgccgACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGcc                                                           >  5:100256/1‑149 (MQ=255)                                                                                         gccgccgACAGGTTGCCGCTCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCa                                                                <  1:146043/142‑1 (MQ=255)                                                                                         gccgccgACAGGTTGCCGCTCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCa                                                                >  2:146043/1‑142 (MQ=255)                                                                                          ccgccgACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTg                                                        <  7:211112/149‑1 (MQ=255)                                                                                                aCAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTg                                                  >  1:338496/1‑149 (MQ=255)                                                                                                 cAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGc                                                 <  7:132911/149‑1 (MQ=255)                                                                                                    gTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCgg                                               >  5:105545/1‑148 (MQ=255)                                                                                                     ttGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGGTTGTCCAGCAGCAGGTGCAGGGCTTTCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGaa                                             <  2:340315/149‑1 (MQ=255)                                                                                                          gccCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCggg                                                                             >  3:245491/1‑112 (MQ=255)                                                                                                          gccCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCggg                                                                             <  4:245491/112‑1 (MQ=255)                                                                                                            cccTGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTACGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGTTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTc                                      <  3:222954/149‑1 (MQ=255)                                                                                                            cccAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTc                                      <  7:276075/149‑1 (MQ=255)                                                                                                              cAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTcga                                    <  4:123879/149‑1 (MQ=255)                                                                                                                      gAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCCGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTAGTAGGCCATGGCCTTGGTCTTGAGGAAGTCCTTCGACGAgccgc                            >  4:49448/1‑149 (MQ=255)                                                                                                                      gAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAgccgc                            >  5:100162/1‑149 (MQ=255)                                                                                                                      gAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAgccgc                            >  7:13809/1‑149 (MQ=255)                                                                                                                            cTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGa                      >  3:20795/1‑149 (MQ=255)                                                                                                                                   gCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgta               <  5:144446/149‑1 (MQ=255)                                                                                                                                    ccTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgtag              >  6:356398/1‑149 (MQ=255)                                                                                                                                    ccTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgtag              <  6:105545/149‑1 (MQ=255)                                                                                                                                      tGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgtagg             >  2:359629/1‑148 (MQ=255)                                                                                                                                            acGGCCAGATCCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAg      >  7:105568/2‑149 (MQ=255)                                                                                                                                              ggCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTAGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAggg    >  2:145944/1‑149 (MQ=255)                                                                                                                                                 cAGGATCTGGCCGTTAAGGTAGCTGGTGCCCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAGGGTc  >  8:185038/1‑148 (MQ=255)                                                                                                                                                    gATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTcg                                     >  6:321007/1‑110 (MQ=255)                                                                                                                                                     |                                                                                                                                                GATCACCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAGGGTC  >  NC_002947/5790980‑5791271

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