Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,638,813	G→C	5.9%	R270P (CGC→CCC)	PP_3205 →	fumarylacetoacetate hydrolase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,638,813	0	G	C	5.9%	217.4 / 3.0	67	R270P (CGC→CCC)	PP_3205	fumarylacetoacetate hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base G (31/32);  new base C (2/2);  total (33/34)
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.05e-01

GTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCATTGAC  >  NC_002947/3638668‑3638961                                                                                                                                                   |                                                                                                                                                      gTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAg                                                                                                                                                    >  3:77200/1‑149 (MQ=255)   gCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCg                                                                                                                                                  >  7:293173/1‑149 (MQ=255)     aCTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCACCGCCGAGCGCAGCGTg                                                                                                                                                >  5:273756/1‑149 (MQ=255)       tGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTggg                                                                                                                                              >  7:90792/1‑149 (MQ=255)         gAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTggg                                                                                                                                              <  8:293173/147‑1 (MQ=255)           aCGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCg                                                                                                                                          >  4:181536/1‑149 (MQ=255)           aCGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCg                                                                                                                                          >  2:109229/1‑149 (MQ=255)            cGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCgg                                                                                                                                         >  6:256567/1‑149 (MQ=255)             ggCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCCCACGCTGCGCTCGGc                                                                                                                                        <  6:153209/149‑1 (MQ=255)               cGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑t                                                                                                                                     <  4:284481/149‑1 (MQ=255)               cGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑t                                                                                                                                     <  7:218543/149‑1 (MQ=255)                  gTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCgatgc                                                                                                                                   <  8:95530/149‑5 (MQ=255)                     gTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGc                                                                                                                                                   <  2:35670/130‑1 (MQ=255)                     gTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTg                                                                                                                                                <  7:79451/133‑1 (MQ=255)                     gTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTg                                                                                                                                                >  8:79451/1‑133 (MQ=255)                        cGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCg                                                                                                                            >  6:93739/1‑149 (MQ=255)                           ccATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGGCGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAgc                                                                                                                         >  7:237658/1‑149 (MQ=255)                           ccATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGATGCA‑CGCCGAgc                                                                                                                         <  8:114639/149‑1 (MQ=255)                           ccATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAgc                                                                                                                         >  1:53543/1‑149 (MQ=255)                              tCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAgcgcc                                                                                                                      <  1:330677/149‑1 (MQ=255)                                 ccATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCATCGCCGAGCCCACGCTGCGCTCGGCG‑tt                                                                                                                                    >  7:162743/1‑131 (MQ=255)                                 ccATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCATCGCCGAGCCCACGCTGCGCTCGGCG‑tt                                                                                                                                    <  8:162743/132‑2 (MQ=255)                                 ccATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAgcgccgcg                                                                                                                   >  6:185950/1‑149 (MQ=255)                                          cgcCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGCTCAGTTGAGTGCTGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGaa                                                                                                           <  1:279463/149‑1 (MQ=255)                                             cATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCa                                                                                                       >  5:44828/1‑149 (MQ=255)                                                  aCTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGcc                                                                                                  >  6:87262/1‑149 (MQ=255)                                                    ttCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTg                                                                                                >  6:330886/1‑149 (MQ=255)                                                      cGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTg                                                                                                <  8:90792/147‑1 (MQ=255)                                                              accaGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGc                                                                                               <  5:251581/140‑1 (MQ=255)                                                              accaGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGc                                                                                               >  6:251581/1‑140 (MQ=255)                                                              accaGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCca                                                                                      >  7:76511/1‑149 (MQ=255)                                                              accaGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCca                                                                                      >  7:158155/1‑149 (MQ=255)                                                              accaGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCca                                                                                      >  3:26459/1‑149 (MQ=255)                                                              accaGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCca                                                                                      >  8:44427/1‑149 (MQ=255)                                                              accaGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGACCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCca                                                                                      <  1:109229/149‑1 (MQ=255)                                                                  gTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCcacgca                                                                                  >  5:271471/1‑149 (MQ=255)                                                                          cccATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCCCGCACGGCCtt                                                                           >  2:248718/1‑148 (MQ=255)                                                                          cccATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCtt                                                                           <  6:44828/148‑1 (MQ=255)                                                                             aTGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGGGGTGCGCCACGCACGGCCTTCATg                                                                       <  1:304428/149‑1 (MQ=255)                                                                             aTGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGGCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGGGGTGCGCCACGCACGGCCTTCATg                                                                       <  1:308424/149‑1 (MQ=255)                                                                              tGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATgc                                                                      <  7:44427/149‑1 (MQ=255)                                                                                cggcgCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATgcgc                                                                    <  5:185950/149‑1 (MQ=255)                                                                                   cgcTGGCGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCttt                                                                 <  8:237658/149‑1 (MQ=255)                                                                                        acgcg‑acgcTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGc                                                              <  8:76511/147‑1 (MQ=255)                                                                                        acgcg‑acgcTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGc                                                              <  8:158155/147‑1 (MQ=255)                                                                                                gcTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTg                                                     >  3:125727/1‑149 (MQ=255)                                                                                                gcTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTg                                                     >  4:24213/1‑149 (MQ=255)                                                                                                  ttGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACAGGGTgcg                                                   >  1:2815/1‑149 (MQ=255)                                                                                                    gAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTgcgc                                                  <  4:26459/148‑1 (MQ=255)                                                                                                         cGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGTCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGACTTCATGCGCTTTGGCGACCGGGTGCGCATGGa                                             <  5:93739/148‑1 (MQ=255)                                                                                                          gggCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAgg                                           <  5:90691/149‑1 (MQ=255)                                                                                                                gCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGCGCTCGGCGATGCA‑CGCCGAgccca                                                                                                                      <  3:273120/68‑4 (MQ=25)                                                                                                                gCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGCGCTCGGCGATGCA‑CGCCGAgccca                                                                                                                      >  4:273120/1‑65 (MQ=25)                                                                                                                    cTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGc                                  <  5:256567/148‑1 (MQ=255)                                                                                                                         cTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGTGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCAGGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAgg                            <  2:2815/149‑1 (MQ=255)                                                                                                                         cTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAgg                            <  6:330509/149‑1 (MQ=255)                                                                                                                         cTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAgg                            >  5:330509/1‑149 (MQ=255)                                                                                                                              gCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATCGAGGCTCGCGGGCCCGAGGGCGAg                       >  5:193860/1‑149 (MQ=255)                                                                                                                                 ccGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTg                    <  4:77200/149‑1 (MQ=255)                                                                                                                                  cGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGc                   >  5:313931/1‑149 (MQ=255)                                                                                                                                        cAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTg             >  2:331554/1‑149 (MQ=255)                                                                                                                                         aaCGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTgg            >  8:179464/1‑149 (MQ=255)                                                                                                                                           cGCCGAGCCCACGCTGCGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTgggg          >  4:310483/1‑149 (MQ=255)                                                                                                                                            gCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGc         <  5:330886/149‑1 (MQ=255)                                                                                                                                               gAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCAt      <  2:53543/149‑1 (MQ=255)                                                                                                                                               gAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCAt      <  7:179464/149‑1 (MQ=255)                                                                                                                                                aGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGGGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCAtt     <  3:24213/149‑1 (MQ=255)                                                                                                                                                   gcAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCATTGAc  <  4:125727/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCATTGAC  >  NC_002947/3638668‑3638961

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