Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	607,378	T→A	6.0%	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	6.0%	203.7 / 2.4	67	R94R (CGT→CGA)	ribA	GTP cyclohydrolase II
Reads supporting (aligned to +/- strand):  ref base T (28/35);  new base A (2/2);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

CGGTACTGGGGCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGT  >  NC_002947/607230‑607521                                                                                                                                                     |                                                                                                                                                cGGTACTGGGGCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGt                                                                                                                                                 <  2:372766/149‑1 (MQ=255)  ggTACTGGGGCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTg                                                                                                                                                <  1:161825/149‑1 (MQ=255)  ggTACTGGGGCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTg                                                                                                                                                <  6:197397/149‑1 (MQ=255)       tGGGGCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATc                                                                                                                                           >  4:104180/1‑149 (MQ=255)         gggCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCgg                                                                                                                                         >  7:281468/1‑149 (MQ=255)             gcCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCct                                                                                                                                      >  8:265570/1‑148 (MQ=255)                 gCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCctgctg                                                                                                                                  >  6:119535/1‑148 (MQ=255)                 gCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGATGCCACGGCcttcc                                                                                                                                   <  3:60333/147‑4 (MQ=255)                 gCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGATGCCACGGCcttcc                                                                                                                                   >  4:60333/1‑144 (MQ=255)                  cATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGATGCCACGGCcttcct                                                                                                                                  <  7:320786/147‑5 (MQ=255)                      ccGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACaa                                                                                                                             <  8:281468/148‑1 (MQ=255)                       cGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAAc                                                                                                                               >  5:97733/1‑145 (MQ=255)                        gAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAt                                                                                                                          <  6:206274/149‑1 (MQ=255)                          gTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATcc                                                                                                                        >  8:343196/1‑149 (MQ=255)                              cTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCgcgc                                                                                                                    <  8:257289/149‑1 (MQ=255)                               tGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTCGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGATGTTGTACCTGCGTCAGGAAGGCCGAGGCATCCGCCTGCTGAACAAGATCCgcgc                                                                                                                    >  5:341586/1‑148 (MQ=255)                                         tGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGc                                                                                                         <  7:343196/149‑1 (MQ=255)                                           ctcTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGGCGCTCTGCAGGCCATAGCCCTCGTAGGCCGTGGTGTTTTGTTGTACCTGCGTCAGGAAGGCGGTGGCATCGTCCTGCTGAACATGATCCGCTCCTACGAGCTGCa                                                                                                        <  8:9537/146‑1 (MQ=255)                                           cACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCa                                                                                                        <  6:272050/148‑1 (MQ=255)                                            aCTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAg                                                                                                       >  2:343472/1‑148 (MQ=255)                                                ttCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATgg                                                                                                  >  6:150542/1‑149 (MQ=255)                                                ttCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATgg                                                                                                  >  4:196627/1‑149 (MQ=255)                                                  cagccagcGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATggcg                                                                                                >  7:183161/1‑149 (MQ=255)                                                   agccagcGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATggcgg                                                                                               <  7:289236/149‑1 (MQ=255)                                                     ccagcGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGaa                                                                                                                                <  3:189863/114‑1 (MQ=255)                                                     ccagcGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGaa                                                                                                                                >  4:189863/1‑114 (MQ=255)                                                         cgcTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGa                                                                                         <  4:122182/149‑1 (MQ=255)                                                         cgcTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGa                                                                                         >  5:302470/1‑149 (MQ=255)                                                           cTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGaca                                                                                       <  5:97098/149‑1 (MQ=255)                                                                    cGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGaa                                                                               <  6:346961/148‑1 (MQ=255)                                                                    cGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGaa                                                                               <  4:226327/148‑1 (MQ=255)                                                                    cGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCATGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGaa                                                                               <  7:54558/148‑1 (MQ=255)                                                                     ggCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGc                                                                             <  5:232536/149‑1 (MQ=255)                                                                          gCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAAc                                                                         <  3:14739/148‑1 (MQ=255)                                                                             gTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAgc                                                                     <  5:150542/149‑1 (MQ=255)                                                                                ggAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAgcgc                                                                   <  3:104180/148‑1 (MQ=255)                                                                                  aaGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTg                                                                <  1:252749/149‑1 (MQ=255)                                                                                     ccGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGACGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGc                                                             >  4:366794/1‑149 (MQ=255)                                                                                       gcgcTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCtt                                                           >  2:227990/1‑149 (MQ=255)                                                                                         gcTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGACGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCg                                                         <  3:366794/149‑1 (MQ=255)                                                                                              aGGCCATCGCCCGCGAAGGCCGTGGTGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGc                                                     <  5:234886/148‑1 (MQ=255)                                                                                                gCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGa                                                  >  8:91970/1‑149 (MQ=255)                                                                                                gCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGa                                                  >  4:319596/1‑149 (MQ=255)                                                                                                  cATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCCACGAGCGCCTGGGCTTCGCTGCCGAc                                                 >  5:248675/1‑148 (MQ=255)                                                                                                   aTCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCa                                               <  7:265570/149‑1 (MQ=255)                                                                                                    tCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAg                                              >  6:366991/1‑149 (MQ=255)                                                                                                      gCCCGCGAAGGCCGTGGGGTGTTGTTGTTCCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCg                                            >  1:210367/1‑149 (MQ=255)                                                                                                         cgcgAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGa                                         <  3:319596/149‑1 (MQ=255)                                                                                                         cgcgAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGa                                         >  1:80018/1‑149 (MQ=255)                                                                                                         cgcgAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGa                                         <  6:97733/149‑1 (MQ=255)                                                                                                             aaGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTAc                                     <  5:119535/149‑1 (MQ=255)                                                                                                                    tagggTGTTGTTGTACCTGCGGCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCatat                              <  6:312127/147‑1 (MQ=255)                                                                                                                    tGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCatat                              <  3:196627/149‑1 (MQ=255)                                                                                                                         tgttgttgtACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCt                                                            <  2:222588/114‑1 (MQ=255)                                                                                                                          gttgttgtACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGATCAAGATCCGCGCGTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAAGGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGGCATAtgcctgc                        <  6:341586/149‑1 (MQ=255)                                                                                                                             gttgtACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGa                     >  1:326962/1‑149 (MQ=255)                                                                                                                              ttgtACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGAt                    >  3:110866/1‑149 (MQ=255)                                                                                                                                   ccTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGc                                                             >  3:86170/1‑103 (MQ=255)                                                                                                                                   ccTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGc                                                             <  4:86170/103‑1 (MQ=255)                                                                                                                                   ccTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTc                <  6:302470/148‑1 (MQ=255)                                                                                                                                   ccTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTc                <  7:91970/148‑1 (MQ=255)                                                                                                                                    cTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGAACCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTc                >  2:351595/1‑147 (MQ=255)                                                                                                                                      gCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGCCCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGc            >  2:59739/1‑149 (MQ=255)                                                                                                                                      gCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGc            >  2:20021/1‑149 (MQ=255)                                                                                                                                            ggAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTgg      <  1:343472/149‑1 (MQ=255)                                                                                                                                             gAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTggg     >  6:208517/1‑149 (MQ=255)                                                                                                                                               aGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCg   <  5:72847/149‑1 (MQ=255)                                                                                                                                                 gCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGt  >  4:329814/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                CGGTACTGGGGCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGGAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGT  >  NC_002947/607230‑607521

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