Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,857,961	T→A	5.1%	L170L (CTT→CTA) ‡	PP_5134 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,857,961	0	T	A	5.1%	248.4 / 2.0	78	L170L (CTT→CTA) ‡	PP_5134	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/38);  new base A (2/2);  total (38/40)
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.94e-01

GCCTGTATGTGATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGCAC  >  NC_002947/5857813‑5858106                                                                                                                                                     |                                                                                                                                                  gCCTGTATGTGATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCtt                                                                                                                                                   >  5:3036/1‑149 (MQ=255)          tgATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCtt                                                                                                                                          <  2:110307/149‑1 (MQ=255)           gATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCCGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTg                                                                                                                                         >  8:218562/1‑149 (MQ=255)            atagATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTTGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGAACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTg                                                                                                                                         <  5:231811/145‑1 (MQ=255)            aTCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTg                                                                                                                                         <  6:114159/148‑1 (MQ=255)            aTCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTg                                                                                                                                         <  5:230704/148‑1 (MQ=255)             tCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTggg                                                                                                                                       >  1:25351/1‑149 (MQ=255)              cGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCATGGGc                                                                                                                                      >  2:111545/1‑149 (MQ=255)                aTGCGCGCTAGGCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCACCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCtt                                                                                                                                    <  3:13169/149‑1 (MQ=255)                  gcgcgcGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCATCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTc                                                                                                                                   <  5:25218/148‑1 (MQ=255)                   cgcgcgAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCAc                                                                                                                                 <  4:68033/149‑1 (MQ=255)                   cgcgcgAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCAc                                                                                                                                 <  6:3036/149‑1 (MQ=255)                      gcgAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCa                                                                                                                              >  2:77026/1‑149 (MQ=255)                      gcgAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCa                                                                                                                              >  5:115807/1‑149 (MQ=255)                      gcgAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCAAGCCGGGCAAGGGCTTCa                                                                                                                                  <  8:48110/145‑1 (MQ=255)                      gcgAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCAAGCCGGGCAAGGGCTTCa                                                                                                                                  >  7:48110/1‑145 (MQ=255)                         agagggTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCAt                                                                                                                             <  1:337811/143‑1 (MQ=255)                         aGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTg                                                                                                                           <  5:108030/149‑1 (MQ=255)                            cGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTa                                                                                                                        >  3:68401/1‑149 (MQ=255)                                    gCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGa                                                                                                                <  8:35470/149‑1 (MQ=255)                                      cAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGa                                                                                                                >  5:72586/1‑147 (MQ=255)                                      cAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATCGGTACGGTTCGACa                                                                                                              <  7:334182/149‑1 (MQ=255)                                      cAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCAAGCCGGGcaa                                                                                                                                          <  3:174008/121‑3 (MQ=255)                                      cAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCAAGCCGGGcaa                                                                                                                                          >  4:174008/1‑119 (MQ=255)                                       aaGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAg                                                                                                             >  7:206080/1‑149 (MQ=255)                                          tACCGCGACTAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCGCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGt                                                                                                          <  2:218499/149‑1 (MQ=255)                                            ccGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGACGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTgg                                                                                                        <  5:102416/149‑1 (MQ=255)                                                   gaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGc                                                                                                  >  8:180765/1‑148 (MQ=255)                                                   gaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGc                                                                                                  <  7:180765/148‑1 (MQ=255)                                                     aCTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGc                                                                                                  >  5:36795/1‑146 (MQ=255)                                                         gAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGa                                                                                           >  1:263096/1‑149 (MQ=255)                                                          aaGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGAt                                                                                          <  8:119459/149‑1 (MQ=255)                                                           aGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATg                                                                                         <  1:77026/149‑1 (MQ=255)                                                             tGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGt                                                                                       <  6:115807/149‑1 (MQ=255)                                                             tGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGt                                                                                       <  1:195362/149‑1 (MQ=255)                                                              gCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGt                                                                                       >  8:18182/1‑148 (MQ=255)                                                                 ggCAAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCCAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTAGACAGCGTGGCCCCGCCAGAGGATGGTGAgc                                                                                   >  4:72499/1‑149 (MQ=255)                                                                  gCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAgc                                                                                   <  1:228930/148‑1 (MQ=255)                                                                      aCTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGt                                                                              >  6:35376/1‑149 (MQ=255)                                                                       cTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGtt                                                                             >  5:169503/1‑149 (MQ=255)                                                                         gCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAgg                                                                                            <  4:390018/132‑1 (MQ=255)                                                                         gCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAgg                                                                                            >  3:390018/1‑132 (MQ=255)                                                                          cGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGttt                                                                            <  4:122057/147‑1 (MQ=255)                                                                          cGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGACCCCGCCAGAGGATGGTGAGCGCCGTATg                                                                           >  2:165143/1‑148 (MQ=255)                                                                           ggCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTa                                                                         <  2:25351/149‑1 (MQ=255)                                                                             cAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTAc                                                                        >  1:386547/1‑148 (MQ=255)                                                                                gCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGaaa                                                                    <  7:18182/149‑1 (MQ=255)                                                                                       ctgcCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGgc                                                             <  1:154147/149‑1 (MQ=255)                                                                                        tgcCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGgcg                                                            >  7:19264/1‑149 (MQ=255)                                                                                          ccGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCgct                                                          <  2:263096/149‑1 (MQ=255)                                                                                          ccGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCgct                                                          >  7:71323/1‑149 (MQ=255)                                                                                           cGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCgctg                                                         >  6:29261/1‑149 (MQ=255)                                                                                               gCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCgctgctgc                                                     >  3:38875/1‑149 (MQ=255)                                                                                                 tCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGt                                                                          >  2:178780/1‑126 (MQ=255)                                                                                                 tCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGt                                                                          <  1:178780/126‑1 (MQ=255)                                                                                                  cAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTgg                                                  <  4:68401/149‑1 (MQ=255)                                                                                                        cgtcgCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACTGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAAc                                             <  6:72586/148‑1 (MQ=255)                                                                                                          tcgCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAAcgcc                                          >  8:150023/1‑149 (MQ=255)                                                                                                              cagcagccagcaCCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAAcgcccgcc                                      <  2:259441/149‑1 (MQ=255)                                                                                                              cagcagccagcaCCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAAcgcccgcc                                      >  3:370872/1‑149 (MQ=255)                                                                                                                  agccagcaCCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGgc                                  <  6:36795/149‑1 (MQ=255)                                                                                                                   gccagcaCCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGgcg                                 >  1:293557/1‑149 (MQ=255)                                                                                                                           cgcgAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACg                         >  3:377281/1‑149 (MQ=255)                                                                                                                             cgAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGc                       <  8:71323/149‑1 (MQ=255)                                                                                                                             cgAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGc                       <  4:377281/149‑1 (MQ=255)                                                                                                                             cgAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGc                       <  4:370872/149‑1 (MQ=255)                                                                                                                              gAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCt                      <  8:206080/149‑1 (MQ=255)                                                                                                                                 ccTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAGCGGCTGAt                   >  4:281186/1‑149 (MQ=255)                                                                                                                                 ccTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGAt                   >  4:6786/1‑149 (MQ=255)                                                                                                                                  cTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATc                  >  5:312915/1‑149 (MQ=255)                                                                                                                                   tGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATcg                 >  5:181327/1‑149 (MQ=255)                                                                                                                                     aaTGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATcgac               <  5:35376/149‑1 (MQ=255)                                                                                                                                       tGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATcgacga             <  1:111545/149‑1 (MQ=255)                                                                                                                                           gTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACcag         >  4:176745/1‑149 (MQ=255)                                                                                                                                           gTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACcag         >  6:397888/1‑149 (MQ=255)                                                                                                                                                  ccTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGcac  <  7:150023/149‑1 (MQ=255)                                                                                                                                                  ccTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGATGCAACGGCTGATCGACGACCAGCAGGca   >  3:139205/1‑148 (MQ=255)                                                                                                                                                   cTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGca   >  3:96292/1‑147 (MQ=255)                                                                                                                                                   cTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGca   >  1:69581/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GCCTGTATGTGATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGCAC  >  NC_002947/5857813‑5858106

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