Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	607,378	T→A	5.9%	R94R (CGT→CGA)	ribA →	GTP cyclohydrolase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	607,378	0	T	A	5.9%	236.3 / 2.6	79	R94R (CGT→CGA)	ribA	GTP cyclohydrolase II
Reads supporting (aligned to +/- strand):  ref base T (35/39);  new base A (2/3);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

CGGTACTGGGGCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAG  >  NC_002947/607230‑607525                                                                                                                                                     |                                                                                                                                                     cGGTACTGGGGCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGt                                                                                                                                                      >  7:142168/1‑149 (MQ=255)  ggTACTGGGGCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTg                                                                                                                                                     >  4:124290/1‑149 (MQ=255)        ggggCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCg                                                                                                                                               >  3:137583/1‑149 (MQ=255)         gggCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATc                                                                                                                                                <  3:220737/147‑1 (MQ=255)         gggCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCg                                                                                                                                               <  7:208264/148‑1 (MQ=255)          ggCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGc                                                                                                                                             <  6:318399/149‑1 (MQ=255)              ccTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGATGCCACGGCcttc                                                                                                                                         <  1:66975/149‑3 (MQ=255)              ccTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGGGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCctgc                                                                                                                                         <  6:285508/149‑1 (MQ=255)                tGCATTCCGAGTGCCTGACCGGCGATGCACTGTTAAGCCAGCGCTGCGACTGCGGCTCGCACTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGCGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCctgct                                                                                                                                        >  1:29780/1‑148 (MQ=255)                  cATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGaa                                                                                                                                     <  4:137583/149‑1 (MQ=255)                   atcccGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGaa                                                                                                                                     >  3:181060/4‑148 (MQ=255)                             ccTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑Ccgc                                                                                                                          <  7:237230/148‑1 (MQ=255)                             ccTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑Ccgc                                                                                                                          <  2:25470/148‑1 (MQ=255)                               tGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCctgct                                                                                                                                        <  8:20061/133‑1 (MQ=255)                               tGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCctgct                                                                                                                                        >  7:20061/1‑133 (MQ=255)                                 aCCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGATGCCACGGCcttcc                                                                                                                                        >  3:161186/1‑128 (MQ=255)                                 aCCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGATGCCACGGCcttcc                                                                                                                                        <  4:161186/131‑4 (MQ=255)                                   cGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTAc                                                                                                                    >  3:348278/1‑148 (MQ=255)                                     gCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAg                                                                                                                 >  7:331334/1‑149 (MQ=255)                                     gCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGATGCCACGGCcttcctg                                                                                                                                      <  6:387338/129‑6 (MQ=255)                                     gCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGATGCCACGGCcttcctg                                                                                                                                      >  5:387338/1‑124 (MQ=255)                                      cGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGTCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGc                                                                                                                >  4:168895/1‑149 (MQ=255)                                       gATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGTCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCt                                                                                                               <  3:168895/149‑1 (MQ=255)                                            aCTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGc                                                                                                                <  5:324021/143‑1 (MQ=255)                                            aCTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGc                                                                                                                >  6:324021/1‑143 (MQ=255)                                                  cagccagcGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATggc                                                                                                     <  4:226328/148‑1 (MQ=255)                                                       agcGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGcc                                                                                               <  2:124559/149‑1 (MQ=255)                                                        gcgcTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCTTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCg                                                                                              <  1:268174/149‑1 (MQ=255)                                                          gcTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGac                                                                                            >  8:262402/1‑149 (MQ=255)                                                           cTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGaca                                                                                           <  3:322119/149‑1 (MQ=255)                                                              cGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCg                                                                                        >  6:174232/1‑149 (MQ=255)                                                                 ctccGGCTAGCAGTTGGAAGCCGCGCTGCAGGCCATCGCACGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACGACAT‑CCGAGCCTACGAGCTGCAGGATGGCGGCGCCGACACAGtag                                                                                     >  6:282926/4‑147 (MQ=255)                                                                 cTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCa                                                                                                            >  6:91668/1‑126 (MQ=255)                                                                 cTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCa                                                                                                            <  5:91668/126‑1 (MQ=255)                                                                 cTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTc                                                                                      >  2:351612/1‑148 (MQ=255)                                                                     ggCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGc                                                                                 <  8:331334/149‑1 (MQ=255)                                                                          gCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAAc                                                                             >  4:99184/1‑148 (MQ=255)                                                                           cAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGa                                                                           <  3:124290/149‑1 (MQ=255)                                                                              ttGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAgcg                                                                        >  8:112439/1‑149 (MQ=255)                                                                                  aaGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTTTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTg                                                                    <  4:360290/149‑1 (MQ=255)                                                                                      cgcgCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCt                                                                >  3:17048/1‑149 (MQ=255)                                                                                        cgcTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTCGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTc                                                              >  5:25491/1‑149 (MQ=255)                                                                                             cAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGc                                                         <  4:17048/149‑1 (MQ=255)                                                                                             cAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGc                                                         >  3:158945/1‑149 (MQ=255)                                                                                              aGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCg                                                                                              >  3:122794/1‑111 (MQ=255)                                                                                              aGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCg                                                                                              <  4:122794/111‑1 (MQ=255)                                                                                               ggCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCg                                                       <  7:112439/149‑1 (MQ=255)                                                                                                gCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCCTCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCg                                                       <  8:142168/148‑1 (MQ=255)                                                                                                 ccATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGAc                                                     >  7:54744/1‑149 (MQ=255)                                                                                                 ccATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGAc                                                     <  2:29780/149‑1 (MQ=255)                                                                                                   aTCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCCGCGCCTACGAGCTGCAGGATggcg                                                                                                    <  1:4566/101‑1 (MQ=255)                                                                                                   aTCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATggcg                                                                                                    >  2:4566/1‑100 (MQ=255)                                                                                                   aTCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCa                                                   >  4:348749/1‑149 (MQ=255)                                                                                                      gCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTg                                                          >  8:218206/1‑139 (MQ=255)                                                                                                      gCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTg                                                          <  7:218206/139‑1 (MQ=255)                                                                                                       cccGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGt                                               >  4:49098/1‑149 (MQ=255)                                                                                                       cccGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGt                                               >  4:334842/1‑149 (MQ=255)                                                                                                        ccGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGc                                                 <  4:198822/146‑1 (MQ=255)                                                                                                        ccGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGc                                                 >  3:198822/1‑146 (MQ=255)                                                                                                            gAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTa                                          >  8:245468/1‑149 (MQ=255)                                                                                                                gCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGcc                                      <  7:245468/149‑1 (MQ=255)                                                                                                                  cGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCat                                    <  4:337677/149‑1 (MQ=255)                                                                                                                      gggTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATAtgc                                <  7:262402/149‑1 (MQ=255)                                                                                                                      gggTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATAtgc                                <  8:54744/149‑1 (MQ=255)                                                                                                                        gtgtTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATAtgcct                              >  2:257021/1‑149 (MQ=255)                                                                                                                         tgttgttgtACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATAtgcctg                             <  3:348749/149‑1 (MQ=255)                                                                                                                              ttgtACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGAt                        <  4:158945/149‑1 (MQ=255)                                                                                                                                gtACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCAGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCg                                                             <  6:179140/110‑1 (MQ=255)                                                                                                                                gtACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCAGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCg                                                             >  5:179140/1‑110 (MQ=255)                                                                                                                                   ccTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCat                                    <  7:11394/132‑1 (MQ=255)                                                                                                                                   ccTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCat                                    >  8:11394/1‑132 (MQ=255)                                                                                                                                    cTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCACCGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTc                    >  8:354493/1‑147 (MQ=255)                                                                                                                                    cTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTc                    >  8:272161/1‑147 (MQ=255)                                                                                                                                    cTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGa                  <  1:351612/149‑1 (MQ=255)                                                                                                                                       cGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAg                 <  5:174232/147‑1 (MQ=255)                                                                                                                                         tCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCg                                                <  6:300688/114‑1 (MQ=255)                                                                                                                                         tCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCg                                                >  5:300688/1‑114 (MQ=255)                                                                                                                                                 gCCGTGGCATCGGCCTGCTGAACTAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCTAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTc     <  5:282926/149‑1 (MQ=255)                                                                                                                                                 gCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTc     <  3:49098/149‑1 (MQ=255)                                                                                                                                                    gTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAg  >  7:266005/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CGGTACTGGGGCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAG  >  NC_002947/607230‑607525

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