Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,186,052	C→A	5.3%	E136D (GAG→GAT) ‡	aceK ←	isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,186,052	0	C	A	5.3%	249.7 / 2.7	76	E136D (GAG→GAT) ‡	aceK	isocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (31/41);  new base A (2/2);  total (33/43)
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.49e-01

GCATTCGCGCAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATCAGGTG  >  NC_002947/5185908‑5186192                                                                                                                                                  |                                                                                                                                              gcaTTCGCGCAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTcgc                                                                                                                                           >  7:141611/1‑149 (MQ=255)       gcgcAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCTCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCcgc                                                                                                                                    <  8:287903/149‑1 (MQ=255)       gcgcAGTTGTTCTTCCAGCCGCGCCAGGGCGCCATCCAGGTCTCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCcgc                                                                                                                                    <  8:282754/149‑1 (MQ=255)          cAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCa                                                                                                                                 <  6:191669/149‑1 (MQ=255)           aGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCAt                                                                                                                                >  8:298830/1‑149 (MQ=255)               gttCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCgcgc                                                                                                                                  >  5:162982/1‑143 (MQ=255)               gttCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCgcgc                                                                                                                                  <  6:162982/143‑1 (MQ=255)                ttcttcCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGaa                                                                                                                            <  5:313440/148‑1 (MQ=255)                   ttcCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGcc                                                                                                                        >  2:161214/1‑149 (MQ=255)                   ttcCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGcc                                                                                                                        >  4:110112/1‑149 (MQ=255)                     ccAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCg                                                                                                                       <  1:161214/148‑1 (MQ=255)                      cAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCggg                                                                                                                     >  6:386021/1‑149 (MQ=255)                      cAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCgg                                                                                                                      <  5:19081/148‑1 (MQ=255)                        gCCGCGCCAGGTCGCCTTCCCGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCgcgc                                                                                                                                  <  1:23783/134‑1 (MQ=255)                        gCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCgcgc                                                                                                                                  >  2:23783/1‑134 (MQ=255)                        gCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTa                                                                                                                   <  3:110112/149‑1 (MQ=255)                         ccGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAg                                                                                                                  <  5:84576/149‑1 (MQ=255)                               cAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGa                                                                                                            <  2:282819/149‑1 (MQ=255)                                ctgcagCCCTCCCGCTCGCCCTAGCGCATGTCCAACGGGGACTCGGCGCACCCCGCGCGCAGCACGTTCCTGCCACCGCCGTCCAGCCGATCGGTCCGGGTGTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAACGCCGGGTAGTGCGGAt                                                                                                           <  8:348848/144‑1 (MQ=255)                                atgTCGCCGTCCTGGTCGCCATAGGGCATGTCGAACGGGTAATCTGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGAt                                                                                                           <  4:55078/147‑1 (MQ=255)                                  gTCGCCTTCCAGGTAGCCATAGGGCATGTCGTACGGGTAATCGGCGAACAGCGCGCGCAGCAAGTTCCTGAGACTGCCGTCTAGCCGATAGGTACGGGTTTGCGGGG‑CGCGCTCGCG‑GCCGCGGATGGAAGGCCGGGTAGTGTGGATGa                                                                                                         <  4:145931/149‑1 (MQ=255)                                   tCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAgg                                                                                                                          >  5:231420/1‑131 (MQ=255)                                   tCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAgg                                                                                                                          <  6:231420/131‑1 (MQ=255)                                   tCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTATTGTGGAt                                                                                                           >  7:348848/1‑146 (MQ=255)                                     gCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCgg                                                                                                                      <  3:372437/133‑1 (MQ=255)                                     gCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCgg                                                                                                                      >  4:372437/1‑133 (MQ=255)                                     gCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAAc                                                                                                       <  3:145566/148‑1 (MQ=255)                                            aGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑C‑CGCGAGCGCGCCGCGCaaaccc                                                                                                                            >  7:186985/1‑115 (MQ=38)                                            aGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑C‑CGCGAGCGCGCCGCGCaaaccc                                                                                                                            <  8:186985/120‑6 (MQ=38)                                                          cATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAgg                                                                                  >  7:190729/1‑148 (MQ=255)                                                           aagTCGAACGTGTAATCGTCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCTTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGt                                                                                 <  8:141611/146‑1 (MQ=255)                                                                gAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATgg                                                                            <  7:298830/148‑1 (MQ=255)                                                                gAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGc                                                                           >  7:347792/1‑149 (MQ=255)                                                                  aCGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTg                                                                         <  8:190729/149‑1 (MQ=255)                                                                        aaTCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCACATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATc                                                                                     <  7:281352/131‑1 (MQ=255)                                                                        aaTCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCACATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATc                                                                                     >  8:281352/1‑131 (MQ=255)                                                                          tCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTa                                                                                                                   >  4:31133/1‑99 (MQ=255)                                                                          tCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCt                                                                 >  4:95285/1‑149 (MQ=255)                                                                           cGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAgg                                                                                                                          <  4:83963/91‑1 (MQ=255)                                                                           cGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAgg                                                                                                                          >  3:83963/1‑91 (MQ=255)                                                                             gCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGc                                                               >  1:8328/1‑148 (MQ=255)                                                                             gCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCAGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCTTCGCTGATCAGGTCATGGCTGAAGAGGCTGc                                                               >  8:173549/1‑148 (MQ=255)                                                                               gAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCaaa                                                            <  5:386021/149‑1 (MQ=255)                                                                                aaCACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCaaaa                                                           >  1:15265/1‑149 (MQ=255)                                                                                 acacCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCaaaa                                                           >  7:108891/1‑148 (MQ=255)                                                                                  cacCGCGCGCAGCAAGTTCCTGAGACAGCCGTCAAGCCGATAGCTACGGGTTTGCGCGG‑CGCGCTCGAG‑GCCGCGCATCCAAGGCCCGGTAGAGTGGAAGAACAAGCAGCCGTCGCTGAACAGGTCAAGGCTGAAGACGCTGCAAAACa                                                         >  6:69482/1‑149 (MQ=255)                                                                                   acCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATgg                                                                            <  6:258069/129‑1 (MQ=255)                                                                                   acCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAg                                                        >  4:211704/1‑149 (MQ=255)                                                                                     cgcgcgCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGAAAAACAGc                                                       >  3:91391/1‑148 (MQ=255)                                                                                              caAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTAcc                                             <  6:388913/149‑1 (MQ=255)                                                                                              caAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTAc                                              >  2:56595/1‑148 (MQ=255)                                                                                              caAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTAc                                              >  4:31716/1‑148 (MQ=255)                                                                                               aaGTTCCTGAGACTGCCGTCAAGCCGATAGGTTCGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCa                                            <  5:69482/149‑1 (MQ=255)                                                                                               aaGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCa                                            <  2:129919/149‑1 (MQ=255)                                                                                                aGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAg                                           <  5:127973/149‑1 (MQ=255)                                                                                                  ttCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGt                                                   <  4:22105/139‑1 (MQ=255)                                                                                                  ttCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGt                                                   >  3:22105/1‑139 (MQ=255)                                                                                                      tGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTctct                                     >  1:249755/1‑149 (MQ=255)                                                                                                        cgaCTGCCGTCATGGCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GGCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGc                                   <  4:56008/148‑1 (MQ=255)                                                                                                        agaCTTCCGTCAAGCCTAGAGTTACGTGTTTGCGCGG‑CGAGCTCGCG‑TCCGCTCATGGAAGGCCGGGTAGTGTGTATTAACATGCAGCCGTCGCTGATCAGGTCATGTCTGAAGAGGCTGCAAAACAGCTAGTTGTACCAGGTctct                                     <  4:91391/147‑1 (MQ=255)                                                                                                        agaCTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTctct                                     <  3:95285/147‑1 (MQ=255)                                                                                                         gaCTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGcc                                  <  3:211704/149‑1 (MQ=255)                                                                                                          aCTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGcc                                  >  5:257426/1‑148 (MQ=255)                                                                                                              ccGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACtc                             <  2:8328/149‑1 (MQ=255)                                                                                                              ccGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACtc                             <  8:347792/149‑1 (MQ=255)                                                                                                               cGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACtcg                            <  2:15265/149‑1 (MQ=255)                                                                                                                      ccGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCa                     <  2:249755/149‑1 (MQ=255)                                                                                                                        gATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAGACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAgg                   >  5:342879/1‑149 (MQ=255)                                                                                                                               aCGCGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGACGGGTAGTGTGGAAGAACATGCAGCCGTCGCTGATCATGTCATGGATGAAGAGGCTACAAAACAGCGAATTGTACGAGGTATCTGACAACTCGACGTCCAGGCGGGGGt            >  2:268989/1‑149 (MQ=255)                                                                                                                                cGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGt            <  3:31716/148‑1 (MQ=255)                                                                                                                                    tttGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAATACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATc       <  6:342879/149‑1 (MQ=255)                                                                                                                                     atGCGCGGCCGCGAGCGC‑‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGtt                                                  <  2:135566/104‑1 (MQ=255)                                                                                                                                     atGCGCGGCCGCGAGCGC‑‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGtt                                                  >  1:135566/2‑105 (MQ=255)                                                                                                                                      tGCGCGGCCGCGAGCGC‑‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCgggg              >  4:193199/1‑140 (MQ=255)                                                                                                                                      tGCGCGGCCGCGAGCGC‑‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCgggg              <  3:193199/140‑1 (MQ=255)                                                                                                                                         gcgg‑cgcgCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATCAGGTg  <  6:257426/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                              GCATTCGCGCAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGG‑CGCGCTCGCG‑GCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATCAGGTG  >  NC_002947/5185908‑5186192

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