Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,942,315	G→C	6.3%	P67R (CCG→CGG)	PP_5209 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,942,315	0	G	C	6.3%	286.5 / 8.9	95	P67R (CCG→CGG)	PP_5209	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (43/45);  new base C (2/4);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

GACGTTGCTCATGCTGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATCAAGATCCACGCTTTCACGAA  >  NC_002947/5942170‑5942460                                                                                                                                                    |                                                                                                                                                  gACGTTGCTCATGCTGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCgg                                                                                                                                                  <  6:382532/147‑1 (MQ=255)     ttGCTCATGCTGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGcc                                                                                                                                            <  2:319426/149‑1 (MQ=255)       gCTCATGCTGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCCGCGGCTGcgc                                                                                                                                           <  8:304691/148‑3 (MQ=255)       gCTCATGCTGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCCGCGGCTGcgc                                                                                                                                           <  7:217354/148‑3 (MQ=255)        cTCATGCTGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGAGCCGCGGTTCGGGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGGCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGt                                                                                                                                          <  7:318369/148‑1 (MQ=255)         tCATGCTGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCACTTTTGAGGGGCCACGTACGCATTTAACACGCGCACCCCCGGGCTGCCGt                                                                                                                                          >  1:147322/1‑149 (MQ=255)         tCATGCTGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCCGCCGCGGGCTGCCGTc                                                                                                                                         >  6:67913/1‑148 (MQ=255)         tCATGCTGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTc                                                                                                                                         >  8:291484/1‑148 (MQ=255)          cATGCTGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGa                                                                                                                                       <  4:4138/149‑1 (MQ=255)           aTGCTGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCCGCGGCTGcgcg                                                                                                                                          >  3:297935/1‑142 (MQ=255)           aTGCTGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCCGCGGCTGcgcg                                                                                                                                          <  4:297935/145‑4 (MQ=255)            tGCTGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAgc                                                                                                                                     <  4:379631/149‑1 (MQ=255)            tGCTGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAgc                                                                                                                                     <  4:377593/149‑1 (MQ=255)            tGCTGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAgc                                                                                                                                     <  6:356735/149‑1 (MQ=255)             gCTGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTcgcgcg                                                                                                                                    >  4:346145/1‑149 (MQ=255)               tGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCCC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCa                                                                                                                                  >  4:416708/1‑149 (MQ=255)                gTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTAGTACACGCGCAGCCGCTGGCTGACGTCGAGCGCa                                                                                                                                  >  4:315006/1‑148 (MQ=255)                 tCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCa                                                                                                                                  >  6:372219/1‑147 (MQ=255)                 tCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACACGCAGCCGCGGGCTGCCGTCGAGCGCATa                                                                                                                                <  6:256937/149‑1 (MQ=255)                     ccGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTg                                                                                                                            <  2:68340/149‑1 (MQ=255)                       gTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGcc                                                                                                                          <  3:217522/149‑1 (MQ=255)                             tgGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACg                                                                                                                      >  1:414901/1‑147 (MQ=255)                               ggTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCagccatag                                                                                                                  >  1:372495/1‑142 (MQ=255)                               ggTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGaa                                                                                                                    <  3:50032/147‑1 (MQ=255)                               ggTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAgg                                                                                                                  <  5:154075/149‑1 (MQ=255)                               ggTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAgg                                                                                                                  <  7:291513/149‑1 (MQ=255)                                 tGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAgggc                                                                                                                >  7:225581/1‑149 (MQ=255)                                 tGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCCGCGGCTGcgcg                                                                                                                                          <  2:116576/123‑4 (MQ=255)                                 tGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCCGCGGCTGcgcg                                                                                                                                          >  1:116576/1‑120 (MQ=255)                                      aGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCt                                                                                                           >  5:426853/1‑149 (MQ=255)                                       gTGCCACCAGCTGGTTACGGATCAGCGTTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCATCGAAGGGCGGGCTc                                                                                                          <  7:224244/149‑1 (MQ=255)                                        tGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAAGGAAggggtgccc                                                                                                           >  5:224720/1‑141 (MQ=255)                                        tGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCt                                                                                                           >  3:62955/1‑147 (MQ=255)                                          ccaccaGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGa                                                                                                       >  2:273576/1‑149 (MQ=255)                                           caccaGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGAc                                                                                                      <  8:225581/149‑1 (MQ=255)                                            accaGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACTAAGGGCGGGCTCAGACt                                                                                                     <  2:414901/149‑1 (MQ=255)                                               aGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGAt                                                                                                  <  6:324475/149‑1 (MQ=255)                                                 cTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCATGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCCGGCTCAGACTGATg                                                                                                 >  3:332481/1‑148 (MQ=255)                                                  tGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGAGGGCTCAGACTGATg                                                                                                 >  5:256406/1‑147 (MQ=255)                                                    gTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTa                                                                                               <  5:333641/147‑1 (MQ=255)                                                         ggATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGAcc                                                                                        >  8:299639/1‑149 (MQ=255)                                                          gATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACt                                                                                                     <  5:63816/135‑1 (MQ=255)                                                          gATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACt                                                                                                     >  6:63816/1‑135 (MQ=255)                                                                 ggTTCGTTGGGCTTTACGGTCTTGGCCGCTTGGTCGCCCGTCAC‑‑TCGCAGGGCCACGTCAGCGTTGTAGACGCGCAGCCGCGGGCTGCCGTCTAGCGCATAGTTGCCAACGATGGGCGGGCTCAGACTGATGTATGCGACCTTGTGCAc                                                                                <  8:346099/149‑1 (MQ=255)                                                                  gTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCAcc                                                                               >  1:215622/1‑149 (MQ=255)                                                                     cGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCcct                                                                            >  4:428234/1‑149 (MQ=255)                                                                      gtggtgCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGACGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCcctt                                                                           >  5:309901/1‑149 (MQ=255)                                                                        ggtgCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCccttcc                                                                         >  4:202506/1‑149 (MQ=255)                                                                         gtgCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAGCGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCccttcct                                                                        >  6:100904/1‑149 (MQ=255)                                                                           gCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCAcccc                                                                             >  6:129658/1‑142 (MQ=255)                                                                           gCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCAcccc                                                                             <  5:129658/142‑1 (MQ=255)                                                                            cTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTg                                                                      <  1:273576/148‑1 (MQ=255)                                                                             tttACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGc                                                                    >  1:99368/1‑149 (MQ=255)                                                                              ttACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCGACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCt                                                                   >  7:377511/1‑149 (MQ=255)                                                                               tACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGc                                                                    <  5:372219/147‑1 (MQ=255)                                                                                aCCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCtt                                                                  <  6:135685/148‑1 (MQ=255)                                                                                aCCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTc                                                                 >  1:365676/1‑149 (MQ=255)                                                                                   gTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTc                                                              >  7:869/1‑149 (MQ=255)                                                                                        ggCCGCTTCGTCGCCTGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCa                                                         >  3:80634/1‑149 (MQ=255)                                                                                        ggCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCGCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGc                                                           <  4:156517/147‑1 (MQ=255)                                                                                           cGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATgg                                                      <  8:377511/149‑1 (MQ=255)                                                                                            gCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGc                                                     >  7:305652/1‑149 (MQ=255)                                                                                             cTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGcc                                                    <  7:299639/149‑1 (MQ=255)                                                                                                   cgCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGcc                                              >  8:293587/1‑149 (MQ=255)                                                                                                    gCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCCCCCCTTCCTTGGCTTCCTCGGCCCTGGccgccgccg                                             >  3:397166/1‑149 (MQ=255)                                                                                                     cccGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCgg                                            >  1:267627/1‑149 (MQ=255)                                                                                                     cccGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCgg                                            <  2:365676/149‑1 (MQ=255)                                                                                                     cccGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCATTGGCTTCCTCGGCCATGGCCGCAGCCgg                                            <  7:291484/149‑1 (MQ=255)                                                                                                       cGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCgg                                            >  1:379204/1‑147 (MQ=255)                                                                                                       cGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGgca                                          <  6:426853/149‑1 (MQ=255)                                                                                                        gTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGgcag                                         <  3:315006/149‑1 (MQ=255)                                                                                                        gTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGgcag                                         <  7:280521/149‑1 (MQ=255)                                                                                                                 gcAGGGCCACGTCATCCGTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCagcagcagc                                  <  6:68893/149‑1 (MQ=255)                                                                                                                 gcAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCagcagcagc                                  <  2:379204/149‑1 (MQ=255)                                                                                                                  cAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAgcg                                 <  3:202506/149‑1 (MQ=255)                                                                                                                   aGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAgcgc                                >  4:203479/1‑149 (MQ=255)                                                                                                                    gggCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGcc                               <  4:80634/149‑1 (MQ=255)                                                                                                                    gggCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGcc                               <  3:416708/149‑1 (MQ=255)                                                                                                                     ggCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAgcgc                                >  6:158490/1‑147 (MQ=255)                                                                                                                            tCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCAt                       >  4:341776/1‑149 (MQ=255)                                                                                                                              aGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATc                      <  3:341776/148‑1 (MQ=255)                                                                                                                              aGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATc                      <  3:203479/148‑1 (MQ=255)                                                                                                                               gCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATCaa                    <  2:147322/149‑1 (MQ=255)                                                                                                                                ccTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATc                      >  5:294817/1‑146 (MQ=255)                                                                                                                                ccTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATCAAg                   <  7:293587/149‑1 (MQ=255)                                                                                                                                ccTTGTACACGCGCAGCCGCGGGCTACCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATCAAg                   <  4:62955/149‑1 (MQ=255)                                                                                                                                        acGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCATCACCATCAAGATCCACGc           >  8:319649/1‑149 (MQ=255)                                                                                                                                        acGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATCAAGATCCACGc           >  1:426006/1‑149 (MQ=255)                                                                                                                                          gcgcAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATc                      <  5:147122/136‑1 (MQ=255)                                                                                                                                          gcgcAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATc                      >  6:147122/1‑136 (MQ=255)                                                                                                                                          gcgcAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATCAAGATCCACGCtt         >  3:110753/1‑149 (MQ=255)                                                                                                                                          gcgcAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATCAAGATCCACGCtt         >  3:88753/1‑149 (MQ=255)                                                                                                                                          gcgcAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATCAAGATCCACGCtt         >  2:95923/1‑149 (MQ=255)                                                                                                                                           cgcAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATCAAGATCCACGCttt        <  2:426006/149‑1 (MQ=255)                                                                                                                                            gcAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATCAAGATCCACGCTTTc       <  7:319649/149‑1 (MQ=255)                                                                                                                                              aGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATCAAGATCCACGCTTTCAc     <  7:3020/149‑1 (MQ=255)                                                                                                                                                 cgcgTGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATCAAGATCCACGCTTTCACGaa  <  8:305652/149‑1 (MQ=255)                                                                                                                                                 cgcgGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATCAAGATCCACGCTTTCACGaa  <  6:256406/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GACGTTGCTCATGCTGTCCACCGTCTGCTGGGTGAACAGTGCCACCAGCTGGTTACGGATCAGCGGTTCGTGGTGCTTTACCGTCTTGGCCGCTTCGTCGCCCGTCAC‑‑GCGCAGGGCCACGTCAGCCTTGTACACGCGCAGCCGCGGGCTGCCGTCGAGCGCATAGTTGCCAACGAAGGGCGGGCTCAGACTGATGTAGGCGACCTTGGGCACCCCTTCCTTGGCTTCCTCGGCCATGGCCGCAGCCGGCAGCAGCAGCGCCAGCACCATCAAGATCCACGCTTTCACGAA  >  NC_002947/5942170‑5942460

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