Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	988,216	G→A	6.1%	R23H (CGC→CAC)	hisS →	histidine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	988,216	0	G	A	6.1%	243.9 / 2.4	81	R23H (CGC→CAC)	hisS	histidine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (39/37);  new base A (3/2);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.29e-01

CCAGAAAGCGGCCGAAAAGGCCGAGGCCGACGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCGGTGAAGTGACCGACATCGTCGAAAAAGAGATGTACACCT  >  NC_002947/988069‑988356                                                                                                                                                     |                                                                                                                                              ccAGAAAGCGGCCGAAAAGGCCGAGGCCGACGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGgcgc                                                                                                                                              >  6:310959/1‑149 (MQ=255)   aGAAAGCGGCCGAAAAGGCCGAGGCCGACGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTa                                                                                                                                            >  4:300962/1‑149 (MQ=255)   aGAAAGCGGCCGAAAAGGCCGAGGCCGACGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTa                                                                                                                                            >  5:83144/1‑149 (MQ=255)    gAAAGCGGCCGAAAAGGCCGAGGCCGACGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTAc                                                                                                                                           <  3:80002/149‑1 (MQ=255)     aaaGCGGCCGAAAAGGCCGAGGCCGACGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACt                                                                                                                                          >  4:206095/1‑149 (MQ=255)          ggCCGAAAAGTGCGAGGCCGACTCGGCGTTGATCGCTCGTTGCTGACACACATATTTCTTAAGTACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGaa                                                                                                                                     <  7:167922/149‑1 (MQ=255)              gAAAAGGCCGAGGCCGACGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGca                                                                                                                                 <  6:157044/149‑1 (MQ=255)                   ggccgaggccgaCGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTg                                                                                                                            <  5:175411/149‑1 (MQ=255)                    gccgaggccgaCGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTgg                                                                                                                           >  8:206938/1‑149 (MQ=255)                    gccgaggccgaCGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTgg                                                                                                                           >  5:83024/1‑149 (MQ=255)                    gccgaggccgaCGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTgg                                                                                                                           <  5:15965/149‑1 (MQ=255)                     ccgaggccgaCGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAAGGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGc                                                                                                                          >  1:356350/1‑149 (MQ=255)                      cgaggccgaCGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGcc                                                                                                                         >  4:342659/1‑149 (MQ=255)                        aggccgaCGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCgg                                                                                                                       >  6:280669/1‑149 (MQ=255)                          gccgaCGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGgc                                                                                                                     <  8:177389/149‑1 (MQ=255)                            cgaCGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGgctg                                                                                                                   >  7:267631/1‑149 (MQ=255)                              aCGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGgctgct                                                                                                                 >  5:197071/1‑149 (MQ=255)                                    cgTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGacac                                                                                                           <  3:300962/149‑1 (MQ=255)                                     gTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACAcc                                                                                                          <  3:428102/149‑1 (MQ=255)                                         aTCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑g                                                                                                     >  5:337064/1‑149 (MQ=255)                                                gTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACa                                                                                               <  2:356350/148‑1 (MQ=255)                                                     tGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGc                                                                                             >  1:30065/1‑145 (MQ=255)                                                       acacacAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGAt                                                                                        <  3:342659/148‑1 (MQ=255)                                                       acacacAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATc                                                                                       <  4:61083/149‑1 (MQ=255)                                                        cacacaGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATcc                                                                                      <  7:206938/149‑1 (MQ=255)                                                         acacaGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCg                                                                                     >  6:54656/1‑149 (MQ=255)                                                         acacaGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCg                                                                                     >  1:156160/1‑149 (MQ=255)                                                         acacaGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCg                                                                                     >  5:152860/1‑149 (MQ=255)                                                           acaGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTTCCAGAACCGTC‑GCCGCTGTGGCGCTACTGCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAG‑GGTTGCAGCCAGATCCGCa                                                                                   >  4:296773/1‑149 (MQ=255)                                                           acaGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑gg                                                                                                    <  7:410996/132‑1 (MQ=255)                                                           acaGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑gg                                                                                                    >  8:410996/1‑132 (MQ=255)                                                              gAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGc                                                                                <  2:30065/149‑1 (MQ=255)                                                                   cGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCACTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGAtc                                                                           >  4:132163/1‑149 (MQ=255)                                                                   cGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCACTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGAtc                                                                           >  4:126196/1‑149 (MQ=255)                                                                    gTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGAtc                                                                           <  8:267631/148‑1 (MQ=255)                                                                     tAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGAtcgt                                                                         >  5:296092/1‑149 (MQ=255)                                                                       aGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTATGGGGTACAGCCAGATCCGCACGCCGAtcgtc                                                                        <  1:19999/149‑1 (MQ=255)                                                                       aGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGAtcgtcg                                                                       <  5:280669/149‑1 (MQ=255)                                                                       aGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGAtcgtc                                                                        >  3:22147/1‑148 (MQ=255)                                                                       aGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGAtcgtc                                                                        >  3:16306/1‑148 (MQ=255)                                                                         gACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAg                                                                     >  4:236283/1‑149 (MQ=255)                                                                         gACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAg                                                                     >  3:57719/1‑149 (MQ=255)                                                                          aCTTTTCTTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCGGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAg                                                                     <  3:296773/148‑1 (MQ=255)                                                                           cTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCACTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGtt                                                                   <  3:132163/149‑1 (MQ=255)                                                                           cTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCACTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGtt                                                                   <  3:126196/149‑1 (MQ=255)                                                                               tCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCAc                                                                <  4:57719/148‑1 (MQ=255)                                                                                cGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCAcc                                                               <  6:337064/148‑1 (MQ=255)                                                                                 gTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCAcc                                                               >  2:406468/1‑147 (MQ=255)                                                                                 gTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCAcc                                                               >  1:335108/1‑147 (MQ=255)                                                                                  tGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAg                                                            <  3:18881/149‑1 (MQ=255)                                                                                  tGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAg                                                            >  3:296311/1‑149 (MQ=255)                                                                                    aGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGc                                                           <  4:296311/148‑1 (MQ=255)                                                                                    aGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCt                                                          <  1:406468/149‑1 (MQ=255)                                                                                      cAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGt                                                        >  5:411523/1‑149 (MQ=255)                                                                                        aaTCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTg                                                         >  4:257556/1‑146 (MQ=255)                                                                                        aaTCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTg                                                         <  3:257556/146‑1 (MQ=255)                                                                                        aaTCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTc                                                      <  6:83144/149‑1 (MQ=255)                                                                                          tCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCa                                                     >  2:66092/1‑148 (MQ=255)                                                                                            gCTGCAATCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAAc                                                  <  1:414993/149‑1 (MQ=255)                                                                                            gCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCACTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGAGCGTCGAGTTCACCGAGCTGTTCATAc                                                  >  6:215332/1‑149 (MQ=255)                                                                                              tGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGc                                                >  4:176961/1‑149 (MQ=255)                                                                                                cAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTc                                              >  8:408568/1‑149 (MQ=255)                                                                                                 aaGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCg                                             <  3:236283/149‑1 (MQ=255)                                                                                                      aTCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCgg                                        >  1:435757/1‑149 (MQ=255)                                                                                                       tCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCGGt                                       <  2:435757/149‑1 (MQ=255)                                                                                                                tGAACGACATCCTGCCAGAACAGTCGGCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCGGTGAAGTGAc                               >  1:225058/1‑149 (MQ=255)                                                                                                                tGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCGGTGAAGTGAcc                              >  1:61865/1‑149 (MQ=255)                                                                                                                 gAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCGGTGAAGTGAcc                              >  8:73717/1‑148 (MQ=255)                                                                                                                           cTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGGTGTTCAAACGCGCGAGCGGTGACGTGACCGCCAtcgtcg                    >  3:299957/1‑148 (MQ=255)                                                                                                                           cTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCGGTGAAGTGACCGACATCGTCGa                   <  3:206095/149‑1 (MQ=255)                                                                                                                              ccAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGa                                                                      <  2:21664/95‑1 (MQ=255)                                                                                                                              ccAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGa                                                                      >  1:21664/1‑95 (MQ=255)                                                                                                                               cAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCGGTGAAGTGACCGACATCGTCGaaaaa               <  6:197071/149‑1 (MQ=255)                                                                                                                                aGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCGGTGAAGTGACCGACATCGTCGaaa                 >  7:178570/1‑146 (MQ=255)                                                                                                                                aGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCGGTGAAGTGACCGACATCGTCGAAAAag              <  3:176961/149‑1 (MQ=255)                                                                                                                                  aaCAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCGGTGAAGTGACCGACATCGTCGAAAAagag            <  5:310959/149‑1 (MQ=255)                                                                                                                                   aCAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCGGTGAAGTGACCGACATCGTCGAAAAagag            <  2:156160/148‑1 (MQ=255)                                                                                                                                    cAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCGGTGAAGTGACCGACATCGTCGAAAAAGAGAt          <  6:296092/149‑1 (MQ=255)                                                                                                                                    cAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCGGTGAAGTGACCGACATCGTCGAAAAAGAGAt          <  2:335108/149‑1 (MQ=255)                                                                                                                                            cgcTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCGGTGAAGTGACCGACATCGTCGAAAAAGAGATGTACAcc   <  5:54656/149‑1 (MQ=255)                                                                                                                                             gcTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCGGTGAAGTGACCGACATCGTCGAAAAAGAGATGTACACCt  >  4:270415/1‑149 (MQ=255)                                                                                                                                             gcTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCGGTGAAGTGACCGACATCGTCGAAAAAGAGATGTACACCt  >  4:418804/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                              CCAGAAAGCGGCCGAAAAGGCCGAGGCCGACGCGGCGTTGATCGCTCGTGGCTGACACACAGAATTCGTAAGGACTTTTCGTGAGCAAATCGCTGCAAGCCATCCGTGGCATGAACGACATCCTGCCAGAACAGTC‑GCCGCTGTGGCGCTACTTCGAAGGCACAGTGGCCGGGCTGCTGGACACCTAT‑GGGTACAGCCAGATCCGCACGCCGATCGTCGAGTTCACCGAGCTGTTCAAACGCTCGATCGGTGAAGTGACCGACATCGTCGAAAAAGAGATGTACACCT  >  NC_002947/988069‑988356

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