Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,986,048	C→G	7.3%	P97R (CCG→CGG)	flgA →	flagella basal body P‑ring formation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,986,048	0	C	G	7.3%	219.4 / 7.2	68	P97R (CCG→CGG)	flgA	flagella basal body P‑ring formation protein
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base G (2/3);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

GGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGG‑TGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGT  >  NC_002947/4985900‑4986192                                                                                                                                                     |                                                                                                                                                   gggTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGcc                                                                                                                                                    >  4:356457/1‑149 (MQ=255)      tCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTg                                                                                                                                                <  2:286877/148‑1 (MQ=255)          gAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTggg                                                                                                                                              <  1:198945/146‑1 (MQ=255)          gAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGc                                                                                                                                             <  6:414111/147‑1 (MQ=255)             ttCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCGGCTGGGCCgggc                                                                                                                                        <  2:271964/149‑4 (MQ=255)             ttCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCgtgt                                                                                                                                        >  6:124532/1‑149 (MQ=255)             ttCACCGTTGAAGAATATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCgtgt                                                                                                                                        <  3:211883/149‑1 (MQ=255)                 ccGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACg                                                                                                                                    <  5:180130/149‑1 (MQ=255)                    ttGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTa                                                                                                                                 >  7:44956/1‑149 (MQ=255)                       aaGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACg                                                                                                                               <  2:227350/148‑1 (MQ=255)                       aaGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACg                                                                                                                               >  1:227350/1‑148 (MQ=255)                             aTCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCGGCTGGGCCgg                                                                                                                                          >  6:63154/1‑130 (MQ=255)                             aTCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCGGCTGGGCCgg                                                                                                                                          <  5:63154/131‑2 (MQ=255)                              tCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTg                                                                                                                                                <  8:365735/124‑1 (MQ=255)                              tCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTg                                                                                                                                                >  7:365735/1‑124 (MQ=255)                                    caccacCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACg                                                                                                                >  4:100519/1‑149 (MQ=255)                                      ccaccCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTACGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGgc                                                                                                              <  3:68027/149‑1 (MQ=255)                                       caccCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGgcg                                                                                                             <  5:124532/149‑1 (MQ=255)                                        accCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGgcgc                                                                                                            <  8:153488/149‑1 (MQ=255)                                          ccAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCgcag                                                                                                          >  2:212843/1‑149 (MQ=255)                                            aGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGcc                                                                                                                                                    <  6:1388/106‑1 (MQ=255)                                            aGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGcc                                                                                                                                                    >  5:1388/1‑106 (MQ=255)                                               cTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCg                                                                                                     <  5:93894/149‑1 (MQ=255)                                               cTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCg                                                                                                     <  5:88016/149‑1 (MQ=255)                                                     gccgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACg                                                                                               >  5:62680/1‑149 (MQ=255)                                                      ccgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACTGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACgg                                                                                              >  7:246924/1‑149 (MQ=255)                                                       cgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGt                                                                                             <  4:83387/149‑1 (MQ=255)                                                        gCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTa                                                                                            >  8:251809/1‑149 (MQ=255)                                                          tATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTAt                                                                                           >  8:264529/1‑148 (MQ=255)                                                            tGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTc                                                                                         <  3:100519/148‑1 (MQ=255)                                                            tGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTc                                                                                         <  6:62680/148‑1 (MQ=255)                                                            tGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCg                                                                                        >  4:92886/1‑149 (MQ=255)                                                            tGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCg                                                                                        >  3:33802/1‑149 (MQ=255)                                                              aaaTCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTg                                                                                      >  3:9912/1‑149 (MQ=255)                                                                    aGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCAGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGcc                                                                                >  6:424607/1‑149 (MQ=255)                                                                       tCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTg                                                                                      >  5:302974/1‑140 (MQ=255)                                                                       tCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTg                                                                                      <  6:302974/140‑1 (MQ=255)                                                                        cAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCg                                                                            >  2:247687/1‑149 (MQ=255)                                                                         aaCCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGt                                                                           >  4:95765/1‑149 (MQ=255)                                                                         aaCCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGt                                                                           >  4:90247/1‑149 (MQ=255)                                                                            ccGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGCGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGc                                                                         <  1:106990/149‑1 (MQ=255)                                                                            ccGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCg                                                                                                     <  7:14270/120‑1 (MQ=255)                                                                            ccGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCg                                                                                                     >  8:14270/1‑120 (MQ=255)                                                                              gCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGAc                                                                      >  4:160958/1‑149 (MQ=255)                                                                               cTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACt                                                                     <  3:356457/149‑1 (MQ=255)                                                                                    cccACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTcc                                                                <  8:44956/149‑1 (MQ=255)                                                                                     ccACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCcg                                                               >  1:165140/1‑149 (MQ=255)                                                                                     ccACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCcg                                                               <  7:251809/149‑1 (MQ=255)                                                                                       aCGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGCCTTTTCcgcg                                                             >  1:149605/1‑149 (MQ=255)                                                                                       aCGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCAAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGAGCCGGCCACCGTGCGACTTTTCcgc                                                              >  4:270867/1‑148 (MQ=255)                                                                                               gcATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtgg‑tg                                                    <  3:92886/149‑1 (MQ=255)                                                                                                  tGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtgg‑tggtg                                                 <  1:247687/149‑1 (MQ=255)                                                                                                    ccGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtgg‑tggtggt                                               <  7:264529/149‑1 (MQ=255)                                                                                                    ccGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtgg‑tggtgg                                                <  2:165140/148‑1 (MQ=255)                                                                                                                 agcagcTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGG‑TGGTGGTGACCCGCCCGCt                                   >  5:278073/1‑148 (MQ=255)                                                                                                                      cTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTgg                                                                                                  >  7:223557/1‑81 (MQ=255)                                                                                                                      cTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTgg                                                                                                  <  8:223557/81‑1 (MQ=255)                                                                                                                         gACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTTGGTGGTGACCCGCCCGCTCAAgcgc                            <  1:212843/148‑1 (MQ=255)                                                                                                                         gACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGG‑TGGTGGTGAcccgcc                                       <  6:394582/136‑1 (MQ=255)                                                                                                                         gACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGG‑TGGTGGTGAcccgcc                                       >  5:394582/1‑136 (MQ=255)                                                                                                                         gACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGG‑TGGTGGTGACCCGCCCGCTCAAgcgc                            <  1:362546/147‑1 (MQ=255)                                                                                                                              cTCACTTGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTTCG‑GATACGATGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCAACCTTGCGACTGTTCAGCGCAGTTG‑TGGGGGTGACCCGCCCGCGCAAGCCCGacaaca                     >  5:293777/1‑149 (MQ=255)                                                                                                                               tCACTGGAAAGTCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGG‑TGGTGGTGACCCGCCCGCTCAAGCGCGACAacac                    >  6:259621/1‑149 (MQ=255)                                                                                                                               tCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGG‑TGGTGGTGACCCGCCCGCTCAAGCGCGACAacac                    >  7:149885/1‑149 (MQ=255)                                                                                                                                cACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGG‑TGGTGGTGACCCGCCCGCTCAAGCGCGACAACACg                   >  4:395909/1‑149 (MQ=255)                                                                                                                                 aCTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGG‑TGGTGGTGACCCGCCCGCTCAAGCGCGACAACACg                   <  4:9912/148‑1 (MQ=255)                                                                                                                                     gAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGG‑TGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTTgg              >  1:214397/1‑149 (MQ=255)                                                                                                                                       tcacacGGCCCAGCGGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGG‑TGGTGGTGAcccgcc                                       >  5:2716/6‑122 (MQ=255)                                                                                                                                       tcacacGGCCCAGCGGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGG‑TGGTGGTGAcccgcc                                       <  6:2716/117‑1 (MQ=255)                                                                                                                                         gCCCGGCCCAGCCGGTGGGCCGTGGGCCGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCCTGCCGGCCACCGTGCGACTTTTCCGCGACGTGG‑TGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGaa          <  4:33802/149‑1 (MQ=255)                                                                                                                                                 cAGCCGCTGGGCGGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGG‑TGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGt  <  5:424607/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACG‑GATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGG‑TGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGT  >  NC_002947/4985900‑4986192

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