Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,970,611	G→C	7.5%	intergenic (‑369/+25)	fliC ← / ← PP_4379	flagellin/3‑oxoacyl‑ACP synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,970,611	0	G	C	7.5%	201.0 / 4.5	65	intergenic (‑369/+25)	fliC/PP_4379	flagellin/3‑oxoacyl‑ACP synthase
Reads supporting (aligned to +/- strand):  ref base G (32/28);  new base C (3/2);  total (35/30)
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.60e-01

TCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAACACTGGAG  >  NC_002947/4970463‑4970753                                                                                                                                                      |                                                                                                                                               tCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGccgccg                                                                                                                                                <  3:39495/149‑1 (MQ=255)     tGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGgcg                                                                                                                                            <  4:88859/149‑1 (MQ=255)       gAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGgcgc                                                                                                                                           <  6:252016/148‑1 (MQ=255)       gAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCt                                                                                                                                          >  3:212724/1‑149 (MQ=255)             gTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAgg                                                                                                                                    <  8:124170/149‑1 (MQ=255)               gAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAgggg                                                                                                                                  >  1:146107/1‑149 (MQ=255)                          gaagGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGa                                                                                                                        >  7:1302/1‑148 (MQ=255)                            agGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCa                                                                                                                     <  6:280854/149‑1 (MQ=255)                             ggAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAg                                                                                                                    >  6:295100/1‑149 (MQ=255)                               agaCGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTc                                                                                                                  <  8:102592/149‑1 (MQ=255)                                  cgcgATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTcgcg                                                                                                               <  2:21472/149‑1 (MQ=255)                                  cgcgATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTcgcg                                                                                                               <  5:295100/149‑1 (MQ=255)                                   gcgATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCgg                                                                                                              >  6:298389/1‑149 (MQ=255)                                       tGGGTGGCACGGGGTTTGCGGGTGTTCGGGGGTTAACCCGCTCCTACATTGGGCTGCGTTAGGAATAAGCTTATTTTCGGAAT‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTAc                                                                                                           <  2:146107/148‑1 (MQ=255)                                               aCCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTgg                                                                                                                          <  6:253538/125‑1 (MQ=255)                                               aCCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTgg                                                                                                                          >  5:253538/1‑125 (MQ=255)                                               aCCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGcc                                                                                                  >  3:257975/1‑149 (MQ=255)                                               aCCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGcc                                                                                                  >  6:84731/1‑149 (MQ=255)                                                  ggCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAg                                                                                               <  4:101901/149‑1 (MQ=255)                                                  ggCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAg                                                                                               >  3:101901/1‑149 (MQ=255)                                                     tttGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGccc                                                                                            <  8:287532/149‑1 (MQ=255)                                                        gCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCAc                                                                                         >  4:15551/1‑149 (MQ=255)                                                         ccGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCAc                                                                                         >  4:97263/1‑148 (MQ=255)                                                         ccGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCAc                                                                                         >  2:210005/1‑148 (MQ=255)                                                               ttCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTcgcg                                                                                                               <  2:221679/120‑1 (MQ=255)                                                               ttCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTcgcg                                                                                                               >  1:221679/1‑120 (MQ=255)                                                                  gcgGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCa                                                                               >  3:184088/1‑149 (MQ=255)                                                                  gcgGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCa                                                                               >  4:182080/1‑149 (MQ=255)                                                                    gggTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACg                                                                             <  4:212724/149‑1 (MQ=255)                                                                       tGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGcc                                                                           >  6:321281/1‑148 (MQ=255)                                                                                  cTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGAt                                                               <  3:97263/149‑1 (MQ=255)                                                                                  cTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCAACGCCAAAGCCACTGAt                                                               >  7:286463/1‑149 (MQ=255)                                                                                      aaTGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGa                                                                >  3:120453/1‑144 (MQ=255)                                                                                      aaTGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGa                                                                <  4:120453/144‑1 (MQ=255)                                                                                          ggCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCAcgc                                                       <  5:84731/149‑1 (MQ=255)                                                                                           gCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCAcgc                                                       >  3:230723/1‑149 (MQ=255)                                                                                           gCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCAcgcg                                                      >  5:149560/1‑149 (MQ=255)                                                                                            cTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACgcgc                                                     >  2:256022/1‑149 (MQ=255)                                                                                                   aGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAg                                               <  3:182080/148‑1 (MQ=255)                                                                                                     gCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTg                                            >  7:203848/1‑149 (MQ=255)                                                                                                        gAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCa                                         <  5:321281/149‑1 (MQ=255)                                                                                                                ttttCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCAt                                 >  1:212759/1‑149 (MQ=255)                                                                                                                      gCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCAcgcggc                                                    >  4:100586/1‑122 (MQ=255)                                                                                                                      gCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCAcgcggc                                                    <  3:100586/124‑3 (MQ=255)                                                                                                                        aG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCt                         >  3:154289/1‑149 (MQ=255)                                                                                                                        aG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCt                         >  7:176642/1‑149 (MQ=255)                                                                                                                            cctaaaaGCGCCCGGCGGCCACCGCCCGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGa                                     <  5:55380/131‑1 (MQ=38)                                                                                                                            cctaaaaGCGCCCGGCGGCCACCGCCCGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGa                                     >  6:55380/4‑134 (MQ=38)                                                                                                                             aaaaaaGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTgcag                     >  3:211381/1‑149 (MQ=255)                                                                                                                                 aaGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTgcagcagc                 <  1:256022/149‑1 (MQ=255)                                                                                                                                  aGCGCCCGGCGGCCACCGCCCGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTgcagcag                  >  8:100564/1‑147 (MQ=255)                                                                                                                                    cgcCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCagcagcagc              <  4:257975/149‑1 (MQ=255)                                                                                                                                    cgcCCGGCGGCCACCGCCCGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCa                                  <  8:129530/129‑1 (MQ=255)                                                                                                                                    cgcCCGGCGGCCACCGCCCGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCa                                  >  7:129530/1‑129 (MQ=255)                                                                                                                                     gcCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCg             <  5:298389/149‑1 (MQ=255)                                                                                                                                     gcCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCg             <  4:184088/149‑1 (MQ=255)                                                                                                                                       cGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTgc                       >  2:76899/1‑137 (MQ=255)                                                                                                                                       cGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTgc                       <  1:76899/137‑1 (MQ=255)                                                                                                                                       cGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCgg            >  7:169557/1‑148 (MQ=255)                                                                                                                                        gggCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGaa          <  3:225646/149‑1 (MQ=255)                                                                                                                                            ggTGGCCGCCGGACGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGGTTCTGCAGCAGCAGCGGAACACt      >  8:96903/1‑149 (MQ=255)                                                                                                                                             gTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTc                                                  >  7:147406/1‑104 (MQ=255)                                                                                                                                             gTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTc                                                  <  8:147406/104‑1 (MQ=255)                                                                                                                                               ggCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAACACTgga   <  8:203848/149‑1 (MQ=255)                                                                                                                                                gccgccGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAACACTggag  >  1:90853/1‑149 (MQ=255)                                                                                                                                                gccgccGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAACACTggag  <  8:176642/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                               TCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAG‑AAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAACACTGGAG  >  NC_002947/4970463‑4970753

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑