Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,212,492	C→T	6.8%	intergenic (+56/+38)	PP_4593 → / ← PP_4594	hypothetical protein/cystathionine gamma‑synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,212,492	0	C	T	6.8%	215.3 / 4.5	73	intergenic (+56/+38)	PP_4593/PP_4594	hypothetical protein/cystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base C (38/30);  new base T (2/3);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCAT  >  NC_002947/5212344‑5212634                                                                                                                                                         |                                                                                                                                               ggCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACcgtc                                                                                                                                                <  7:399707/149‑1 (MQ=255)  gCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACcgtcg                                                                                                                                               <  7:284359/149‑1 (MQ=255)     ccAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAg                                                                                                                                            <  3:93188/149‑1 (MQ=255)      cAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAg                                                                                                                                            >  5:75303/1‑148 (MQ=255)       aaCGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAg                                                                                                                                            >  6:397916/1‑147 (MQ=255)           acgaAGCCGTCAACCGCCTGCGCAAGGAACGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCgg                                                                                                                                      >  8:289208/1‑149 (MQ=255)           acgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑CTTTCGCGGGTAAACC‑CGCTCCTACGACGGTGTAGGAGCg                                                                                                                                       <  4:139776/149‑1 (MQ=38)           acgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACGACGGTGTAGGAGCgg                                                                                                                                      <  7:298211/149‑1 (MQ=255)             gaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGt                                                                                                                                    <  2:128288/149‑1 (MQ=255)                   gTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTAccc                                                                                                                              >  4:17680/1‑149 (MQ=255)                   gTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTAccc                                                                                                                              <  5:177689/149‑1 (MQ=255)                   gTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGACCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTAccc                                                                                                                              >  1:148865/1‑149 (MQ=255)                      aaCCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCcgc                                                                                                                            <  7:97050/148‑1 (MQ=255)                             tGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGc                                                                                                                    >  2:91250/1‑149 (MQ=255)                               cgcAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCt                                                                                                                  >  8:146820/1‑149 (MQ=255)                                  aaGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGCCAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGTAAACC‑CGCTCCTACGACGGTGTAGGAGCGGGCTTGCCcgc                                                                                                                            >  5:433808/1‑136 (MQ=25)                                  aaGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGTAAACC‑CGCTCCTACGACGGTGTAGGAGCGGGCTTGCCcgc                                                                                                                            <  6:433808/136‑1 (MQ=37)                                    ggAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGc                                                                                                             >  5:89588/1‑149 (MQ=255)                                      aaaGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGccc                                                                                                           >  2:303536/1‑149 (MQ=255)                                             gCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAGAGGGCCGCAGGCCCtttcgc                                                                                                     <  3:27059/148‑6 (MQ=38)                                             gCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGCAGGCCCtttcgc                                                                                                     <  1:207785/148‑6 (MQ=255)                                                 aTCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGa                                                                                                >  2:29785/1‑149 (MQ=255)                                                  tCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAc                                                                                               <  2:251486/149‑1 (MQ=255)                                                  tCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAc                                                                                               >  2:342305/1‑149 (MQ=255)                                                    aaGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCa                                                                                             <  6:75303/149‑1 (MQ=255)                                                       gACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgcg                                                                                          >  4:389745/1‑149 (MQ=255)                                                       gACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgcg                                                                                          >  4:206382/1‑149 (MQ=255)                                                        aCAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgcgc                                                                                         <  1:29785/149‑1 (MQ=255)                                                          aaGATCGAACGGATCATTCATCCGCCGAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCtt                                                                                       >  1:12516/1‑149 (MQ=255)                                                            gATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCTGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCatt                                                                                      >  7:367055/1‑148 (MQ=255)                                                              ttgtaCGGTTTATTCCTCCGCCCAAATGCTTTCCTTCATTTTTATGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTc                                                                                     <  1:147522/144‑1 (MQ=25)                                                                gAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAAc                                                                                  <  5:343017/148‑1 (MQ=255)                                                                 aaCGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGc                                                                                >  2:109961/1‑149 (MQ=255)                                                                   cGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGccc                                                                              <  7:289208/149‑1 (MQ=255)                                                                    ggATCATTCATCCGCCCAGATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGcccc                                                                             <  5:93751/149‑1 (MQ=255)                                                                       tcattcatCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCTAA‑TACGCGGGCAAGCC‑CGATCCTAAACCGACGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCACACACAGCGACCAGCGCTTTCAACGCCCCAt                                                                           >  1:103106/1‑149 (MQ=255)                                                                       tcattcatCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATc                                                                          >  1:419868/1‑149 (MQ=255)                                                                         attcatCCGCCCAAATCCTGCCAGGGATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGTGAGGGTTTAGCCGCGAAAGGGCCTGCGGCCCTCACACGGCGACCAGCGCTTTCAAGGCCCCATc                                                                          >  6:362627/1‑147 (MQ=255)                                                                             atCCGCCCAAATCCTAACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTTAACGCCCCATCCAGGTc                                                                    >  8:300409/1‑149 (MQ=255)                                                                               ccGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCt                                                                  >  2:116022/1‑149 (MQ=255)                                                                                cGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCt                                                                  <  2:12516/148‑1 (MQ=255)                                                                                  cccAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGa                                                               >  2:165590/1‑149 (MQ=255)                                                                                   ccAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCCGCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGa                                                               >  1:108789/1‑149 (MQ=255)                                                                                   ccAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCCGCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGa                                                               >  1:114638/1‑149 (MQ=255)                                                                                     aaaTCCTGACAGGCATCTAT‑TGAATCAACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGTGCCTGCGACACTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGAc                                                              >  1:131908/1‑147 (MQ=255)                                                                                             aCAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACgg                                                    <  4:11788/149‑1 (MQ=255)                                                                                             aCAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACgg                                                    <  1:277815/149‑1 (MQ=255)                                                                                                 gCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATc                                                                          >  5:409654/1‑123 (MQ=255)                                                                                                 gCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATc                                                                          <  6:409654/123‑1 (MQ=255)                                                                                                 gCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAg                                                           <  5:393896/138‑1 (MQ=255)                                                                                                 gCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAg                                                           >  6:393896/1‑138 (MQ=255)                                                                                                 gCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAAt                                                >  7:186156/1‑149 (MQ=255)                                                                                                       aT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCCGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGa                                          >  8:108706/1‑149 (MQ=255)                                                                                                        t‑tGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGa                                          <  3:17680/148‑1 (MQ=255)                                                                                                         ttcaaTCCACCGGCCT‑C‑TTCGCGTGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGa                                          <  7:192566/145‑1 (MQ=255)                                                                                                              tCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTa                                    >  6:107065/1‑149 (MQ=255)                                                                                                               ccACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTca                                                                                                        <  5:148658/80‑1 (MQ=255)                                                                                                                cACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAATGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTAc                                   <  7:357389/148‑1 (MQ=255)                                                                                                                cACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTAc                                   <  2:71164/148‑1 (MQ=255)                                                                                                                   cGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGAc                               <  2:148865/149‑1 (MQ=255)                                                                                                                        t‑c‑tTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCg                          >  4:226810/1‑149 (MQ=255)                                                                                                                              cgcgGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCgcag                       >  2:406552/1‑148 (MQ=255)                                                                                                                              cgcgGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCgcag                       <  1:303536/148‑1 (MQ=255)                                                                                                                                 gggCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAgg                   <  2:114638/149‑1 (MQ=255)                                                                                                                                 gggCAAGCC‑CGCTCCTACACCGGCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAgg                   <  2:108789/149‑1 (MQ=255)                                                                                                                                        cc‑cGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGc            <  3:206382/149‑1 (MQ=255)                                                                                                                                         c‑cGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCt           >  5:8703/1‑149 (MQ=255)                                                                                                                                           cGCTCCTACGACGGTGTAGGAGCGGGCTTGCCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCTGCAGGTTCTCGcg           >  4:347946/1‑147 (MQ=38)                                                                                                                                            gCTCCTACACCGTCGTAGGAGAGGTATTACCAGCGAAAGGGCCTGCGGACCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGTTCACGGTAATCCTCCATCCCTACCGACAAACGCACCAGGTTCACGATg          >  4:134311/1‑148 (MQ=255)                                                                                                                                              tCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATc       >  3:222470/1‑149 (MQ=255)                                                                                                                                                cTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATccc     >  4:65793/1‑149 (MQ=255)                                                                                                                                                cTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGACACTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGTTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTAGCTGATccc     <  2:131908/149‑1 (MQ=255)                                                                                                                                                 tACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATccc     >  4:233484/1‑148 (MQ=255)                                                                                                                                                   cacCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACATCCGCAGCAGGTTCTCGCTGATCGCCa   >  4:308286/1‑148 (MQ=255)                                                                                                                                                   cacCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCAt  <  6:8703/149‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                               GGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTAT‑TGAATCCACCGGCCT‑C‑TTCGCGGGCAAGCC‑CGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCAT  >  NC_002947/5212344‑5212634

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