Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,638,821	G→C	9.4%	G273R (GGC→CGC)	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,821	0	G	C	9.4%	193.3 / 9.0	64	G273R (GGC→CGC)	PP_3205	fumarylacetoacetate hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base C (3/3);  total (30/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 = 6.43e-01

TGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCATTGACCAGCGCG  >  NC_002947/3638674‑3638968                                                                                                                                                     |                                                                                                                                                     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)                  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                                                                                                      >  6:251581/1‑140 (MQ=255)                                                        accaGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGc                                                                                                      <  5:251581/140‑1 (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‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCca                                                                                             >  7:76511/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)                                                                                                                                                gCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCATTGACCAg      >  1:284586/1‑149 (MQ=255)                                                                                                                                                 cGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCgg                                           <  8:130704/111‑1 (MQ=255)                                                                                                                                                 cGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCgg                                           >  7:130704/1‑111 (MQ=255)                                                                                                                                                   tGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCATTGACCAgcgc   <  6:313931/149‑1 (MQ=255)                                                                                                                                                    gggCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCATTGACCAGcgcg  >  3:104537/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCG‑ACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCG‑TGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCATTGACCAGCGCG  >  NC_002947/3638674‑3638968

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