Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	772,750	T→G	5.4%	L138R (CTC→CGC)	PP_0664 →	homoserine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	772,750	0	T	G	5.4%	286.5 / 4.1	88	L138R (CTC→CGC)	PP_0664	homoserine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (39/44);  new base G (2/3);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCGCCGACATCATCGCCGAGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACCTCTAACTAT  >  NC_002947/772607‑772898                                                                                                                                                |                                                                                                                                                      aGGCGCCGACATCATCGCCGAGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAgg                                                                                                                                                  >  5:105555/1‑149 (MQ=255)   gcgcCGACATCATCGCCGAGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGcc                                                                                                                                                <  5:11286/149‑1 (MQ=255)    cgcCGACATCATCGCCGAGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGGGGGCAAGCACGTTGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCt                                                                                                                                               <  8:30894/149‑1 (MQ=255)      ccGACTTCATCTCAGAGGCCACATTCACCATCCCCAACGACGGTGAGCCGGCGGCAACCTTCTGCCTATGGGCGCTGGAGGCGTGCAAGCACGTGGTAACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTg                                                                                                                                             <  6:33358/149‑1 (MQ=255)      ccGACATCATCGCCGAGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTg                                                                                                                                             <  2:357112/149‑1 (MQ=255)            tcatcGCCGAGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAg                                                                                                                                       <  1:229216/149‑1 (MQ=255)             catcGCCGAGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAgg                                                                                                                                      >  2:106520/1‑149 (MQ=255)                    gAGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAg                                                                                                                                                   <  5:227149/129‑1 (MQ=255)                    gAGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAg                                                                                                                                                   >  6:227149/1‑129 (MQ=255)                     aGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAAt                                                                                                                              >  4:201815/1‑149 (MQ=255)                     aGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAAt                                                                                                                              <  1:250166/149‑1 (MQ=255)                     aGGCCACAATCACCACCCCCCCGGACGGGGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGCGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCGTGGCTAAGGCCAACAAt                                                                                                                              <  8:167737/149‑1 (MQ=255)                      ggCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAAtg                                                                                                                             >  6:312501/1‑149 (MQ=255)                        ccACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAAtgtg                                                                                                                           <  2:102499/149‑1 (MQ=255)                           caTTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTaa                                                                                                                                        <  7:36878/133‑1 (MQ=255)                           caTTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTaa                                                                                                                                        >  8:36878/1‑133 (MQ=255)                           caTTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCg                                                                                                                        <  2:141581/149‑1 (MQ=255)                            aTTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGt                                                                                                                       >  3:16148/1‑149 (MQ=255)                              tCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑Tc                                                                                                                    >  8:253918/1‑149 (MQ=255)                               cACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑Tc                                                                                                                    <  3:132921/148‑1 (MQ=255)                                aCCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑Tc                                                                                                                    >  7:297056/1‑147 (MQ=255)                                aCCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑Tc                                                                                                                    <  8:297056/147‑1 (MQ=255)                                 ccAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGaa                                                                                                                 >  5:320016/1‑149 (MQ=255)                                 ccAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGaa                                                                                                                 <  7:167958/149‑1 (MQ=255)                                 ccAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGAGTGCCTGGGCATGCAAGGcaat                                                                                                                                              <  8:36051/121‑4 (MQ=255)                                 ccAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGAGTGCCTGGGCATGCAAGGcaat                                                                                                                                              >  7:36051/1‑118 (MQ=255)                                  cAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCaa                                                                                                                                                    <  2:374111/114‑1 (MQ=255)                                  cAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCaa                                                                                                                                                    >  1:374111/1‑114 (MQ=255)                                    aCCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGt‑t                                                                                                                     >  3:277520/1‑142 (MQ=255)                                    aCCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGt‑t                                                                                                                     <  4:277520/142‑1 (MQ=255)                                       ccACGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGaa                                                                                                           <  6:299128/149‑1 (MQ=255)                                        gaaGGTCGGTGAGCCGGCGGCAACGTTGTGGCGCTGGGCGCTGGAGGCGGGCAAGCACTTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGGCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAg                                                                                                          <  5:215103/148‑1 (MQ=255)                                        cAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAg                                                                                                          <  3:145467/149‑1 (MQ=255)                                        cAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAg                                                                                                          <  1:106520/149‑1 (MQ=255)                                             aCGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGcc                                                                                                     <  4:268010/149‑1 (MQ=255)                                                   aGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCAt                                                                                                <  3:201815/148‑1 (MQ=255)                                                     ccGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAg                                                                                             >  2:371286/1‑149 (MQ=255)                                                     ccGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAg                                                                                             >  1:202426/1‑149 (MQ=255)                                                      cggcggCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGc                                                                                            <  4:109427/149‑1 (MQ=255)                                                           gcAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGca                                                                                        >  5:240198/1‑148 (MQ=255)                                                             aaCCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCacac                                                                                     >  1:288206/1‑149 (MQ=255)                                                             aaCCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGGCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCacac                                                                                     <  6:31511/149‑1 (MQ=255)                                                              aCCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACAcc                                                                                    >  6:387164/1‑149 (MQ=255)                                                                      ccGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGtgctgc                                                                            >  4:232668/1‑149 (MQ=255)                                                                      ccGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGtgctgc                                                                            >  5:227750/1‑149 (MQ=255)                                                                       cGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCg                                                                           >  7:294125/1‑149 (MQ=255)                                                                           gggCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATg                                                                       >  5:175697/1‑149 (MQ=255)                                                                            ggCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATgg                                                                      <  6:105555/149‑1 (MQ=255)                                                                              cgcTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGgcg                                                                    >  1:222093/1‑149 (MQ=255)                                                                                 tGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGgcgcgc                                                                 >  8:363088/1‑149 (MQ=255)                                                                                        gggCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGAccc                                                          <  6:240198/149‑1 (MQ=255)                                                                                           cAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTgg                                                       <  3:232668/149‑1 (MQ=255)                                                                                            aaGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTggc                                                      <  8:294125/149‑1 (MQ=255)                                                                                              gCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTggcgg                                                    <  2:89486/149‑1 (MQ=255)                                                                                              gCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGACCAGGAAATAAAGGCCTTGGATAATGACAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCAAACCGGTGCTGCGTATGGCGAGCCAGACACTggcgg                                                    >  4:231033/1‑149 (MQ=255)                                                                                                 cGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTggcgggcg                                                 >  8:336495/1‑149 (MQ=255)                                                                                                 cGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTggcgggcg                                                 <  6:320016/149‑1 (MQ=255)                                                                                                  gTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGgcgggcgc                                                <  7:253918/149‑1 (MQ=255)                                                                                                    ggTCACCACCAACAAGGGGACGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGATCCGCACACCGCTGATGCGTATGGCGCGCCAGACCCTGgcgggcgctc                                              >  6:168072/1‑147 (MQ=255)                                                                                                       caccacCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCt                                           >  3:356444/1‑149 (MQ=255)                                                                                                        accaccAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCt                                           <  4:16148/148‑1 (MQ=255)                                                                                                        accaccAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTACGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCtt                                          >  6:283827/1‑149 (MQ=255)                                                                                                         ccaccaACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTg                                         >  1:233376/1‑149 (MQ=255)                                                                                                          ctcctacaaGGGGCCGTTTGCCTTGCATGCGCAGGCTCTGAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCTCCGTCATTAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGAc                                                            <  8:147721/124‑1 (MQ=255)                                                                                                          caccaACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTgg                                        <  6:227750/149‑1 (MQ=255)                                                                                                          caccaACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTgg                                        >  7:362293/1‑149 (MQ=255)                                                                                                          caccaACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTAGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCCGTGCTGCGTATGGCGCGCCAGAc                                                            >  7:147721/1‑129 (MQ=255)                                                                                                             caacaaGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAg                                     <  4:356444/149‑1 (MQ=255)                                                                                                                    gggCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGa                              <  2:288206/149‑1 (MQ=255)                                                                                                                     ggCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAg                             <  5:312501/149‑1 (MQ=255)                                                                                                                      gCCGATTGCCTTTCATGCGCAGGCACGCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGAACCTGGCGGGCGCGGGCTTGGTAGGTTTCGAgg                            <  5:168072/149‑1 (MQ=255)                                                                                                                      gCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAgg                            <  2:222093/149‑1 (MQ=255)                                                                                                                            ggCCTTGAGTGCCTGGGCATGCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTc                                >  4:132386/2‑139 (MQ=255)                                                                                                                            ggCCTTGAGTGCCTGGGCATGCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTc                                <  3:132386/138‑1 (MQ=255)                                                                                                                              ccTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATa                                                                                                               <  8:73996/58‑1 (MQ=255)                                                                                                                              ccTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATa                                                                                                               >  7:73996/1‑58 (MQ=255)                                                                                                                               cTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATtctc                    >  4:87300/1‑148 (MQ=255)                                                                                                                                ttGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTGTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATtctc                    >  1:89346/1‑148 (MQ=255)                                                                                                                                ttGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCaa                  >  4:56730/1‑149 (MQ=255)                                                                                                                                ttGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCaa                  >  8:149071/1‑149 (MQ=255)                                                                                                                                   cATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACg                <  2:202426/148‑1 (MQ=255)                                                                                                                                      gCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTAc            >  6:25536/1‑149 (MQ=255)                                                                                                                                       cccAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTAc            <  1:371286/148‑1 (MQ=255)                                                                                                                                        ccAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACct          >  5:231708/1‑149 (MQ=255)                                                                                                                                         cAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTACGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACctc         <  5:283827/149‑1 (MQ=255)                                                                                                                                          aGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACctct        <  7:336495/149‑1 (MQ=255)                                                                                                                                            gCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACCTCTaa      <  6:231708/149‑1 (MQ=255)                                                                                                                                               cTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACCTCTAACTa   >  5:262954/1‑149 (MQ=255)                                                                                                                                                tCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACCTCTAACTAt  >  5:210454/1‑149 (MQ=255)                                                                                                                                                tCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACCTCTAACTAt  >  3:366510/1‑149 (MQ=255)                                                                                                                                                tCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACCTCTAACTAt  >  1:125545/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                      AGGCGCCGACATCATCGCCGAGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGT‑TCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACCCTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACCTCTAACTAT  >  NC_002947/772607‑772898

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