Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,091,344	A→G	9.2%	intergenic (+30/‑48)	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,344	0	A	G	9.2%	178.0 / 7.3	65	intergenic (+30/‑48)	gudP/garD	glucarate transporter/galactarate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base G (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 = 9.85e-01

CGGTATTCAGCTATCTGGTTATCGTCGGCCCCATCAAGCGTATCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCG  >  NC_002947/4091196‑4091491                                                                                                                                                     |                                                                                                                                                    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                                                                                                                   <  5:266645/134‑1 (MQ=255)                                                  gCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACaa                                                                                                                   >  6:266645/1‑134 (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:59445/1‑149 (MQ=255)                                                                  aaGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACAt                                                                                    >  6:295498/1‑149 (MQ=255)                                                                    tccACACGCTGATCCTGCCACCGACGGCGTACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTc                                                                                  <  4:231082/148‑1 (MQ=255)                                                                     ccACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATc                                                                                          >  4:227479/1‑140 (MQ=255)                                                                     ccACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATc                                                                                          <  3:227479/140‑1 (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                                                                                                                   <  3:329432/107‑1 (MQ=255)                                                                             tGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACaa                                                                                                                   >  4:329432/1‑107 (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                                                                                          <  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)                                                                                                                                           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)                                                                                                                                                     |                                                                                                                                                    CGGTATTCAGCTATCTGGTTATCGTCGGCCCCATCAAGCGTATCGAATTGCGTGATGCTGCCCCGAAGCCACACGCTGAGCCTGCCACCGACGGCGAACTGGCCAGCTCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCG  >  NC_002947/4091196‑4091491

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