Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,091,345	T→C	9.5%	intergenic (+31/‑47)	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,345	0	T	C	9.5%	171.2 / 7.7	65	intergenic (+31/‑47)	gudP/garD	glucarate transporter/galactarate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (27/30);  new base C (3/3);  total (31/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 = 4.95e-01

CAGCTATCTGGTTATCGTCGGCCCCATCAAGCGTATCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCG  >  NC_002947/4091203‑4091491                                                                                                                                               |                                                                                                                                                   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                                                                                                                                          >  1:314280/1‑100 (MQ=255)                                                      ccccGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGt                                                                                                                                          <  2:314280/100‑1 (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                                                <  5:217588/144‑1 (MQ=255)                                                                                                    cTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCacgacgac                                                >  6:217588/1‑144 (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                                                                                          >  5:338590/1‑93 (MQ=255)                                                                                                             cTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATc                                                                                          <  6:338590/93‑1 (MQ=255)                                                                                                              tGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCacgac                                                   >  8:284267/1‑131 (MQ=255)                                                                                                              tGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCacgac                                                   <  7:284267/131‑1 (MQ=255)                                                                                                                aaaaaaCACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAATCATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGgtcgt                               >  2:283631/1‑149 (MQ=255)                                                                                                                aaaaaaCACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAAcgtcg                                      >  2:93917/1‑142 (MQ=255)                                                                                                                aaaaaaCACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAAcgtcg                                      <  1:93917/142‑1 (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)                                                                                                                                    gTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACACTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTgg           <  6:362051/149‑1 (MQ=255)                                                                                                                                      gacccGGCGCGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGc         <  3:155526/147‑1 (MQ=255)                                                                                                                                         ccGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGCCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCTTCAACGACGGCGGCTTGGGCGaa      <  6:358849/149‑1 (MQ=255)                                                                                                                                          cGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACCGCGGCTTGGGCGAAg     <  7:61693/149‑1 (MQ=255)                                                                                                                                             cATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGgcg  <  3:334049/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                   CAGCTATCTGGTTATCGTCGGCCCCATCAAGCGTATCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCG  >  NC_002947/4091203‑4091491

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