Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,866,023	G→T	8.9%	R171R (CGC→CGA)	PP_4277 ←	GntR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,866,023	0	G	T	8.9%	206.6 / 7.0	66	R171R (CGC→CGA)	PP_4277	GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (29/31);  new base T (2/4);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

CCTCTGGCCAGTGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGTCGAGCTG  >  NC_002947/4865880‑4866169                                                                                                                                                |                                                                                                                                                   ccTCTGGCCAGTGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGTCGCTg                                                                                                                                               >  8:331664/1‑149 (MQ=255) ccTCTGGCCAGTGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTg                                                                                                                                               <  2:166923/149‑1 (MQ=255)       gCCAGTGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAg                                                                                                                                         >  3:87978/1‑149 (MQ=255)          actGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTCAGCACCTCGCCAATGCGCGGCAGCTCCCGCGAGACCACGCCTTCCAGCAAGCGGCCGAACAGGTCGACATAGTCGATCATGTTAGCCTCGCGCTGGTCCCGCGa                                                                                                                                      >  8:197592/3‑149 (MQ=255)          aGTGCCGCCAGGACCAGTTCGCAGCGGCGCTGGCGATAGGCGGTGAGCACCTCGCGATTGCGCGGCAGGTCAGGCGAGACCACGGCTGCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGa                                                                                                                                      <  7:197592/149‑1 (MQ=255)            tGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGTCGCTGGACCAGCGCCa                                                                                                                                    <  4:332690/149‑1 (MQ=255)            tGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGTCGCTGGACCAGCGCCa                                                                                                                                    <  7:331664/149‑1 (MQ=255)              ccgccAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCTGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCTGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGTCGCTGGACCAGCGCCAgg                                                                                                                                  <  2:73607/149‑1 (MQ=255)              ccgccAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGTCGCTGGACCAGCGCCAg                                                                                                                                   <  3:317081/148‑1 (MQ=255)              ccgccAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAgg                                                                                                                                  >  5:303377/1‑149 (MQ=255)              ccgccAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAgg                                                                                                                                  >  5:163503/1‑149 (MQ=255)                     cACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAgg                                                                                                                           >  3:153360/1‑149 (MQ=255)                               cAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCAt                                                                                                                 >  4:150339/1‑149 (MQ=255)                                aGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGAGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGTCACGGCTCATg                                                                                                                >  8:225696/1‑149 (MQ=255)                                       cTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCt                                                                                                         >  6:36044/1‑149 (MQ=255)                                       cTGGCCAGAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCt                                                                                                         <  6:165023/149‑1 (MQ=255)                                           ccATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAg                                                                                                     <  1:48554/149‑1 (MQ=255)                                           ccATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAg                                                                                                     <  1:43053/149‑1 (MQ=255)                                                    gTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCAc                                                                                            <  3:297294/149‑1 (MQ=255)                                                    gTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGGCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCAc                                                                                            <  3:292262/149‑1 (MQ=255)                                                     tGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACg                                                                                           <  5:36044/149‑1 (MQ=255)                                                      gAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACGACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCt                                                                                                         <  3:328023/134‑1 (MQ=255)                                                      gAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCt                                                                                                         >  4:328023/1‑134 (MQ=255)                                                         cACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGttt                                                                                       >  8:294818/1‑149 (MQ=255)                                                          acgtCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTc                                                                                      <  5:164241/146‑1 (MQ=255)                                                          aCCTCGCGAATGCGCGGCAGGTCACGGGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTc                                                                                      <  8:335409/149‑1 (MQ=255)                                                          aCCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTc                                                                                      <  6:303377/149‑1 (MQ=255)                                                          aCCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTc                                                                                      <  6:163503/149‑1 (MQ=255)                                                          aCCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTc                                                                                      <  1:180170/149‑1 (MQ=255)                                                          aCCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTc                                                                                      <  3:324443/149‑1 (MQ=255)                                                             tCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTg                                                                                   >  8:107557/1‑149 (MQ=255)                                                              cgcgAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGACCAAGCGGGCGAACAGGTCGAGTTAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGTCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGc                                                                                  <  7:57057/149‑1 (MQ=255)                                                                cgAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAg                                                                                <  4:200256/149‑1 (MQ=255)                                                                 gAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAg                                                                                >  5:126018/1‑148 (MQ=255)                                                                      cgcgGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTc                                                                                                                   >  5:158013/1‑108 (MQ=255)                                                                      cgcgGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTc                                                                                                                   <  6:158013/108‑1 (MQ=255)                                                                      cgcgGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACgg                                                                          <  4:87978/149‑1 (MQ=255)                                                                       gcgGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACggg                                                                         <  6:166683/149‑1 (MQ=255)                                                                       gcgGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACggg                                                                         <  3:150339/149‑1 (MQ=255)                                                                             ggTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCg                                                                   <  4:153360/149‑1 (MQ=255)                                                                              gTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTACACGGTTTCCTGCCGGAAGGGGTTGTCgg                                                                  >  8:179985/1‑149 (MQ=255)                                                                               tCACGCGAGATCACGGCTTCGAGCAAGCGGGCGAACAGGTCTAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGc                                                                 <  6:207876/149‑1 (MQ=255)                                                                               tCACGCGAGACCACGACTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGc                                                                 <  6:4353/149‑1 (MQ=255)                                                                                 aCGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCg                                                                >  7:45590/1‑148 (MQ=255)                                                                                   gcgAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAAc                                                             >  3:306712/1‑149 (MQ=255)                                                                                       gaCCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGAt                                                         <  7:37752/149‑1 (MQ=255)                                                                                         ccACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAg                                                       >  2:124193/1‑149 (MQ=255)                                                                                         ccACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAg                                                       >  2:321317/1‑149 (MQ=255)                                                                                          cACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAgg                                                      >  8:2243/1‑149 (MQ=255)                                                                                                           gggCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGCTTTACCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGc                                      >  5:85249/1‑148 (MQ=255)                                                                                                                aaCAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACt                                <  8:45590/149‑1 (MQ=255)                                                                                                                    ggTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGTCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCg                            <  7:179985/149‑1 (MQ=255)                                                                                                                     gTCGAGATAGTCGATCATGTTGGCCTGTCGCTGTTCCAGCGACAGGTAATAGGCACCGCTCAAGGCCGCATGCAGGATTTCCACTGTTTCCTGAACGAACGCGTTGTCGGCGAACGGATAGTCCGCACGCATGACAGCAAAACTgactt                           >  7:247087/1‑145 (MQ=255)                                                                                                                     gTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCgg                           >  1:54044/1‑149 (MQ=255)                                                                                                                     gTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCgg                           >  4:372129/1‑149 (MQ=255)                                                                                                                               tCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATa                                                        >  1:34092/1‑110 (MQ=255)                                                                                                                               tCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATa                                                        <  2:34092/110‑1 (MQ=255)                                                                                                                                 gATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCt               >  4:69196/1‑149 (MQ=255)                                                                                                                                    cATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGa            <  4:306712/149‑1 (MQ=255)                                                                                                                                      tGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATg          >  4:9142/1‑149 (MQ=255)                                                                                                                                       gTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAg                   >  8:219887/1‑139 (MQ=255)                                                                                                                                       gTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAg                   <  7:219887/139‑1 (MQ=255)                                                                                                                                       gTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATg          >  6:332875/1‑148 (MQ=255)                                                                                                                                        ttGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGTc        <  1:124193/149‑1 (MQ=255)                                                                                                                                          ggCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGTc        >  2:355904/1‑147 (MQ=255)                                                                                                                                           gCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGTCGAg     >  6:210255/1‑149 (MQ=255)                                                                                                                                                gcgcTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGTCGAGCTg  <  7:107557/147‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   CCTCTGGCCAGTGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGTCGAGCTG  >  NC_002947/4865880‑4866169

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