Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,576,052	T→A	5.5%	V150D (GTC→GAC)	ybhP‑II →	phosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,576,052	0	T	A	5.5%	223.7 / 2.3	73	V150D (GTC→GAC)	ybhP‑II	phosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (33/36);  new base A (2/2);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

CCTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCTTGG  >  NC_002947/4575907‑4576199                                                                                                                                                   |                                                                                                                                                    ccTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTg                                                                                                                                                  >  2:167538/1‑149 (MQ=255)   tGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCg                                                                                                                                                >  8:176478/1‑149 (MQ=255)     aCGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTg                                                                                                                                              >  4:318561/1‑149 (MQ=255)        gCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCAc                                                                                                                                           <  8:61571/149‑1 (MQ=255)        gCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCACGCAGACCGCATGCAc                                                                                                                                           <  2:334360/149‑1 (MQ=255)             cacGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGGGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTc                                                                                                                                      <  1:205342/149‑1 (MQ=255)               cGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGACAGGTGCACGCAGATCGCATGCACCTggcct                                                                                                                                    <  4:332358/149‑5 (MQ=255)                ggCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGtt                                                                                                                                   >  5:318660/1‑149 (MQ=255)                       cGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGa                                                                                                                             <  5:126175/148‑1 (MQ=255)                        gTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGaaa                                                                                                                           <  4:311090/149‑1 (MQ=255)                           gCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGc                                                                                                                         >  2:281112/1‑148 (MQ=255)                              gTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCAcc                                                                                                                     >  3:336895/1‑149 (MQ=255)                               tCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCAc                                                                                                                      <  1:167538/147‑1 (MQ=255)                               tCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGA‑ACTGTGTGCTGGAGGTGCCGGTGCAAGGCCAGGTGCATGCTGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCAc                                                                                                                      <  6:182345/147‑1 (MQ=255)                                 cAAATACACGATCATCGAGCACCGCAACATGGACCTGTCGATCACCGGCCCCGAACGCCCCGGCCTGTTCC‑ACTGTGTGCTGGAGGTGCAGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGATCCTGGAAAGCAACCgg                                                                                                                   >  4:151386/1‑148 (MQ=255)                                  aaaTACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGc                                                                                                                  >  7:103125/1‑148 (MQ=255)                                    aTACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCaa                                                                                                                >  1:4469/1‑148 (MQ=255)                                     tACGCGATCATCGAGCACCGCAACCTTGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGCTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCaaaa                                                                                                              <  4:30703/149‑1 (MQ=255)                                       cccGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCTGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAag                                                                                                             <  4:336895/148‑1 (MQ=255)                                        ccGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAag                                                                                                             >  7:114922/1‑147 (MQ=255)                                               tcGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCAGCAAAAGCAGCTGCAg                                                                                                    <  6:318660/149‑1 (MQ=255)                                                cGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCTACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAGAAGCAGCTGCAg                                                                                                    >  1:102809/1‑149 (MQ=255)                                                 gAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGt                                                                                                   <  8:114922/148‑1 (MQ=255)                                                       cGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAagcagc                                                                                                         >  3:211326/1‑136 (MQ=255)                                                       cGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAagcagc                                                                                                         <  4:211326/136‑1 (MQ=255)                                                           aCCTGGACGTGTCGATCACCGGCCCCGATCGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTaaa                                                                                        <  5:201593/149‑1 (MQ=255)                                                           aCCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTaaa                                                                                        <  6:126368/149‑1 (MQ=255)                                                            ccTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAAc                                                                                       >  5:367342/1‑149 (MQ=255)                                                                gACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAActgc                                                                                    <  7:176478/148‑1 (MQ=255)                                                                        gATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGccc                                                                           >  4:268829/1‑149 (MQ=255)                                                                        gATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCCGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGccc                                                                           >  6:358499/1‑149 (MQ=255)                                                                         aTCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCt                                                                          >  3:13572/1‑149 (MQ=255)                                                                           cACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGc                                                                        >  6:241617/1‑149 (MQ=255)                                                                              cGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAActgc                                                                                    <  6:370972/134‑1 (MQ=255)                                                                              cGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAActgc                                                                                    >  5:370972/1‑134 (MQ=255)                                                                              cGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCgg                                                                     <  2:4469/149‑1 (MQ=255)                                                                                   ccGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAAc                                                                >  4:80582/1‑149 (MQ=255)                                                                                    cGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACg                                                               >  6:264504/1‑149 (MQ=255)                                                                                                    ttGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGc                                               >  2:202351/1‑149 (MQ=255)                                                                                                         aCTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAGAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGc                                           <  2:102809/149‑1 (MQ=255)                                                                                                         aCTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGc                                           <  4:267096/149‑1 (MQ=255)                                                                                                           tgtgtgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGa                                         >  3:30238/1‑149 (MQ=255)                                                                                                           tgtgtgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGCGATGATCACCggggc                                         >  1:25256/1‑146 (MQ=255)                                                                                                             tgtgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGa                                         <  1:281112/147‑1 (MQ=255)                                                                                                                gCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTc                                     >  8:292243/1‑148 (MQ=255)                                                                                                                   ggAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAAc                                  >  4:368944/1‑148 (MQ=255)                                                                                                                      ggTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACt                              <  5:264504/149‑1 (MQ=255)                                                                                                                      ggTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACt                              <  1:202351/149‑1 (MQ=255)                                                                                                                        tGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCGCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTgg                            <  5:53842/149‑1 (MQ=255)                                                                                                                        tGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTgg                            <  7:292243/149‑1 (MQ=255)                                                                                                                         tttccagCAACGACAGGTGCACGCAGACCGCATGCACCTGGCCTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGGGACTTCAACGACTgg                            >  3:88995/7‑148 (MQ=25)                                                                                                                          ccGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGat                                                     <  4:206461/122‑1 (MQ=255)                                                                                                                            ggggCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCCCCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGaa                          >  2:70770/1‑147 (MQ=255)                                                                                                                            ggggCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGaa                          >  5:375118/1‑147 (MQ=255)                                                                                                                             gggCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAg                       <  3:80582/149‑1 (MQ=255)                                                                                                                             gggCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAg                       >  4:129118/1‑149 (MQ=255)                                                                                                                             gggCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAg                       >  2:332849/1‑149 (MQ=255)                                                                                                                                 aaGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCAt                   <  3:368944/149‑1 (MQ=255)                                                                                                                                 aaGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCAt                   <  8:103125/149‑1 (MQ=255)                                                                                                                                 aaGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCAt                   <  6:367342/149‑1 (MQ=255)                                                                                                                                 aaGGCCAGGTGCATGCCGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGTCAAAAGCATCTGCAGTTGCTGCGTAAACTTCTGGACGCCCTGCCCGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCAt                   <  4:296019/149‑1 (MQ=255)                                                                                                                                  aGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGCCTTCAACGACTGGAAACGCCATg                  >  6:332379/1‑149 (MQ=255)                                                                                                                                  aGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATg                  <  4:30238/149‑1 (MQ=255)                                                                                                                                   ggCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATgg                 <  4:13572/149‑1 (MQ=255)                                                                                                                                   ggCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATgg                 <  1:332849/149‑1 (MQ=255)                                                                                                                                    gacAGGTGCACGCAGACCGCATGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATgg                 >  1:174429/3‑148 (MQ=255)                                                                                                                                    gCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGa                                                      >  1:331434/1‑111 (MQ=255)                                                                                                                                    gCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGa                                                      <  2:331434/111‑1 (MQ=255)                                                                                                                                       aGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTTCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAAc             <  5:241617/149‑1 (MQ=255)                                                                                                                                              tGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCAccc      >  6:205299/1‑149 (MQ=255)                                                                                                                                               gaggTCTGAGTGCACCTGTCGTTGCTGGAAAGACACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTCCTGGCCGCCCTCCCCGCCAACCCCACGGAGATCATCGCCGGCGCCATCAACGACTGGAAAAGCCATGGCAACCGCTccca     >  6:188599/3‑148 (MQ=255)                                                                                                                                                  gTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCTTgg  <  3:318561/149‑1 (MQ=255)                                                                                                                                                  gTCTGCGTGCACCCGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCTTgg  <  5:358499/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CCTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGC‑ACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCTTGG  >  NC_002947/4575907‑4576199

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