Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	648,287	A→C	13.7%	intergenic (+56/+35)	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,287	0	A	C	13.7%	165.7 / 5.1	61	intergenic (+56/+35)	acoR/accC	acetoin catabolism regulatory protein/acetyl‑CoA carboxylase biotin carboxylase subunit
Reads supporting (aligned to +/- strand):  ref base A (14/38);  new base C (7/2);  total (21/40)
Fisher's exact test for biased strand distribution p-value = 5.72e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGC  >  NC_002947/648147‑648431                                                                                                                                             |                                                                                                                                                   ggCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCTGCGGATGAAGGGGTAGAGCCCTGGGGGCGCGACGCGGCCCCAGCTTCTGGGAGGCAATAAAAAAACCCCGCCCAC‑GGccc                                                                                                                                           >  4:132050/1‑148 (MQ=255)    aCGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCGGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCgggg                                                                                                                                        <  3:132050/149‑1 (MQ=255)    aCGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCgggg                                                                                                                                        <  3:199545/149‑1 (MQ=255)       ggCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTc                                                                                                                                     <  7:273652/149‑1 (MQ=255)        gCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGTATGAAGGGGTAGAGCCCTGTGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCCCTG‑TGCCGGGGtttt                                                                                                                                    <  7:246239/149‑3 (MQ=255)            tgcgattttCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCAGCGACCTGCCGCGGATTAAGGGGTAGAGCCCGGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGt                                                                                                                                 <  8:323792/142‑1 (MQ=255)            tGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCAGCGACCTGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTg                                                                                                                                   <  7:324161/146‑1 (MQ=255)                    tCGCGGGCTACGCTTTAGCGGCGAGTGCAGCACCACCGAATCCACATCCCGCGTTTGATGGGTTAGAGCCCGGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATcaca                                                                                                                         <  2:188024/148‑1 (MQ=255)                      gcgGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGCCGCGGCCCCAGGGGCTGGCAGGCCAAAAAAAAACCCCGCCCAG‑GGCCCGGGGTCAGGTAGAAAACACt                                                                                                                       >  5:96779/1‑148 (MQ=255)                            gACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGAt                                                                                                               <  1:82375/149‑1 (MQ=255)                              cGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCa                                                                                                             <  7:319725/149‑1 (MQ=255)                                   tATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTc                                                                                                          <  4:151280/147‑1 (MQ=255)                                     tCGGCGTGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTg                                                                                                       <  7:128857/148‑1 (MQ=255)                                     tCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTg                                                                                                       <  8:64439/148‑1 (MQ=255)                                      cGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGt                                                                                                                                 <  5:276666/122‑1 (MQ=255)                                       ggCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGtt                                                                                                    <  2:11525/149‑1 (MQ=255)                                       ggCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGt                                                                                                     <  5:80592/148‑1 (MQ=255)                                        gCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTg                                                                                                   <  3:188570/149‑1 (MQ=255)                                        gCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACCCCGCCCCCAAATATGGCAGGCAATAAAAAAACCCCGCACAC‑GGGCCGGGTTCTGGTAAAACAAACAA‑AAAGATCACACCAGCTTg                                                                                                   >  4:109062/1‑149 (MQ=255)                                                atccgCACCGCATCGACATGCCGCGGATGACGGGGCAGCGCCCTGGGGCCGCGCCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGttt                                                                                           <  4:150475/145‑1 (MQ=255)                                                  gAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGGCGCGGCCCCAGCTTATGGCAGGCAATAAAAAACCCCGGCCCAG‑GGGCGGGGGTGTGGTACATCACCCTA‑CACGATGACTCCGGGGTGGGCAGTTtgg                                                                                         >  5:370327/1‑148 (MQ=255)                                                     cACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGCCGCGGGCCCCCCATGTGGGAGGCAATAAAAAACCCCCCCACAG‑GGCCGGGGGTGGGGTAAATCACACAA‑CACGAAGACGCCGGGGTGGCCCGGTtgttgg                                                                                      >  4:79063/1‑146 (MQ=255)                                                        cGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGa                                                                                   <  5:307648/149‑1 (MQ=255)                                                          aatgtcgaTGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAgg                                                                                 <  7:11340/142‑1 (MQ=255)                                                               aCATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTgg                                                                                                  <  3:334027/127‑1 (MQ=255)                                                                   gCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGCCGCGGCCCCACATTATGGGAGGCAAAAAAAAAACCCCGCCCAC‑GGCCGGGGTTGTGGTACATAACACAA‑CACCATCACACCGGGTTGGCCAGTTTGTGCTGGAGGTTGGGaaag                                                                        >  4:277899/1‑147 (MQ=255)                                                                    ccGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGt                                                                       <  3:109062/149‑1 (MQ=255)                                                                       cgGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGaa                                                                          <  5:281853/143‑1 (MQ=255)                                                                       cgGATGAAGGGGTAGAGCCCTGGGGCCCCGCCGCGGCCCCAGAATATGGCAGGCAATAAAAAAACCCCGCCCAG‑GGCCGGGGGTGTGGTAAAAAACACAA‑AACCAACACAACTGGGTTGGCACTTTGTGCTCGAGGTAGTGAAtagaac                                                                    >  6:281853/1‑144 (MQ=255)                                                                                 ggTAGAGCCCTGGGGCCGCGACGCGTCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCttt                                                          <  6:96779/149‑1 (MQ=255)                                                                                  gTAGAGCCCTGGGGCCGCGGCGCGGGCCCAGAATCTTGCAGGCAATAAAAAACCCCGGCCCCGCGGCCGGGG‑TGTGGTACATCACACAA‑CACGCTCACACCGGGGTGGCCGGGTTTTGGTCGGGGGTGGGGAAGGTGTCGCCGCCGTTg                                                         >  2:51332/1‑149 (MQ=255)                                                                                     gagCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCa                                                       <  5:125171/148‑1 (MQ=255)                                                                                         tcTGGGGCCGCGTCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGc                                                  >  1:287299/2‑149 (MQ=255)                                                                                           tGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCt                                                >  6:160221/1‑149 (MQ=255)                                                                                                     aCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCa                                      <  2:327831/149‑1 (MQ=255)                                                                                                        cgGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCa                                   <  5:263289/149‑1 (MQ=255)                                                                                                            cccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCtt                                               >  6:54867/1‑133 (MQ=255)                                                                                                              cAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCATTTGCAGAAGCCTTCATCACGCACCAGGTcgc                             >  8:204681/1‑149 (MQ=255)                                                                                                               aGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑TGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcgcg                            <  1:174069/149‑1 (MQ=255)                                                                                                               aGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcgcg                            >  5:125963/1‑149 (MQ=255)                                                                                                               aGATGCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcgcg                            <  6:125963/149‑1 (MQ=255)                                                                                                                gATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGTGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCgg                           <  8:219008/149‑1 (MQ=255)                                                                                                                 aTTCTGGCAGGCAATATAAAAGCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGTCCTGGTTGTTTTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTtgctg                           <  5:267486/148‑3 (MQ=255)                                                                                                                     tGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAg                      >  3:209763/1‑149 (MQ=255)                                                                                                                       gcaggcaATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGCAGCCTTCATCCCGCACCAGGTCGCGGTGCAGCg                    >  4:95408/1‑149 (MQ=255)                                                                                                                             aaTAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTAGCACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGCCGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATg               >  1:274430/1‑149 (MQ=255)                                                                                                                             aaTAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGaa                                                    >  6:318579/1‑111 (MQ=255)                                                                                                                             aaTAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGaa                                                    <  5:318579/111‑1 (MQ=255)                                                                                                                             aaTAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGt              <  6:370327/149‑1 (MQ=255)                                                                                                                               tAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGcc                                                 >  8:202903/1‑112 (MQ=255)                                                                                                                               tAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGcc                                                 <  7:202903/112‑1 (MQ=255)                                                                                                                                 aaaaaCCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGt          <  3:277899/149‑1 (MQ=255)                                                                                                                                    aaCCCCGGCCCTG‑TGCCGGGGTTTTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCtt       >  5:294533/1‑149 (MQ=255)                                                                                                                                    aaCCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCtt       <  2:287299/149‑1 (MQ=255)                                                                                                                                     aCCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGt          <  7:220535/145‑1 (MQ=255)                                                                                                                                     aCCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGt          >  8:220535/1‑145 (MQ=255)                                                                                                                                      ccccGGCACAG‑GGCCGGTGTTCTGGTACATCACTCTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGATGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGTGATGTTGGTCTTGa     <  5:264406/149‑1 (MQ=255)                                                                                                                                       cccGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGAt    <  7:10239/149‑1 (MQ=255)                                                                                                                                        ccGGCCCTG‑TGCCGGGGTTTTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATg   <  6:294533/149‑1 (MQ=255)                                                                                                                                         cGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTTCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTTCAGCGGGATGTTGGTCt        <  2:245847/143‑1 (MQ=255)                                                                                                                                         cGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCttagatc  >  1:245847/1‑144 (MQ=255)                                                                                                                                             |                                                                                                                                                   GGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAG‑GGCCGGGGTTCTGGTACATCACACTA‑CACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGC  >  NC_002947/648147‑648431

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑