Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,297,462	T→G	7.1%	Q155P (CAG→CCG)	PP_4671 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,297,462	0	T	G	7.1%	265.6 / 9.6	85	Q155P (CAG→CCG)	PP_4671	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (41/38);  new base G (2/4);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.55e-01

GCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTAGAGGGCA  >  NC_002947/5297314‑5297602                                                                                                                                                     |                                                                                                                                             gCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACctgct                                                                                                                                              <  7:201800/149‑1 (MQ=255)   aCCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTgg                                                                                                                                            <  5:267750/149‑1 (MQ=255)    ccGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGt                                                                                                                                           >  5:369271/1‑149 (MQ=255)    ccGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGt                                                                                                                                           >  7:13366/1‑149 (MQ=255)     cGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTc                                                                                                                                          >  5:205312/1‑149 (MQ=255)        gcgcAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTcgcg                                                                                                                                       <  8:159827/149‑1 (MQ=255)          gcAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGa                                                                                                                                     <  4:55407/149‑1 (MQ=255)          gcAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGa                                                                                                                                     >  6:15306/1‑149 (MQ=255)               gTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCg                                                                                                                                 >  5:252005/1‑148 (MQ=255)                  gggCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCgcagg                                                                                                                             >  2:223808/1‑149 (MQ=255)                  gggCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcagg                                                                                                                             <  1:71545/149‑1 (MQ=255)                    gcgcACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCgcaggca                                                                                                                           <  5:124752/149‑1 (MQ=255)                        aCGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAAc                                                                                                                       <  7:90839/149‑1 (MQ=255)                        aCGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcaggtgacct                                                                                                                       <  8:131444/149‑7 (MQ=255)                              cGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTg                                                                                                                 <  3:98248/149‑1 (MQ=255)                               gCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGc                                                                                                                >  6:52628/1‑149 (MQ=255)                                 ttGCGGGTC‑GCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGa                                                                                                             >  5:40306/1‑149 (MQ=255)                                    cGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGAcc                                                                                                           <  6:369271/149‑1 (MQ=255)                                      ggTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGt                                                                                                         <  2:340828/149‑1 (MQ=255)                                        tCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTg                                                                                                        >  6:129849/1‑148 (MQ=255)                                            cTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAgg                                                                                                                    <  6:43950/132‑1 (MQ=255)                                            cTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAgg                                                                                                                    >  5:43950/1‑132 (MQ=255)                                            cTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATgg                                                                                                   <  6:40306/149‑1 (MQ=255)                                            cTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATgg                                                                                                   <  1:223808/149‑1 (MQ=255)                                            cTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATgg                                                                                                   <  6:252005/149‑1 (MQ=255)                                                aTGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTACCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCa                                                                                               >  5:62210/1‑149 (MQ=255)                                                      gctCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGc                                                                                          >  2:189237/1‑148 (MQ=255)                                                      gctCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcagg                                                                                                                             >  5:367316/1‑113 (MQ=255)                                                      gctCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcagg                                                                                                                             >  5:362078/1‑113 (MQ=255)                                                      gctCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcagg                                                                                                                             <  6:362078/113‑1 (MQ=255)                                                      gctCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcagg                                                                                                                             <  6:367316/113‑1 (MQ=255)                                                        tCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCAGcgag                                                                                       >  5:76276/1‑146 (MQ=255)                                                          cGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGc                                                                                     <  4:134900/149‑1 (MQ=255)                                                          cGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGc                                                                                     <  2:267547/149‑1 (MQ=255)                                                          cGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGc                                                                                     >  7:51874/1‑149 (MQ=255)                                                               ggCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATg                                                                                <  5:52628/149‑1 (MQ=255)                                                                  gccgcGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCt                                                                              >  3:276124/1‑148 (MQ=255)                                                                  gccgcGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTc                                                                             <  6:171720/149‑1 (MQ=255)                                                                   ccgcGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCa                                                                            <  4:353808/149‑1 (MQ=255)                                                                       gACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCAc                                                                        >  5:74253/1‑149 (MQ=255)                                                                            ttGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGt                                                                   <  8:13366/149‑1 (MQ=255)                                                                                agcaCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGgcc                                                               <  6:320672/149‑1 (MQ=255)                                                                                 gcaCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTc                                                                             >  3:60449/1‑134 (MQ=255)                                                                                 gcaCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTc                                                                             <  4:60449/134‑1 (MQ=255)                                                                                 gcaCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGgccg                                                              >  4:151736/1‑149 (MQ=255)                                                                                  caCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCACAGCACCGGGTGgccgc                                                             >  5:305333/1‑149 (MQ=255)                                                                                   aCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGgccgcc                                                            <  8:51874/149‑1 (MQ=255)                                                                                     gCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCaa                                                          <  6:208112/149‑1 (MQ=255)                                                                                     gCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCaa                                                          <  6:213236/149‑1 (MQ=255)                                                                                     gCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCaa                                                          >  8:98927/1‑149 (MQ=255)                                                                                       tGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATa                                                        <  6:74253/149‑1 (MQ=255)                                                                                         gggTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGATGATGCTCAGCACCGGGTGGCCGCCAATACg                                                      <  3:52133/149‑1 (MQ=255)                                                                                           gTGAATGCGCACGCCAGAGCGTGCCCCCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGCCCGGGCATGGGTCATCCGCCTAGCTGATGCTCAGCACCGGGTGGCCTCCACTACGCt                                                    >  8:79097/1‑149 (MQ=255)                                                                                           gTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCt                                                    >  4:234366/1‑149 (MQ=255)                                                                                            tGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTg                                                   <  3:151736/149‑1 (MQ=255)                                                                                                   gcACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGc                                            >  7:124427/1‑149 (MQ=255)                                                                                                    cACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGcc                                           <  3:208747/149‑1 (MQ=255)                                                                                                    cACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGcc                                           <  2:313043/149‑1 (MQ=255)                                                                                                      cGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGa                                         <  5:15306/149‑1 (MQ=255)                                                                                                      cGCCAGAGCGTACCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGa                                         <  6:62210/149‑1 (MQ=255)                                                                                                         cAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATcg                                      >  6:333990/1‑149 (MQ=255)                                                                                                          agagCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGCCCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATcgc                                     >  4:156291/1‑149 (MQ=255)                                                                                                           gagCGTGCCCCCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATcgcg                                    >  1:163117/1‑149 (MQ=255)                                                                                                              cGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCAGCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTg                                 <  7:79097/149‑1 (MQ=255)                                                                                                                tGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGc                               >  4:131786/1‑149 (MQ=255)                                                                                                                  ccACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCAACCGTGCATGGGCCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGccc                             >  4:90506/1‑149 (MQ=255)                                                                                                                   cACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCt                            >  7:285108/1‑149 (MQ=255)                                                                                                                              gaTTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCTGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTg                 <  7:314789/149‑1 (MQ=255)                                                                                                                                  gACTTCAAAGGTCACCTGCTGGTCGCGGACCCgcaggc                                                                                                                            >  2:236687/1‑38 (MQ=255)                                                                                                                                  gACTTCAAAGGTCACCTGCTGGTCGCGGACCCgcaggc                                                                                                                            <  1:236687/38‑1 (MQ=255)                                                                                                                                  gACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCaa             >  8:82612/1‑149 (MQ=255)                                                                                                                                    cTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATg           >  1:120344/1‑149 (MQ=255)                                                                                                                                     ttCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCTTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGGCAATg           >  7:171400/1‑148 (MQ=255)                                                                                                                                     ttCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATg           >  5:33470/1‑148 (MQ=255)                                                                                                                                     ttCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGAAGAACAGGCTGCAGACCGTGCATGGGTCATACGCCTCGCTGATGATCAGCACCGTGTGGCCGCCAATACGCTGGAACAGCCGATAGCGCTGGCCCTCGTCAAGCTTGCCAAtc           >  4:238422/1‑147 (MQ=255)                                                                                                                                      tCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGGCAATg           >  3:242810/1‑147 (MQ=255)                                                                                                                                       cAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGt          <  3:238422/147‑1 (MQ=255)                                                                                                                                        aaaGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTaga       <  2:120344/149‑1 (MQ=255)                                                                                                                                        aaaGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTa         >  5:275699/1‑147 (MQ=255)                                                                                                                                         aaGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTagag      <  3:353608/149‑1 (MQ=255)                                                                                                                                         aaGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTagag      <  6:33470/149‑1 (MQ=255)                                                                                                                                            gTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTAGAggg    >  7:248779/1‑148 (MQ=255)                                                                                                                                             tCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCTCTGGTCAAGCTTGCCAATGTAGAGGGCa  >  1:359187/1‑149 (MQ=255)                                                                                                                                             tCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTAGAGGGc   >  4:81272/1‑148 (MQ=255)                                                                                                                                               aCCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCt                                                                                    <  8:150774/65‑1 (MQ=255)                                                                                                                                               aCCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCt                                                                                    >  7:150774/1‑65 (MQ=255)                                                                                                                                                  tgctgGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTgg                                >  5:120290/1‑114 (MQ=255)                                                                                                                                                  tgctgGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTgg                                <  6:120290/114‑1 (MQ=255)                                                                                                                                                    ctgGTCGCGGACCCGCAGGCAGAACAGGCTTCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGTCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGcc               <  3:367777/129‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                             GCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTAGAGGGCA  >  NC_002947/5297314‑5297602

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