Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,928,022	C→T	9.7%	G717D (GGC→GAC)	cheA ←	chemotaxis histidine kinase CheA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,928,022	0	C	T	9.7%	164.9 / 6.5	62	G717D (GGC→GAC)	cheA	chemotaxis histidine kinase CheA
Reads supporting (aligned to +/- strand):  ref base C (24/31);  new base T (3/3);  total (27/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 = 1.00e+00

CTGCCATAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTCGTGCACCTGGCCCT  >  NC_002947/4927875‑4928169                                                                                                                                                    |                                                                                                                                                       cTGCCATAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCGGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGcc                                                                                                                                                       <  5:153820/149‑1 (MQ=255)     cATAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCggg                                                                                                                                                   >  5:195416/1‑149 (MQ=255)     cATAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCggg                                                                                                                                                   >  4:213929/1‑149 (MQ=255)       tAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGt                                                                                                                                                 >  5:235209/1‑149 (MQ=255)       tAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGt                                                                                                                                                 >  6:178905/1‑149 (MQ=255)         aaCACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGc                                                                                                                                               <  1:174776/149‑1 (MQ=255)         aaCACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGc                                                                                                                                               >  2:267660/1‑149 (MQ=255)         aaCACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGc                                                                                                                                               <  2:175513/149‑1 (MQ=255)          acacTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGt                                                                                                                                                 <  2:82462/146‑1 (MQ=255)          acacTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGt                                                                                                                                                 >  1:82462/1‑146 (MQ=255)          acacTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGcc                                                                                                                                              <  4:273165/149‑1 (MQ=255)            acTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGGCATGTCCggggcca                                                                                                                                               >  4:209457/1‑143 (MQ=255)            acTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGGCATGTCCggggcca                                                                                                                                               <  3:209457/146‑4 (MQ=255)            acTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGGCATGTCCggggcca                                                                                                                                               <  7:100120/146‑4 (MQ=255)            acTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGGCATGTCCggggcca                                                                                                                                               >  8:100120/1‑143 (MQ=255)            acTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCtt                                                                                                                                            >  8:276976/1‑149 (MQ=255)             cTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTg                                                                                                                                           <  4:296987/149‑1 (MQ=255)              tCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTg                                                                                                                                           >  4:263895/1‑148 (MQ=255)               ccATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTgca                                                                                                                                         <  5:83818/149‑1 (MQ=255)               ccATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTgca                                                                                                                                         <  3:5582/149‑1 (MQ=255)                cATTAGTCGGCGCAACCCGCCAGGGCGAGCCACCGAAATCAAATACGCCGCGCGCCGTAACGCTTGAGAATGAACGGGACGTAGAGAATCAACGCGATGCCACCGTCGCCGGTAATGGTGCCCACGGACAAACCCGCGGTCCCTTgcag                                                                                                                                        >  2:139052/1‑149 (MQ=255)                cATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTgc                                                                                                                                          <  2:272684/147‑1 (MQ=255)                 aTTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTg                                                                                                                                                >  8:74368/1‑140 (MQ=255)                 aTTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTg                                                                                                                                                <  7:74368/140‑1 (MQ=255)                      gcggcgCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGTCCGTCGCCGGTAATGGTGGCCCCGGACATGCCCg                                                                                                                                                     <  5:331525/130‑1 (MQ=255)                      gcggcgCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCg                                                                                                                                                     >  6:331525/1‑130 (MQ=255)                         gcgcAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGcc                                                                                                                               <  6:105034/149‑1 (MQ=255)                                acgacAGGCCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCAAGCTCCGGACGTCGAGAATCAACGCGATGCGACCGTCTGCGGTAATGTTGGCCCAGGACATGCCCGCGCAGCCATGCAGCATCTTGACCagcagc                                                                                                                        >  8:98697/5‑149 (MQ=255)                                     aGGGCGCGCCACCTAAATCTAATACGCCGCGAGGCGTAACGATTGTGCATGCTCGGGACGTTGAGAATCAAGGCGATGCGACCGTCGCCGGTTATGGTGGCCCCGGAGATGGCCGGGGTGCCTTGCAGAATCTTGCCCAGCGGCTtgat                                                                                                                   <  1:139052/149‑1 (MQ=255)                                       ggCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTtgatga                                                                                                                 >  3:310701/1‑149 (MQ=255)                                        gcgcgcCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGAc                                                                                                                <  3:79184/149‑1 (MQ=255)                                                 cGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACt‑tc‑ttc                                                                                                     >  5:129715/1‑149 (MQ=255)                                                 cGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACt‑tc‑ttc                                                                                                     <  2:172680/149‑1 (MQ=255)                                                 cGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACt‑tc‑ttc                                                                                                     <  5:178905/149‑1 (MQ=255)                                                        aaaTACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGccc                                                                                              <  4:318440/149‑1 (MQ=255)                                                         gttACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCcca                                                                                             >  1:203761/3‑149 (MQ=255)                                                          aTACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGTCCCGGACATGCCCGGGGTGCCTTgca                                                                                                                                         >  3:113352/1‑106 (MQ=255)                                                          aTACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGTCCCGGACATGCCCGGGGTGCCTTgca                                                                                                                                         <  4:113352/106‑1 (MQ=255)                                                             cgccgcGCGGCGTAACGCTTGAGGATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGTTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCccac‑ca                                                                                         <  8:23741/148‑1 (MQ=255)                                                              gccgcGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCggg                                                                                                                                                   <  2:287957/92‑1 (MQ=255)                                                              gccgcGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCggg                                                                                                                                                   >  1:287957/1‑92 (MQ=255)                                                              gccgcGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGt                                                                                       >  6:119522/1‑149 (MQ=255)                                                                cgcgcgGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTg                                                                                     <  4:310701/149‑1 (MQ=255)                                                                      gcgTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTcgac                                                                               <  3:213929/149‑1 (MQ=255)                                                                      gcgTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTcgac                                                                               <  6:235209/149‑1 (MQ=255)                                                                         tAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTcgacga                                                                             >  1:113266/1‑148 (MQ=255)                                                                            cGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACaaa                                                                         <  3:263895/149‑1 (MQ=255)                                                                             gCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACaaaa                                                                        <  2:203761/149‑1 (MQ=255)                                                                                          cGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTggg                                                           >  4:12682/1‑149 (MQ=255)                                                                                           gggACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTgg                                                            <  6:195416/147‑1 (MQ=255)                                                                                             gACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTGGGTGc                                                        <  4:137016/149‑1 (MQ=255)                                                                                                gTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTGGGTGccga                                                     <  7:276976/149‑1 (MQ=255)                                                                                                      aaTCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACa                                               <  2:202131/149‑1 (MQ=255)                                                                                                            cgcgATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATca                                         <  3:12682/149‑1 (MQ=255)                                                                                                             gcgATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATcac                                        >  7:24364/1‑149 (MQ=255)                                                                                                                                    aggcaccccGGGCATGTCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCtt                 >  5:33062/6‑149 (MQ=255)                                                                                                                                     ggcaccccGGGCATGTCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGAAGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTc                <  6:33062/145‑1 (MQ=255)                                                                                                                                        gCCCCGGAGATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGTCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTCGt              <  7:98697/148‑1 (MQ=255)                                                                                                                                        gCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTCTCTTTCCTGGCCCCCACCGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTc                <  1:267660/149‑1 (MQ=255)                                                                                                                                        gCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCCC‑CAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTCGTg             >  1:220580/1‑149 (MQ=255)                                                                                                                                          cccGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCCGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTCGTGCa           >  1:262928/1‑149 (MQ=255)                                                                                                                                               aCATGCCCGGGGTGTCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTCGTGCACCTgg      <  7:155490/149‑1 (MQ=255)                                                                                                                                                   gCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTCGTGCACCTGGCCCt  >  2:222101/1‑149 (MQ=255)                                                                                                                                                    cccGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTCGTGCACCTGGCCCt  >  7:12063/1‑148 (MQ=255)                                                                                                                                                    cccGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTCGTGCACCTGGCCCt  <  8:12063/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                       CTGCCATAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTCGTGCACCTGGCCCT  >  NC_002947/4927875‑4928169

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