Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,091,328	C→T	7.5%	intergenic (+14/‑64)	gudP → / → garD	glucarate transporter/galactarate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,091,328	0	C	T	7.5%	193.7 / 3.8	66	intergenic (+14/‑64)	gudP/garD	glucarate transporter/galactarate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base T (3/2);  total (34/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 = 3.34e-01

AACGCGCTGGTCGCGGTATTCAGCTATCTGGTTATCGTCGGCCCCATCAAGCGTATCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACG  >  NC_002947/4091183‑4091473                                                                                                                                                  |                                                                                                                                                  aaCGCGCTGGTCGCGGTATTCAGCTATCTGGTTATCGTCGGCCCCATCAAGCGTATCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGcc                                                                                                                                                >  7:186227/1‑149 (MQ=255) aaCGCGCTGGTCGCGGTATTCAGCTATCTGGTTATCGTCGGCCCCATCAAGCGTATCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGcc                                                                                                                                                >  5:359439/1‑149 (MQ=255)    gcgcTGGTCGCGGTATTCAGCTATCTGGTTATCGTCGGCCCCATCAAGCGTATCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCggg                                                                                                                                             <  6:342284/149‑1 (MQ=255)              cgGTATTCAGCTATCTGGTTATCGTCGGCCCCATCAAGCGTATCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACATGCCCGCGCCGGGTCGCCCGGCa                                                                                                                                   <  7:53367/149‑1 (MQ=255)                     cAGCTATCTGGTTATCGTCGGCCCCATCAAGCGTATCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCg                                                                                                                             <  8:186227/148‑1 (MQ=255)                      aGCTATCTGGTTATCGTCGGCCCCATCAAGCGTATCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCATGCCGGGCGACCCGGCGCGGGCgt                                                                                                                            <  5:207140/148‑2 (MQ=255)                                  aTCGTCGGCCCCATCAAGCGTATCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGc                                                                                                                <  6:359439/148‑1 (MQ=255)                                   tcgtcgGCCCCATCAAGCGTATCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCt                                                                                                              <  5:147293/149‑1 (MQ=255)                                             cATCAAGCGTATCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATggg                                                                                                                               <  4:326832/122‑1 (MQ=255)                                             cATCAAGCGTATCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATggg                                                                                                                               >  3:326832/1‑122 (MQ=255)                                                     gTAACGAATTGCGTGATCCTGCCCCGAAGCCACACGCAGACCCACCCACCGACCGCGAACTGGCCAGCTCCCGGCACTGAAAAAACACGCCCGCGCCGGCACCCACGGCAAGGGCGGATCTACAAAAGCCTGAGACCCTGACAAGTAgg                                                                                            >  4:244611/1‑149 (MQ=255)                                                      tATCGAATTGCGTGATGCTGCCCCTAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTTAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCt                                                                                                              <  6:243374/130‑1 (MQ=255)                                                      tATCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCTAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAggg                                                                                           <  3:239474/149‑1 (MQ=255)                                                       aTCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGc                                                                                          >  6:345686/1‑149 (MQ=255)                                                          gAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAGGCATGAGAACCTGACAAGTAGGGccc                                                                                        >  3:231082/1‑148 (MQ=255)                                                               gCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACaa                                                                                                 >  6:266645/1‑134 (MQ=255)                                                               gCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACaa                                                                                                 <  5:266645/134‑1 (MQ=255)                                                                cGTGATGCTGCCCCGAAGCCACAAGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGACACTGAAAAAACACGCCCGCGCCGGGTCGCCCCGCATGGGCGGATGTTCAAAAGCCTGAGAACCTGACAAGTAGGGACCTACCATg                                                                                 >  8:61693/1‑149 (MQ=255)                                                                          ccccGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGt                                                                                                                        <  2:314280/100‑1 (MQ=255)                                                                          ccccGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGt                                                                                                                        >  1:314280/1‑100 (MQ=255)                                                                             cGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCa                                                                                                                                   <  1:77238/86‑1 (MQ=255)                                                                             cGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCa                                                                                                                                   >  2:77238/1‑86 (MQ=255)                                                                               aaGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACAt                                                                  >  6:295498/1‑149 (MQ=255)                                                                               aaGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACAt                                                                  >  6:59445/1‑149 (MQ=255)                                                                                 tccACACGCTGATCCTGCCACCGACGGCGTACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTc                                                                <  4:231082/148‑1 (MQ=255)                                                                                  ccACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATc                                                                        <  3:227479/140‑1 (MQ=255)                                                                                  ccACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATc                                                                        >  4:227479/1‑140 (MQ=255)                                                                                  ccACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTcc                                                               <  5:345686/149‑1 (MQ=255)                                                                                    acacGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGcc                                                                                         <  5:54973/121‑1 (MQ=255)                                                                                    acacGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCg                                                              >  4:334049/1‑148 (MQ=255)                                                                                    acacGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTCACAAGTAGGGcc                                                                                         >  6:54973/1‑121 (MQ=255)                                                                                      acGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACt                                                           <  5:295498/149‑1 (MQ=255)                                                                                          tGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACaa                                                                                                 >  4:329432/1‑107 (MQ=255)                                                                                          tGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACaa                                                                                                 <  3:329432/107‑1 (MQ=255)                                                                                                tGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTa                                                 >  5:362051/1‑149 (MQ=255)                                                                                                 gCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTAc                                                <  2:108858/149‑1 (MQ=255)                                                                                                      cGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCg                                           >  7:16463/1‑149 (MQ=255)                                                                                                         cGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCt                                        <  8:292763/149‑1 (MQ=255)                                                                                                             gAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCac                                    >  5:139983/1‑149 (MQ=255)                                                                                                                     cAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCAcgacgacga                            >  5:358849/1‑149 (MQ=255)                                                                                                                      aGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACgacgacgac                           <  6:139983/149‑1 (MQ=255)                                                                                                                        cTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCacgacgac                              >  6:217588/1‑144 (MQ=255)                                                                                                                        cTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCacgacgac                              <  5:217588/144‑1 (MQ=255)                                                                                                                         tCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAAc                        >  8:282928/1‑149 (MQ=255)                                                                                                                           ccGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCa                                     >  6:195996/1‑134 (MQ=255)                                                                                                                           ccGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCa                                     <  5:195996/134‑1 (MQ=255)                                                                                                                           ccGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAAcgt                      >  4:236384/1‑149 (MQ=255)                                                                                                                            cGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAAcgtc                     <  7:282928/149‑1 (MQ=255)                                                                                                                                aCTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTgg                 <  5:59445/149‑1 (MQ=255)                                                                                                                                aCTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTg                  <  8:16463/148‑1 (MQ=255)                                                                                                                                 cTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATc                                                                        <  6:338590/93‑1 (MQ=255)                                                                                                                                 cTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATc                                                                        >  5:338590/1‑93 (MQ=255)                                                                                                                                  tGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCacgac                                 >  8:284267/1‑131 (MQ=255)                                                                                                                                  tGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCacgac                                 <  7:284267/131‑1 (MQ=255)                                                                                                                                    aaaaaaCACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAATCATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGgtcgt             >  2:283631/1‑149 (MQ=255)                                                                                                                                    aaaaaaCACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAAcgtcg                    <  1:93917/142‑1 (MQ=255)                                                                                                                                    aaaaaaCACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAAcgtcg                    >  2:93917/1‑142 (MQ=255)                                                                                                                                        aaCACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGgtcgtcgt          >  6:361881/1‑148 (MQ=255)                                                                                                                                         atccGCCCATGCCGGGCGACCCGGCGCGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGgtcgtcgtc         >  3:173252/4‑148 (MQ=255)                                                                                                                                         atccGCCCATGCCGGGCGACCCGGCGCGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGgtcgtcgtc         >  4:155526/4‑148 (MQ=255)                                                                                                                                          tccGCCCATGCCGGGCGACCCGGCGCGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATg                                                                                 >  7:105124/3‑75 (MQ=18)                                                                                                                                          tccGCCCATGCCGGGCGACCCGGCGCGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATg                                                                                 <  8:105124/73‑1 (MQ=18)                                                                                                                                           acGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCaa       >  6:11923/1‑148 (MQ=255)                                                                                                                                           acGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCaa       <  5:11923/148‑1 (MQ=255)                                                                                                                                              cccGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAacgac   >  7:9474/1‑149 (MQ=255)                                                                                                                                               ccGCGCCGGGTCGCCCGGGAGGTGGGGTTGTACAAAAGCCTGAGAACCTGCCAAGTAGGTCCCAAACATGGAGTTGATCGAACATTCGGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAacgacg  <  3:244611/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  AACGCGCTGGTCGCGGTATTCAGCTATCTGGTTATCGTCGGCCCCATCAAGCGTATCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACG  >  NC_002947/4091183‑4091473

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