Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,784,030	G→A	7.1%	A16V (GCG→GTG)	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,784,030	0	G	A	7.1%	249.9 / 4.8	88	A16V (GCG→GTG)	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (39/40);  new base A (2/4);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

TTGCAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGC  >  NC_002947/1783882‑1784173                                                                                                                                                      |                                                                                                                                                 ttGCAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑g                                                                                                                                                  >  2:344550/1‑149 (MQ=255)   gCAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGACTTGAGGGCACGTTCC‑Acc                                                                                                                                                >  8:323468/1‑149 (MQ=255)    cAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCt                                                                                                                                               <  4:399268/149‑1 (MQ=255)       gATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTgg                                                                                                                                             >  6:347275/1‑148 (MQ=255)             gTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTgg                                                                                                                                       <  1:344550/148‑1 (MQ=255)              tCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGc                                                                                                                                      >  7:193921/1‑148 (MQ=255)              tCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGc                                                                                                                                      >  2:380465/1‑148 (MQ=255)              tCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑AGGCGGAACTGGcc                                                                                                                                     <  6:88642/149‑1 (MQ=255)              tCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑AGGCGGAACTGGcc                                                                                                                                     <  6:2652/149‑1 (MQ=255)              tCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑AGGCGGAACTGGc                                                                                                                                      <  5:204043/148‑1 (MQ=255)                     gCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTc                                                                                                                               >  6:396887/1‑148 (MQ=255)                       cAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCGGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGGGATGGGGCTCGCGTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGaa                                                                                                                            <  5:287033/149‑1 (MQ=255)                        aGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGTGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGaa                                                                                                                            >  6:40311/1‑148 (MQ=255)                             cgGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑AGGCGGAACTGGcc                                                                                                                                     <  7:137304/134‑1 (MQ=255)                             cgGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑AGGCGGAACTGGcc                                                                                                                                     >  8:137304/1‑134 (MQ=255)                              gggCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCg                                                                                                                     <  8:193921/149‑1 (MQ=255)                                gCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCggg                                                                                                                   <  2:145677/149‑1 (MQ=255)                                 ccTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCgggg                                                                                                                  <  5:396887/149‑1 (MQ=255)                                  cTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑t                                                                                                                <  3:295218/149‑1 (MQ=255)                                  cTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑t                                                                                                                <  3:289802/149‑1 (MQ=255)                                  cTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑t                                                                                                                >  5:409659/1‑149 (MQ=255)                                        cAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCt                                                                                                          >  6:280766/1‑149 (MQ=255)                                          gCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTgg                                                                                                        <  7:323468/149‑1 (MQ=255)                                              tCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCa                                                                                                    >  3:350385/1‑149 (MQ=255)                                               cGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCAt                                                                                                   >  3:181779/1‑149 (MQ=255)                                                  aTGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATcgc                                                                                                >  6:121189/1‑149 (MQ=255)                                                     tCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCa                                                                                                                                    >  2:279711/1‑111 (MQ=255)                                                     tCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCa                                                                                                                                    <  1:279711/111‑1 (MQ=255)                                                     tCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGATTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACgg                                                                                                                                 >  8:35609/1‑114 (MQ=255)                                                        cTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCt                                                                                                                                               >  6:199655/1‑97 (MQ=255)                                                        cTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCt                                                                                                                                               <  5:199655/97‑1 (MQ=255)                                                        cTGCCGCCTCGAAAACGATTGACGAGCACCTCGCCGGCCAGGCGGATTGCTATCATGAAGGCTGCGATCGAGCTCGCCTAGAGGCCCAGTTCC‑CCCTGGAACTGGCAACGGTCGAACAGCTCGGGG‑TCCAACAGTGGAATCGCGTCAct                                                                                          >  8:351402/1‑148 (MQ=255)                                                            cGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTcc                                                                                      >  5:157348/1‑149 (MQ=255)                                                            cGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTcc                                                                                      <  5:347275/149‑1 (MQ=255)                                                                gcgAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGTGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATcg                                                                                                 >  5:261320/1‑134 (MQ=255)                                                                gcgAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGTGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATcg                                                                                                 <  6:261320/134‑1 (MQ=255)                                                                gcgAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTAGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTc                                                                                  >  7:151185/1‑149 (MQ=255)                                                                  gAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGt                                                                                >  1:101353/1‑149 (MQ=255)                                                                        gTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTc                                                                                  >  1:220213/1‑141 (MQ=255)                                                                        gTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTc                                                                                  <  2:220213/141‑1 (MQ=255)                                                                        gTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTg                                                                          >  7:338415/1‑149 (MQ=255)                                                                        gTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAACGCGGCGATGGGGCACGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCACCCTGGGGCATCGCGTCACCTTCCTGTCGTGGGGTg                                                                          >  5:38992/1‑149 (MQ=255)                                                                         tttGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTgg                                                                         >  2:183589/1‑149 (MQ=255)                                                                            gTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCt                                                                      >  5:89843/1‑149 (MQ=255)                                                                             ccGAGCACGTCGCCGGCCAGGCGGATGTCTTTCTTGAAGGCGGCGATGGGGCTCGGCTTGAGGGCCAGTTGG‑GCCTGGAAGTGGCGACGGTCGAACAGGTGGGGG‑TCCACCTGGGGCATCGCGTCAGGTTCCTGTCGTGGGGTGGCCTg                                                                     <  6:38992/148‑1 (MQ=255)                                                                               gAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTc                                                                   >  5:241918/1‑149 (MQ=255)                                                                               gAGCACCTCGCCGGCCAGGCGGATGGCTGTCTTGAAGGCGGCGATGGGGCGCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTc                                                                   <  2:184076/149‑1 (MQ=255)                                                                                 gCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCgg                                                                 <  1:380465/149‑1 (MQ=255)                                                                                  gagctctCCGGCCAGGCGGATGGCTTTATTGAAGGGGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGTACTGGCCACGGTCGAACTGCTCGGGG‑TGCAGCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCggg                                                                <  2:145060/146‑1 (MQ=255)                                                                                   aCCTCGCCGGCCAGGCGGATGGGTTTCTTGAAGGCGGCGATTGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCggg                                                                <  1:300019/149‑1 (MQ=255)                                                                                    ccTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGccgttt                                                                   <  8:318291/144‑5 (MQ=255)                                                                                    ccTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGccgttt                                                                   >  7:318291/1‑140 (MQ=255)                                                                                     cTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTg                                                             >  2:371595/1‑149 (MQ=255)                                                                                     cTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTg                                                             >  3:111532/1‑149 (MQ=255)                                                                                      gccccGGCCAGGCGTATGGCTTTCTTGAAGGCGGCGATGGGGCTGGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑GCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGc                                                            <  7:167927/146‑1 (MQ=255)                                                                                        gCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCTGTTCC‑TCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCg                                                           <  7:41035/148‑1 (MQ=255)                                                                                        gCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCg                                                           <  5:42556/148‑1 (MQ=255)                                                                                        gCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGGTCC‑GCCTGGAACTGGCGACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGGTCCTGTCGTGGGGTGGCCTGTCGGGGTGCgg                                                          >  5:332645/1‑149 (MQ=255)                                                                                             ccAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAgg                                                     <  6:217479/149‑1 (MQ=255)                                                                                                  cGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGaa                                                <  2:101353/149‑1 (MQ=255)                                                                                                  cGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGaa                                                <  6:409659/149‑1 (MQ=255)                                                                                                    gATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGt                                                              <  4:234415/133‑1 (MQ=255)                                                                                                    gATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGt                                                              >  3:234415/1‑133 (MQ=255)                                                                                                     aTGGCTTTCTTGACGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTg                                             <  5:280766/149‑1 (MQ=255)                                                                                                         cTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGcgcg                                         >  1:168360/1‑149 (MQ=255)                                                                                                            tCTTGCAGGCGGCGATGGGGCTCGCCTTGAGGGGCAGGTCCTTTCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCgg                                        <  1:339885/148‑1 (MQ=255)                                                                                                                gAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGgt                                  >  2:298167/1‑149 (MQ=255)                                                                                                                 aaGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGgt                                  <  8:120299/148‑1 (MQ=255)                                                                                                                  aGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGgtg                                 >  4:320622/1‑148 (MQ=255)                                                                                                                  aGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGgtg                                 <  3:320622/148‑1 (MQ=255)                                                                                                                   ggcggcGTTTGGGGTCGGCTTGTGGGCCATTTCC‑GGCTGGAACTGGCCACGGGCGTACATGTGGGGG‑TCCACCTGGGGGATCGCGTCACGTTCCTGTCGTGGGGTGGGCTGTCGGGGTGGGGTCAGGGCGTAGTGCGCGGGATGGTgtc                               <  6:332645/149‑1 (MQ=255)                                                                                                                   ggcggcGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTgtc                               <  8:151185/149‑1 (MQ=255)                                                                                                                          aTGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTg                        <  6:89843/149‑1 (MQ=255)                                                                                                                           aggtgCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTgc                       <  5:62571/145‑1 (MQ=255)                                                                                                                              ggCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTgcg                      >  2:99292/1‑147 (MQ=255)                                                                                                                               gCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACATGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAg                   >  4:212597/1‑149 (MQ=255)                                                                                                                                 tCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAggg                 >  1:65726/1‑149 (MQ=255)                                                                                                                                 tCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAggg                 >  6:250836/1‑149 (MQ=255)                                                                                                                                   gCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTg               <  8:338415/149‑1 (MQ=255)                                                                                                                                   gCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGAGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGt                >  7:257812/1‑148 (MQ=255)                                                                                                                                    ccTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTg               >  7:303749/1‑148 (MQ=255)                                                                                                                                    ccTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGa              <  1:298167/149‑1 (MQ=255)                                                                                                                                     cTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGaa             <  4:181779/149‑1 (MQ=255)                                                                                                                                        gAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAg            <  4:350385/147‑1 (MQ=255)                                                                                                                                            gCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCg      >  4:26412/1‑149 (MQ=255)                                                                                                                                            gCCAGTTCC‑GCCTGGAACTGGCCACGGGCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGGGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCg      <  6:296123/149‑1 (MQ=255)                                                                                                                                              cAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTa    <  1:311138/149‑1 (MQ=255)                                                                                                                                              cAGTTAC‑GCCAGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCAACCTGGGGCATCGCGTCACGTTACTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGATGTGCGCGGGATGGTGTCGGCCTTGCTCAGGGTGAAGATCTCg      >  8:369991/1‑147 (MQ=255)                                                                                                                                                gTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGc  <  2:168360/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                 TTGCAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCC‑GCCTGGAACTGGCCACGGTCGAACAGCTCGGGG‑TCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGC  >  NC_002947/1783882‑1784173

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