Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,835,980	G→C	5.6%	W156C (TGG→TGC)	PP_3387 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,835,980	0	G	C	5.6%	230.3 / 2.9	71	W156C (TGG→TGC)	PP_3387	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/35);  new base C (2/2);  total (34/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 = 8.38e-01

CACTTGAGCCACGCCAACATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGAGCTTTCAACTT  >  NC_002947/3835834‑3836121                                                                                                                                                   |                                                                                                                                              cacTTGAGCCACGCCAACATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGc                                                                                                                                             <  6:277393/149‑1 (MQ=255)    ttGAGCCACGCCAACATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGcccc                                                                                                                                          >  7:260751/1‑149 (MQ=255)        gCCACGCCAACATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGCGccacgggg                                                                                                                                      >  2:307183/1‑143 (MQ=255)         ccACGCCAACATCGTCACGGGCATGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTgg                                                                                                                                     <  3:19926/149‑1 (MQ=255)         ccACGCCAACATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTgg                                                                                                                                     >  7:226923/1‑149 (MQ=255)         ccACGCCAACATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTgg                                                                                                                                     >  4:199723/1‑149 (MQ=255)         ccACGCCAACATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGCGccacggggg                                                                                                                                     >  6:82805/1‑142 (MQ=255)           aCGCCAACATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGgcg                                                                                                                                   >  7:297293/1‑149 (MQ=255)            cGCCAACATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGgcgc                                                                                                                                  <  3:199723/149‑1 (MQ=255)             gCCAACATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGgcgga                                                                                                                                 >  5:377992/1‑147 (MQ=255)               cAACATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCAcc                                                                                                                               >  7:394050/1‑149 (MQ=255)                aaCATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCAccc                                                                                                                              <  8:260751/149‑1 (MQ=255)                 aCATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGCGCCACGGGGGCCCAccag                                                                                                                             <  1:307183/149‑3 (MQ=255)                 aCATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGCGCCACGGGGGCCCAccag                                                                                                                             <  5:82805/149‑3 (MQ=255)                          cGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCgg                                                                                                                    <  3:10070/149‑1 (MQ=255)                           gggCACGGACGGCCACTCGTACCTGTCGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGt                                                                                                                   >  6:273839/1‑149 (MQ=255)                              cACGGACGGCCACTCGTACGTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTtgt                                                                                                                <  8:226923/149‑1 (MQ=255)                               acggacggCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTtgtg                                                                                                               <  4:227592/149‑1 (MQ=255)                                 ggacggCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGa                                                                                                             >  5:45099/1‑149 (MQ=255)                                  gacggCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGAt                                                                                                            <  3:417267/149‑1 (MQ=255)                                  gacggCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGAt                                                                                                            <  1:425991/149‑1 (MQ=255)                                   acggCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATg                                                                                                           <  1:99891/149‑1 (MQ=255)                                                    ttGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTg                                                                                           <  8:394050/148‑1 (MQ=255)                                                       cgcTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATc                                                                                        <  6:109400/148‑1 (MQ=255)                                                        gcTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAg                                                                                      <  1:37953/149‑1 (MQ=255)                                                        gcTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAg                                                                                      >  6:113567/1‑149 (MQ=255)                                                        gcTTCCACAACGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAg                                                                                      >  1:41304/1‑149 (MQ=255)                                                            ccACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCGGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCt                                                                                   >  1:100935/1‑148 (MQ=255)                                                             cacaCCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTc                                                                                  >  4:183544/1‑148 (MQ=255)                                                              acacCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCaa                                                                                >  7:135423/1‑149 (MQ=255)                                                               cacCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAAt                                                                               >  2:187794/1‑149 (MQ=255)                                                                 ccgccgCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATcc                                                                             <  4:112940/149‑1 (MQ=255)                                                                     cgcgGCTGTGCTCGTCCTGTACAGGCACCGTTAGCTGCCTGGGGTAAGCTAATGGCTGGCCCCAATCCATAAATGGTGGGCCCCCGTGGCGAACCCTAGCCAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGc                                                                         <  1:320008/149‑1 (MQ=255)                                                                           tctGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTc                                                                   >  3:278753/1‑149 (MQ=255)                                                                              gCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGc                                                                >  2:359538/1‑149 (MQ=255)                                                                               cTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGc                                                                >  4:387144/1‑148 (MQ=255)                                                                                  gTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCa                                                              <  1:237492/147‑1 (MQ=255)                                                                                     cTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGat                                                                                                      <  4:7559/104‑1 (MQ=255)                                                                                     cTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGat                                                                                                      >  3:7559/1‑104 (MQ=255)                                                                                        gACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACg                                                        <  6:45099/147‑1 (MQ=255)                                                                                           aGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGaa                                                   >  6:386720/1‑149 (MQ=255)                                                                                           aGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGaa                                                   <  6:197396/149‑1 (MQ=255)                                                                                             gCACGGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACg                                                 <  1:127048/149‑1 (MQ=255)                                                                                             gCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAAc                                                  >  4:83175/1‑148 (MQ=255)                                                                                             gCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAAc                                                  >  3:168422/1‑148 (MQ=255)                                                                                                       cTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTg                                       >  2:357274/1‑149 (MQ=255)                                                                                                        tGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTgg                                      <  5:386720/149‑1 (MQ=255)                                                                                                          ccTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGc                                    >  7:363260/1‑149 (MQ=255)                                                                                                            tGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAg                                              <  3:72969/137‑1 (MQ=255)                                                                                                            tGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAg                                              >  4:72969/1‑137 (MQ=255)                                                                                                                   aGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTgg                            <  1:111320/148‑1 (MQ=255)                                                                                                                   aGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTgg                            >  2:111320/1‑148 (MQ=255)                                                                                                                     cGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACg                                                        <  6:270573/118‑1 (MQ=255)                                                                                                                     cGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACg                                                        >  5:270573/1‑118 (MQ=255)                                                                                                                     cGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCgg                         <  8:135423/149‑1 (MQ=255)                                                                                                                     cGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCgg                         <  5:113567/149‑1 (MQ=255)                                                                                                                      gAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCggg                        <  3:83175/149‑1 (MQ=255)                                                                                                                       aaTGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGt                       >  6:9612/1‑149 (MQ=255)                                                                                                                         tggctggcCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGa                     <  5:273839/149‑1 (MQ=255)                                                                                                                                 ccccATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCACGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGa             <  8:297293/149‑1 (MQ=255)                                                                                                                                 ccccATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGAcc                   <  8:360181/143‑1 (MQ=255)                                                                                                                                 ccccATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGAcc                   >  7:360181/1‑143 (MQ=255)                                                                                                                                 ccccATCCATCACTGCTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCAAAACACGTTGAACGAAGCGTGTTGGGCAGAACTGGCGGGTGACCCGCt               >  6:331710/1‑147 (MQ=255)                                                                                                                                   ccatccatCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGAGc           <  5:9612/149‑1 (MQ=255)                                                                                                                                      tccatcACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGAGCttt        <  2:41304/149‑1 (MQ=255)                                                                                                                                       ccatcACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGAGCttt        >  5:128677/1‑148 (MQ=255)                                                                                                                                         ctctcgggtgcGCCACGGGGGCCCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGAGCTTTCaa     <  7:12273/143‑1 (MQ=255)                                                                                                                                         atcACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGAGCTTTCaa     >  6:58192/1‑149 (MQ=255)                                                                                                                                          tcACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTTTGGATGCAGATATCCGGCTGTGATCAGCCTCAATCCACGCACTGTCCGCCAATAACGTTTAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTTAGCTTTCAAc    <  6:377992/149‑1 (MQ=255)                                                                                                                                          tcACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGAGCTTTCAAc    <  4:278753/149‑1 (MQ=255)                                                                                                                                          tcACTGGTGGGCCCCCGTGGCGCACACGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGAGCTTTCAAc    >  5:312302/1‑149 (MQ=255)                                                                                                                                            tcgggtgcGCCACGGGGGCCCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGAGCTTTCaa     >  8:12273/4‑146 (MQ=255)                                                                                                                                            aCTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGAGCTTTCAACtt  <  1:187794/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                              CACTTGAGCCACGCCAACATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGAGCTTTCAACTT  >  NC_002947/3835834‑3836121

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