Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,740,811	(T)7→6	intergenic (+48/+705)	gltA → / ← yeiW	citrate synthase/putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,740,805	0	T	.	90.2%	178.8 / 18.1	52	intergenic (+42/+711)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base . (22/24);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 3.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGC  >  NC_002947/4740659‑4740949                                                                                                                                                   |                                                                                                                                                   cACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttt                                                                                                                                                  <  2:73983/149‑1 (MQ=255)      gaaagaaaTGCTCTCCAGACCGTACAAGATCGGCCGTCCGCGCCAGCTGGACACCGCCGAGCAGAAGCGCGACATCGTAGCCCTGAAAGACCGCAAGAAAGCAGCACTAGCCGAACGCAAGAATGCTCCCCTGGCGCAGCCttttttat                                                                                                                                             >  1:291831/1‑147 (MQ=255)        aagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                                           <  7:269251/149‑5 (MQ=255)           aaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttgtc                                                                                                                                        <  7:35679/149‑8 (MQ=255)                     cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGc                                                                                                                             <  5:227523/149‑1 (MQ=255)                     cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGc                                                                                                                             >  2:221320/1‑149 (MQ=255)                              cAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                                <  3:218068/122‑1 (MQ=255)                              cAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                                >  4:218068/1‑122 (MQ=255)                              cAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                                     <  4:429206/148‑1 (MQ=255)                                 gATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGcc                                                                                                                 >  5:146261/1‑149 (MQ=255)                                     ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTg                                                                                                              <  1:253223/148‑1 (MQ=255)                                       ccGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGccc                                                                                                           <  6:272306/149‑1 (MQ=255)                                            cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGTCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                      <  8:16202/149‑1 (MQ=255)                                            cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                      <  3:215429/149‑1 (MQ=255)                                            cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                      <  6:185386/149‑1 (MQ=255)                                                  aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGcc                                                                                                <  3:216245/149‑1 (MQ=255)                                                   gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCg                                                                                               <  4:106205/149‑1 (MQ=38)                                                   gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCg                                                                                               <  5:425283/149‑1 (MQ=38)                                                           ccGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCGCTTGGCGGGCGTGCCCGCGAAGAGGGCGGGacaca‑‑c                                                                                     >  4:255028/1‑149 (MQ=255)                                                            cGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACTTAGCGGCCTTGCCCGCGAAGAGACCGGCACAGa‑‑ca                                                                                    >  8:202682/1‑149 (MQ=255)                                                                      aaGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAAcc                                                                          >  4:309390/1‑149 (MQ=255)                                                                             aCATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGttt                                                                   >  6:185909/1‑149 (MQ=255)                                                                                   ttGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTTCT‑GTGCCGGCCTCATCGCGGCCTTGCCCGCGAAGAGGCCGGCGAAGA‑‑CAACACAAAACCTCAGgttttatac                                                             >  5:270796/1‑145 (MQ=255)                                                                                    tGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTcc                                                             <  5:185909/148‑1 (MQ=255)                                                                                      cccTCAAAGACCCCAAGTACGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTccgcc                                                          <  6:146261/149‑1 (MQ=255)                                                                                       ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTATccgccc                                                         >  5:257186/1‑149 (MQ=255)                                                                                       ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTATccgccc                                                         >  7:249969/1‑149 (MQ=255)                                                                                        cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACTTAGCGGCCGTGGCCGCGCAGAGGCCGGCCCAGC‑‑CACCACCAAACCGAAGTTTTTGTacgcccc                                                        >  8:223419/1‑148 (MQ=255)                                                                                         tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTccgcccgc                                                       >  6:345434/1‑149 (MQ=255)                                                                                          gAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTGTTTCTCCGCCCGCt                                                      >  5:218562/1‑149 (MQ=255)                                                                                           aaaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGACGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACCGA‑‑CACCACAAAACCTCAGTTTTTCTccgcccgca                                                      >  5:163103/1‑147 (MQ=255)                                                                                            aaGACCGCAAGTAAGCTGCACTAGCCGAAAGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTAGTCCTGTGCCGCCCTCATCGCGGACTTCCCCCCGAAGAGGCCGGCACACA‑‑AACCACAAAACCTAACTTATAATCCGCCCGCTCa                                                    >  4:276714/1‑149 (MQ=255)                                                                                                 cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGcc                                                                                                <  3:130411/102‑1 (MQ=39)                                                                                                 cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGcc                                                                                                >  4:130411/1‑102 (MQ=39)                                                                                                 cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCGCAGGTTTTGTACGCCCGCGCACGCAg                                               >  7:236334/1‑149 (MQ=255)                                                                                                     aGTATGCTGCTGTTGCCTAACGCAAGAAGGCTGCCTTGGGTCAG‑CTTTTTTTATTTGTCTTGTGCCGGGCTCTTCGCTGCCGTGCCCGCGATGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGACCGCTCACGCAGCCCt                                           <  7:202682/149‑1 (MQ=255)                                                                                                     aGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCt                                           <  3:309390/149‑1 (MQ=255)                                                                                                       tAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGC‑‑CACCACAAAACCTCAGGTTTTCTCCGCCCGCTCACCCAGcccat                                          >  5:461182/1‑146 (MQ=255)                                                                                                         cgCTGCAGTAGCCCAACTCAAGACGGGGGCCGGCGGGC‑TTCGTTTTTTATTTGTCGTGGGCCGGCCTCTTCGCGGCCTGGCCGGGGATGTGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                       <  6:163103/148‑1 (MQ=255)                                                                                                          gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCATGCCCGCGAAGAGGCCGGCACAGAGTTGGTGCAAAACAACAGATTTTCTCAGCCCGCTCACTCACCCCtat                                         >  1:387741/1‑146 (MQ=255)                                                                                                           cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                     <  5:345434/149‑1 (MQ=255)                                                                                                           cTGCAGTACCCTAACGCAAGATGGCTTCCCTCGGGCAGCGTTTTTTT‑TTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                     <  2:387741/149‑1 (MQ=255)                                                                                                                 tAGCCTACCGCAATAAGGGTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                                <  8:249969/148‑1 (MQ=255)                                                                                                                       accGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTGGCCCGCTAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                         <  7:223419/147‑1 (MQ=255)                                                                                                                       aaCGCAAGACGGCGGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTCCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                         <  1:276732/149‑1 (MQ=255)                                                                                                                       aaCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGCCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                         <  1:243057/149‑1 (MQ=255)                                                                                                                            aagaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCCCGCGAACAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc                    >  7:31201/1‑149 (MQ=255)                                                                                                                                 ggCTGCCCTCGTGCAGCC‑TTTTTTATTTGTCCGTTCCCGGCCCCTGCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacca               <  2:158104/149‑1 (MQ=255)                                                                                                                                 ggCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                     >  2:56316/1‑127 (MQ=255)                                                                                                                                 ggCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                     <  1:56316/127‑1 (MQ=255)                                                                                                                                  gCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCATGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCACAGTTTTTATACCCCCCCTAACCCAGCCCGTTAACGGTATTGAACCGCGCAaccacccc              >  8:114869/1‑147 (MQ=255)                                                                                                                                      cccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc          <  7:114869/149‑1 (MQ=255)                                                                                                                                      cccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc          >  6:130747/1‑149 (MQ=255)                                                                                                                                           gggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttg       <  5:130747/147‑1 (MQ=255)                                                                                                                                           ggccAGCC‑TTTTTTATTTGTCCTGTCCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt     <  2:328566/146‑1 (MQ=255)                                                                                                                                            ggccagcTTTTTTTATTTGTCGTGTGCCGGCCTCTTCGCTGCCTTGCCCGGGAAGTGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCAAGCAGCCCTTTCAGGGTATTGAACGGCTCATCCACCACGAACttg       <  5:234057/141‑1 (MQ=255)                                                                                                                                            ggCAGCTTTTTTTTTTTGCCCTGG‑GCCGCCCTCTTCGCGCCGTGGCCCGCTAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg    <  3:255028/149‑1 (MQ=255)                                                                                                                                             gCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCAAAGTTTTTCTCCGCCCGCTCCCGCAGCCCATTCAGGGTATTGAACGGCGCATCCACCACCAACTTGTTgg   >  8:10567/1‑149 (MQ=255)                                                                                                                                               ttggttttttttttttTCCTGTGCGGCCCTTTTCGGGGCCTTGCCCGCTAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGc  <  7:173984/145‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGC  >  NC_002947/4740659‑4740949

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