Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,974,358	T→C	6.4%	K416E (AAG→GAG)	flgK ←	flagellar hook‑associated protein FlgK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,974,358	0	T	C	6.4%	184.8 / 2.3	62	K416E (AAG→GAG)	flgK	flagellar hook‑associated protein FlgK
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base C (2/2);  total (32/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 = 5.66e-01

TTCGCCCAGCGTAGGTTGGGTGATGGTGCCGGTACCGCTGTTGCCACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCACACTGTACTGGCTGGCGCTGGTGAACT  >  NC_002947/4974212‑4974504                                                                                                                                                   |                                                                                                                                                     ttCGCCCAGCGTAGGTTGGGTGATGGTGCCGGTACCGCTGTTGCCACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTg                                                                                                                                                     <  8:234752/148‑1 (MQ=255)    gCCCAGCGTAGGTTGGGTGATGGTGCCGGTACCGCTGTTGCCACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGc                                                                                                                                                 >  4:43525/1‑149 (MQ=255)          cGTAGGTTGGGTGATGGTGCCGGTACCGCTGTTGCCACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTctc                                                                                                                                           <  4:173152/149‑3 (MQ=255)            tAGGTTGGGTGATGGTGCCGGTACCGCTGTTGCCACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCgccg                                                                                                                                         <  1:98637/149‑1 (MQ=255)              ggTTGGGTGATGGTGCCGGTACCGCTGTTGCCACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGGGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCgccggc                                                                                                                                       <  5:65853/149‑1 (MQ=255)                   ggTGATGGTGCCGGTACCGCTGTTGCCACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGgccgc                                                                                                                                    >  4:263063/1‑147 (MQ=255)                   ggTGATGGTGCCGGTACCGCTGTTGCCACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCa                                                                                                                                  <  3:38893/149‑1 (MQ=255)                     tGATGGTGTCGGTACCGCTGTTGCCACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAgg                                                                                                                                >  2:1090/1‑149 (MQ=255)                          gTGCCGGTACCGCTGTTGCCTCTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGgccgc                                                                                                                            <  3:204531/148‑1 (MQ=255)                            gCCGGTACCGCTGTTGCCACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTTAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGgccgc‑cgc                                                                                                                        >  8:168441/1‑149 (MQ=255)                             ccGGTACCGCTGTGGCCACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCAGCGGTCAGGGTGGTGGTGAGGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGgccgc‑cgcc                                                                                                                       <  1:127531/149‑1 (MQ=255)                                gTACCGCTGTTGCCACTGCCTGCGGCAGCACTGATCGGCCCGGCGAAGGCGACGTTGTTGGCATCGGTCAGGGGGGGGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCCCCTTGACGCAGTCGCCGGCCGCCAGGCCGC‑CGCCttt                                                                                                                    <  1:231686/149‑1 (MQ=255)                                gTACCGCTGTTGCCACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCttt                                                                                                                    >  7:208731/1‑149 (MQ=255)                                    cGCTGTTGCCACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAAc                                                                                                                <  1:1090/149‑1 (MQ=255)                                    cGCTGTTGCCACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAAc                                                                                                                >  2:96124/1‑149 (MQ=255)                                           gCCACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAgg                                                                                                          <  8:208731/148‑1 (MQ=255)                                            ccACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCt                                                                                                        <  3:263063/149‑1 (MQ=255)                                            ccACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCt                                                                                                        >  4:30861/1‑149 (MQ=255)                                              aCTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTg                                                                                                       >  7:23232/1‑148 (MQ=255)                                              aCTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTg                                                                                                       >  1:230302/1‑148 (MQ=255)                                                      ccGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑t                                                                                               <  7:37966/147‑1 (MQ=255)                                                                 aTCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCCGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCAcc                                                                                       >  5:152551/1‑144 (MQ=255)                                                                 aTCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCCGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCAcc                                                                                       <  6:152551/144‑1 (MQ=255)                                                                    ggCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCgg                                                                               <  3:43525/149‑1 (MQ=255)                                                                    ggCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCgg                                                                               <  2:230302/149‑1 (MQ=255)                                                                    ggCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCgg                                                                               >  4:116679/1‑149 (MQ=255)                                                                      cccGGCGAAGGCCAGTTTGTTGGCATCGGTTAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGt                                                                             <  7:168441/149‑1 (MQ=255)                                                                        cGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGgccgc‑cgcc                                                                                                                       <  6:243321/106‑1 (MQ=255)                                                                        cGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGgccgc‑cgcc                                                                                                                       >  5:243321/1‑106 (MQ=255)                                                                         ggCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGgccgcc                                                                                                                                   <  5:317266/94‑1 (MQ=255)                                                                         ggCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGgccgcc                                                                                                                                   >  6:317266/1‑94 (MQ=255)                                                                           cGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCa                                                                        <  6:112441/149‑1 (MQ=255)                                                                               ggCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCAtt                                                                      >  6:287743/1‑147 (MQ=255)                                                                                  cAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGt                                                                 <  3:30861/149‑1 (MQ=255)                                                                                   aGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTc                                                                <  5:287743/149‑1 (MQ=255)                                                                                            gCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGc                                                        <  4:67437/148‑1 (MQ=255)                                                                                                  gTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTCGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCN‑TCTATCACCGGTGCCGGGTCGGCATTCAgg                                                                  <  1:26921/132‑1 (MQ=255)                                                                                                  gTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTCGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAgg                                                                  >  2:22060/1‑132 (MQ=255)                                                                                                  gTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTCGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAgg                                                                  >  2:26921/1‑132 (MQ=255)                                                                                                  gTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTCGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAgg                                                                  <  1:22060/132‑1 (MQ=255)                                                                                                           gtggtgATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCAcc                                                                                       <  6:47880/102‑1 (MQ=255)                                                                                                           gtggtgATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCAcc                                                                                       >  5:47880/1‑102 (MQ=255)                                                                                                           gtggtgATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTc                                        <  8:23232/149‑1 (MQ=255)                                                                                                            tggtgATGCTGCCGGCTGACGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGACTTTCAAAGCAAGGCTGAAACCG‑TCTATCACCGGTGCCCGGTCGGCATTCAGGTCTAAGCTGAACATGTCGGTGCCGTCt                                       >  8:279432/1‑149 (MQ=255)                                                                                                                      gccggctgccgCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGca                             >  6:189933/1‑149 (MQ=255)                                                                                                                       ccggctgccgCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCCACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGca                             >  1:214912/1‑148 (MQ=255)                                                                                                                       ccggctgccgCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGca                             >  6:118304/1‑148 (MQ=255)                                                                                                                       ccggctgccgCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCGTAAACCGCACTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAgcggcg                               <  1:96124/147‑1 (MQ=255)                                                                                                                             gccgCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCACACTGTa                      >  2:205584/1‑149 (MQ=255)                                                                                                                             gccgCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCACACTGTa                      >  8:166839/1‑149 (MQ=255)                                                                                                                              ccgCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCACACt                         <  8:265602/145‑1 (MQ=255)                                                                                                                              ccgCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCACACt                         >  7:265602/1‑145 (MQ=255)                                                                                                                                     ggTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCCCACTGTAatggctgg              >  7:146240/1‑149 (MQ=255)                                                                                                                                      gTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGTATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCACACTGTACtggctggc             >  5:303407/1‑149 (MQ=255)                                                                                                                                         gggATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCGAGGCTGAAACCG‑TCTATCACCCGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCACACTGTACTGGCTGGCGCt          <  8:291650/149‑1 (MQ=255)                                                                                                                                              cACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCACACTGTACTGGCTGGCGCTGGTg      >  8:193747/1‑148 (MQ=255)                                                                                                                                              cACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCACACTGTACTGGCTGGCGCTGGTg      >  1:18698/1‑148 (MQ=255)                                                                                                                                                ccTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCACACTGTACTGGCTGGCGCTGGTGAAc   >  4:93335/1‑149 (MQ=255)                                                                                                                                                 cTTGAAGCTGTCGCCGGCCGCCAGGCCGCACGCCTTTCAACGCCAGGCTGAAACCG‑TCTGTCACCGGTGCCGGGTAGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCACACTGTActggc                 >  1:25377/1‑135 (MQ=255)                                                                                                                                                 cTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTGTCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCACACTGTActggc                 <  2:25377/134‑1 (MQ=255)                                                                                                                                                 cTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCACACTGTACTGGCTGGCGCTGGTGAACt  >  7:160975/1‑149 (MQ=255)                                                                                                                                                 cTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCACACTGTACTGGCTGGCGCTGGTGAACt  >  5:305698/1‑149 (MQ=255)                                                                                                                                                  ttGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCAtt                                                                      >  7:236998/1‑80 (MQ=255)                                                                                                                                                  ttGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCAtt                                                                      <  8:236998/80‑1 (MQ=255)                                                                                                                                                  ttGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCACACTGTACTGGCTGGCGCTg         <  7:229389/141‑1 (MQ=255)                                                                                                                                                  ttGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCACACTGTACTGGCTGGCGCTg         >  8:229389/1‑141 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TTCGCCCAGCGTAGGTTGGGTGATGGTGCCGGTACCGCTGTTGCCACTGCCTGCCGCAGCACTGATCGGCCCGGCGAAGGCCAGTTTGTTGGCATCGGTCAGGGTGGTGGTGATGCTGCCGGCTGCCGCACGGGTAGGGATCACCTTGAAGCTGTCGCCGGCCGCCAGGCCGC‑CGCCTTTCAACGCCAGGCTGAAACCG‑TCTATCACCGGTGCCGGGTCGGCATTCAGGTCGAAGCTGCCCATGTCGGTGCCGTCTGAGCGGCGCACACTGTACTGGCTGGCGCTGGTGAACT  >  NC_002947/4974212‑4974504

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