Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	648,289	A→T	10.1%	intergenic (+58/+33)	acoR → / ← accC	acetoin catabolism regulatory protein/acetyl‑CoA carboxylase biotin carboxylase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	648,289	0	A	T	10.1%	157.0 / 7.0	59	intergenic (+58/+33)	acoR/accC	acetoin catabolism regulatory protein/acetyl‑CoA carboxylase biotin carboxylase subunit
Reads supporting (aligned to +/- strand):  ref base A (19/33);  new base T (2/4);  total (22/37)
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.65e-01

GATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGTC  >  NC_002947/648142‑648435                                                                                                                                                     |                                                                                                                                                   gATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGCCGCCGCCCCCCA‑ATATGGGAGGCAATAAAAAAACCCCGcacac                                                                                                                                                   >  8:236971/1‑148 (MQ=255)  aTGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAg                                                                                                                                                   <  3:160963/148‑1 (MQ=255)      ggCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGCCGCGGCCCCACA‑ATATGGGAGGCAATAAAAAAACACCGCACACGGccc                                                                                                                                              >  6:3580/1‑148 (MQ=255)        cACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCggg                                                                                                                                            <  2:261798/149‑1 (MQ=255)          cGGGCTTTGGGCATTTCGCGGGCCACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACCTGCCGCGGATGAAGGGGTCGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAGAACCCCGGCCCTGTGCCGGGGt                                                                                                                                          <  5:346920/149‑1 (MQ=255)           gggCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGGCGCGCCCCGGCCCCAAAGTGAGGGCGGGGAAAAAAAAACCCCCCCCACGGGCCggggg                                                                                                                                          >  3:210402/1‑148 (MQ=255)           gggCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGCCGCGGCCCCCGA‑TTATGGGAGGCAATAAAAAACCCCCGCACAGGGCCgggggt                                                                                                                                         >  4:102801/1‑147 (MQ=255)            ggCTTTGGGTATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATGGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTc                                                                                                                                        <  6:91191/149‑1 (MQ=255)               tttGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGGCCCGCCCCCGCCCCCCC‑ATATGGCAGGCAAAAAAAAAACCCCCCCCACGGGCgggggtggtg                                                                                                                                     >  2:175622/1‑143 (MQ=255)                 tGGGCATTTGGCGGGCGACGCTTTATCGGCGATTGCATGATGACCGAATGGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGTGCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAAAAAAAAACCCCGGCCCTGTGCCGGGGtttttta                                                                                                                                    <  2:8464/148‑6 (MQ=255)                 tGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTgg                                                                                                                                     <  2:252880/147‑1 (MQ=255)                 tGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGt                                                                                                                                    <  3:132058/148‑1 (MQ=255)                 tGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGt                                                                                                                                    <  7:97273/148‑1 (MQ=255)                 tGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGt                                                                                                                                    <  3:306978/148‑1 (MQ=255)                   ggCATTTCGCGGGCGACGCTTTATCGGCTAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAGAACCCCGGCCCTGTGCCGGGGttttttatt                                                                                                                                  <  5:93264/148‑8 (MQ=255)                     cATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATc                                                                                                                               <  3:25352/149‑1 (MQ=255)                            cgGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAGCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTac                                                                                                                        <  1:175622/149‑1 (MQ=255)                                          tCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTgg                                                                                                          <  7:257859/149‑1 (MQ=255)                                                 gTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGGACAGGGCCGGGGGTGAGGTTCATAACACTACACGCTCACTCCGGGTTggg                                                                                                     >  5:334826/1‑146 (MQ=255)                                                                     cATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGa                                                                               <  3:116295/149‑1 (MQ=255)                                                                          cgcgGATGACGGGTGAGAGCCGTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGtt                                                                          <  2:318322/149‑1 (MQ=255)                                                                             ggATGAAGGTGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAc                                                                       <  7:236971/149‑1 (MQ=255)                                                                                  aaGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAgg                                                                                     <  5:315896/130‑1 (MQ=255)                                                                                   aGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAcgccgc                                                                  <  8:87675/148‑1 (MQ=255)                                                                                     gggTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCtt                                                               <  1:320037/149‑1 (MQ=255)                                                                                       gTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGttt                                                                                               <  1:116357/115‑1 (MQ=255)                                                                                          gagCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAg                                                          <  5:3580/149‑1 (MQ=255)                                                                                             cccTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAg                                                       >  5:150413/1‑149 (MQ=255)                                                                                              tcTGGGGCCGCGTCGCGGCCCCAGG‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTc                                                                                                              <  1:230160/92‑1 (MQ=255)                                                                                              tcTGGGGCCGCGTCGCGGCCCCAGG‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTc                                                                                                              >  2:230160/2‑93 (MQ=255)                                                                                              ccTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGc                                                      >  2:89755/1‑149 (MQ=255)                                                                                                  gggCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTTCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCtt                                                   >  6:229548/1‑148 (MQ=255)                                                                                                   ggCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGTCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtca                                                 <  1:89755/149‑1 (MQ=255)                                                                                                                ccccAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCtt                                                   <  6:289913/134‑1 (MQ=255)                                                                                                                ccccAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCtt                                                   >  5:289913/1‑134 (MQ=255)                                                                                                                 cccAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCa                                       >  6:325089/1‑145 (MQ=255)                                                                                                                 cccAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCa                                       <  5:325089/145‑1 (MQ=255)                                                                                                                   cAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcgc                                 <  6:260842/149‑1 (MQ=255)                                                                                                                    aGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcgcg                                <  5:304162/149‑1 (MQ=255)                                                                                                                     gA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCgg                               <  5:123649/149‑1 (MQ=255)                                                                                                                        ttCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTg                             >  7:112799/1‑149 (MQ=255)                                                                                                                               aggcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCgg                       >  3:96568/1‑148 (MQ=255)                                                                                                                                  caATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATg                   >  6:25417/1‑149 (MQ=255)                                                                                                                                   aaTAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCg                                                                                        <  7:10156/79‑1 (MQ=255)                                                                                                                                   aaTAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCg                                                                                        >  8:10156/1‑79 (MQ=255)                                                                                                                                   aaTAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGt                  >  7:268937/1‑149 (MQ=255)                                                                                                                                       cagaaCCCCGGCCCTGTGCCGGGGTTTTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGt              >  6:13708/4‑149 (MQ=255)                                                                                                                                        aaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATg                                                                            <  6:25368/86‑1 (MQ=255)                                                                                                                                        aaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATg                                                                            >  5:25368/1‑86 (MQ=255)                                                                                                                                        aaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAcgcc                                                                    >  6:224860/1‑94 (MQ=255)                                                                                                                                        aaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAcgcc                                                                    <  5:224860/94‑1 (MQ=255)                                                                                                                                        aaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTc             <  4:96568/149‑1 (MQ=255)                                                                                                                                         gaaCCCCGGCCCTGTGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCt            >  6:33525/2‑149 (MQ=255)                                                                                                                                            ccccGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGa         <  8:112799/149‑1 (MQ=255)                                                                                                                                             cccGGCCCTGTGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGAt        <  5:33525/149‑1 (MQ=255)                                                                                                                                             cccGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtc                                                  <  3:321300/107‑1 (MQ=255)                                                                                                                                             cccGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtc                                                  >  4:321300/1‑107 (MQ=255)                                                                                                                                             cccGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGAt        >  1:260126/1‑149 (MQ=255)                                                                                                                                                ggCAGAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGcc     <  6:334826/149‑1 (MQ=255)                                                                                                                                                 gCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCg    <  8:268937/149‑1 (MQ=255)                                                                                                                                                  cacaGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGt   <  3:102801/149‑1 (MQ=255)                                                                                                                                                  cacaGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGt   <  1:311738/149‑1 (MQ=255)                                                                                                                                                     aGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCGGTTGCAGAAGCCTTCATCACGCACCAGGTCGGGGGGGAGGGGGGTGTTGGTCTTGATGCCGTc  >  4:278914/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGA‑TTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGTC  >  NC_002947/648142‑648435

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