Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,662,394	A→C	7.8%	E166D (GAA→GAC)	mmgF →	2‑methylisocitrate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,662,394	0	A	C	7.8%	183.5 / 4.9	61	E166D (GAA→GAC)	mmgF	2‑methylisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base A (27/29);  new base C (2/3);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

CCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTACAAGACTTTCGCTGATCGGGTGAAGGCACCGATCCTGGCCAAC  >  NC_002947/2662247‑2662529                                                                                                                                                     |                                                                                                                                        ccAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAg                                                                                                                                        <  5:90583/149‑1 (MQ=255) ccAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAg                                                                                                                                        <  1:142734/149‑1 (MQ=255) ccAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAg                                                                                                                                        <  5:231257/149‑1 (MQ=255)   aGGTGGGCGCCACGCGGTGCGGCCACCGTCCAACCAAGGAAATGGTCAGCCAGCCCGAATTGGGCGACCGTTTCAAAGCTGCCGGCGACGCCCCC‑AGCGATGACAGCTTCGTGATCTTGGCGCGTACCGATGCGCTGGCGGTCGAAGGc                                                                                                                                      <  6:38946/149‑1 (MQ=255)     gTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCt                                                                                                                                    >  1:197528/1‑149 (MQ=255)      ttGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCTTCGACGGccag                                                                                                                                   <  6:257836/149‑3 (MQ=255)         gcgcCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCTTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAAc                                                                                                                                <  4:244441/149‑1 (MQ=255)         gcgataaGCGTTGCGGCCACAGTACCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAAc                                                                                                                                <  8:138948/144‑1 (MQ=255)           gcCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGCAATCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCTTCGACGGccagcgca                                                                                                                               <  1:267139/149‑7 (MQ=255)                     gCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGa                                                                                                                     <  1:281184/148‑1 (MQ=255)                         ccACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCTTCGACGGccagc                                                                                                                                  <  7:242155/131‑4 (MQ=255)                         ccACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCTTCGACGGccagc                                                                                                                                  >  8:242155/1‑128 (MQ=255)                          cACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGc                                                                                                               >  2:231363/1‑149 (MQ=255)                           aCCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGcc                                                                                                              <  6:248928/149‑1 (MQ=255)                                 ccAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCg                                                                                                        >  8:207205/1‑149 (MQ=255)                                   aacaaGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCTTGAACGCCGCCCTGGACCGTGCCCAGGCgtg                                                                                                      >  3:136345/1‑149 (MQ=255)                                   aacaaGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCt                                                                                                                                    >  4:179289/1‑119 (MQ=255)                                   aacaaGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCt                                                                                                                                    <  3:179289/119‑1 (MQ=255)                                   aacaaGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCgtg                                                                                                      >  2:47266/1‑149 (MQ=255)                                      aaGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCgtgtgt                                                                                                   >  7:161244/1‑149 (MQ=255)                                        ggAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGTCCAGGCGTGTGTc                                                                                                  >  2:30757/1‑148 (MQ=255)                                        ggAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTc                                                                                                  >  6:211099/1‑148 (MQ=255)                                                tCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTgg                                                                                          <  1:47266/148‑1 (MQ=255)                                                          ggAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACgccgcc                                                                                                                          <  6:225318/106‑1 (MQ=255)                                                          ggAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACgccgcc                                                                                                                          >  5:225318/1‑106 (MQ=255)                                                           gAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGa                                                                               <  1:224841/148‑1 (MQ=255)                                                           gAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGa                                                                               >  2:224841/1‑148 (MQ=255)                                                             aaTGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCt                                                                            >  6:220334/1‑149 (MQ=255)                                                               tGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTc                                                                          <  1:155753/149‑1 (MQ=255)                                                                 gTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGgcg                                                                                        <  6:336784/133‑1 (MQ=255)                                                                 gTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGgcg                                                                                        >  5:336784/1‑133 (MQ=255)                                                                 gTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTcc                                                                         <  5:220334/148‑1 (MQ=255)                                                                  tCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCg                                                                       >  4:27466/1‑149 (MQ=255)                                                                    gACCGTATCAAAGCTGCCGTCGCCGCCCGC‑AGCGATGACAGCTTGGTGATGATGGCGCGTACCGATGCGCGGGCCGTCGACGGCCGGAGCGCCGCCCTGGACGGTGCCGAGGCGTGTGTCGAAGGGGGGGGCGCCGTGATCTTCCCgg                                                                      >  8:190322/1‑148 (MQ=255)                                                                       cGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGc                                                                  >  7:67377/1‑149 (MQ=255)                                                                         tgtcatcgCTGCGGGCGTCGACGGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGcc                                                                 >  3:75166/8‑148 (MQ=255)                                                                          aTCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCAt                                                               >  8:6601/1‑149 (MQ=255)                                                                              aaGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCAcc                                                           <  8:161244/149‑1 (MQ=255)                                                                                 cTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGaa                                                        >  5:83939/1‑149 (MQ=255)                                                                                      gTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCa                                                   <  4:11157/149‑1 (MQ=255)                                                                                         gccGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAg                                                                                              >  2:109074/3‑103 (MQ=255)                                                                                         gccGCCCGC‑AGCGATGACAGCTTCGAGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAg                                                                                              <  1:109074/101‑1 (MQ=255)                                                                                           cGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGACCTGCAGATGt                                              <  4:194230/149‑1 (MQ=255)                                                                                              ccGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTAc                                            >  3:247583/1‑148 (MQ=255)                                                                                              ccGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGACGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTAc                                            >  2:1212/1‑148 (MQ=255)                                                                                               cGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTACaa                                          <  4:135406/149‑1 (MQ=255)                                                                                                        tGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTACAAGACTTTc                                   <  4:75166/148‑1 (MQ=255)                                                                                                        tGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTACAAGACTTTCg                                  <  1:231363/149‑1 (MQ=255)                                                                                                           cAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCAACGAACTGCAGATGTACAAGACTTTCGCt                                >  6:77655/1‑148 (MQ=255)                                                                                                              cTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCTTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTACAAGACTTTCGCTGATc                            <  4:136345/149‑1 (MQ=255)                                                                                                                 cGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTACAAGACTTTCGCTGATCgg                          <  7:207205/148‑1 (MQ=255)                                                                                                                 cGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTACAAGACTTTCGCTGATCgg                          <  2:197528/148‑1 (MQ=255)                                                                                                                     aTCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTACAAGACTTTCGCTGATCGGGTg                       >  6:110218/1‑147 (MQ=255)                                                                                                                         tGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTACAAGACTTTCGCTGATCGGGTGAAGGCa                 >  3:94969/1‑149 (MQ=255)                                                                                                                           gcgcGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGACGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTACAAGACTTTCGCTGATCGGGTGAAGGCAc                <  1:1212/148‑1 (MQ=255)                                                                                                                             gcgTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTACAAGACTTTCGCTGATCGGGt                        >  6:292285/1‑138 (MQ=255)                                                                                                                             gcgTACCGATGCGCTGGCCGTCGAAGGCCTCAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTACAAGACTTTCGCTGATCGGGt                        <  5:292285/138‑1 (MQ=255)                                                                                                                                   cGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTACAAGACTTTCGCTGATCGGGTGAAGGCACCGATCCTgg       >  2:180807/1‑149 (MQ=255)                                                                                                                                        cgcTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCTTTTTTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTACAAGACTTTCGCTGATCGGGTGAAGGCACCGATCCTGGCCAAc  <  7:190322/149‑1 (MQ=255)                                                                                                                                         gcTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTACAAGACTTTCGCTGATc                            <  7:250165/122‑1 (MQ=255)                                                                                                                                         gcTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTACAAGACTTTCGCTGATc                            >  8:250165/1‑122 (MQ=255)                                                                                                                                                     |                                                                                                                                        CCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGC‑AGCGATGACAGCTTCGTGATCATGGCGCGTACCGATGCGCTGGCCGTCGAAGGCCTGAACGCCGCCCTGGACCGTGCCCAGGCGTGTGTCGAAGCTGGCGCCGACATGATCTTCCCGGAAGCCATCACCGAACTGCAGATGTACAAGACTTTCGCTGATCGGGTGAAGGCACCGATCCTGGCCAAC  >  NC_002947/2662247‑2662529

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