Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,928,017	A→G	9.5%	S719P (TCC→CCC)	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,017	0	A	G	9.5%	165.3 / 7.8	63	S719P (TCC→CCC)	cheA	chemotaxis histidine kinase CheA
Reads supporting (aligned to +/- strand):  ref base A (23/34);  new base G (3/3);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

CCTTGACTGCCATAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTCGTGCACCTGG  >  NC_002947/4927869‑4928165                                                                                                                                                     |                                                                                                                                                        ccTTGACTGCCATAAACACTCCATTAGGCGGCGCAACCCGCCCGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGa                                                                                                                                                         <  8:50661/149‑1 (MQ=255)  cTTGACTGCCATAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGa                                                                                                                                                         >  8:98921/1‑148 (MQ=255)  cTTGACTGCCATAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGa                                                                                                                                                         >  3:137016/1‑148 (MQ=255)  cTTGACTGCCATAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGAc                                                                                                                                                        <  7:98921/149‑1 (MQ=255)     gACTGCCATAAACACTCCATTAGGCGGTGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACAt                                                                                                                                                      <  1:205996/148‑1 (MQ=255)       cTGCCATAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCGGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGcc                                                                                                                                                   <  5:153820/149‑1 (MQ=255)           cATAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCggg                                                                                                                                               >  4:213929/1‑149 (MQ=255)           cATAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCggg                                                                                                                                               >  5:195416/1‑149 (MQ=255)             tAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGt                                                                                                                                             >  5:235209/1‑149 (MQ=255)             tAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGt                                                                                                                                             >  6:178905/1‑149 (MQ=255)               aaCACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGc                                                                                                                                           <  2:175513/149‑1 (MQ=255)               aaCACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGc                                                                                                                                           >  2:267660/1‑149 (MQ=255)               aaCACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGc                                                                                                                                           <  1:174776/149‑1 (MQ=255)                acacTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGt                                                                                                                                             >  1:82462/1‑146 (MQ=255)                acacTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGt                                                                                                                                             <  2:82462/146‑1 (MQ=255)                acacTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGcc                                                                                                                                          <  4:273165/149‑1 (MQ=255)                  acTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGGCATGTCCggggcca                                                                                                                                           <  7:100120/146‑4 (MQ=255)                  acTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGGCATGTCCggggcca                                                                                                                                           >  8:100120/1‑143 (MQ=255)                  acTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGGCATGTCCggggcca                                                                                                                                           >  4:209457/1‑143 (MQ=255)                  acTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGGCATGTCCggggcca                                                                                                                                           <  3:209457/146‑4 (MQ=255)                  acTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCtt                                                                                                                                        >  8:276976/1‑149 (MQ=255)                   cTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTg                                                                                                                                       <  4:296987/149‑1 (MQ=255)                    tCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTg                                                                                                                                       >  4:263895/1‑148 (MQ=255)                     ccATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTgca                                                                                                                                     <  3:5582/149‑1 (MQ=255)                     ccATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTgca                                                                                                                                     <  5:83818/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                                                                                                 <  5:178905/149‑1 (MQ=255)                                                       cGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACt‑tc‑ttc                                                                                                 <  2:172680/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                                                                                                                                               >  1:287957/1‑92 (MQ=255)                                                                    gccgcGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCggg                                                                                                                                               <  2:287957/92‑1 (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                                                                           <  6:235209/149‑1 (MQ=255)                                                                            gcgTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTcgac                                                                           <  3:213929/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)                                                                                                                                                     |                                                                                                                                                        CCTTGACTGCCATAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGGGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACT‑TC‑TTCCTGGCCCAC‑CAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTCGTGCACCTGG  >  NC_002947/4927869‑4928165

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