Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	607,379	G→T	6.4%	G95C (GGC→TGC)	ribA →	GTP cyclohydrolase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	607,379	0	G	T	6.4%	250.7 / 6.0	78	G95C (GGC→TGC)	ribA	GTP cyclohydrolase II
Reads supporting (aligned to +/- strand):  ref base G (33/39);  new base T (2/3);  total (36/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.47e-01

GGTACTGGGGCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAG  >  NC_002947/607231‑607525                                                                                                                                                     |                                                                                                                                                    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                                                                                                                                        >  7:20061/1‑133 (MQ=255)                              tGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCctgct                                                                                                                                        <  8:20061/133‑1 (MQ=255)                                aCCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGATGCCACGGCcttcc                                                                                                                                        <  4:161186/131‑4 (MQ=255)                                aCCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGATGCCACGGCcttcc                                                                                                                                        >  3:161186/1‑128 (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                                                                                                                >  6:324021/1‑143 (MQ=255)                                           aCTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGc                                                                                                                <  5:324021/143‑1 (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                                                                                                            <  5:91668/126‑1 (MQ=255)                                                                cTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCa                                                                                                            >  6:91668/1‑126 (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                                                         >  3:158945/1‑149 (MQ=255)                                                                                            cAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGc                                                         <  4:17048/149‑1 (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:334842/1‑149 (MQ=255)                                                                                                      cccGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGt                                               >  4:49098/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                                                >  5:300688/1‑114 (MQ=255)                                                                                                                                        tCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCg                                                <  6:300688/114‑1 (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)                                                                                                                                                     |                                                                                                                                                    GGTACTGGGGCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAT‑CCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAG  >  NC_002947/607231‑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: ‑