Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,984,700	G→T	6.7%	A140A (GCC→GCA)	ycaC‑II ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,984,700	0	G	T	6.7%	209.7 / 3.1	68	A140A (GCC→GCA)	ycaC‑II	putative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (34/29);  new base T (2/3);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTC  >  NC_002947/5984553‑5984843                                                                                                                                                    |                                                                                                                                                    cGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAgg                                                                                                                                                    <  8:162427/149‑1 (MQ=255) cGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAg                                                                                                                                                     <  2:34935/148‑1 (MQ=255)       cagGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCgg                                                                                                                                              >  7:267556/1‑149 (MQ=255)       cagGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCgg                                                                                                                                              >  7:7629/1‑149 (MQ=255)         ggTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAgg                                                                                                                                                    <  1:281687/141‑1 (MQ=255)         ggTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAgg                                                                                                                                                    >  2:281687/1‑141 (MQ=255)         ggTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGt                                                                                                                                             >  3:67934/1‑148 (MQ=255)           ttACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGACCGATGCCTCGGGcac                                                                                                                                          <  7:192221/149‑1 (MQ=255)             aCCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcacc                                                                                                                                         <  2:308912/148‑1 (MQ=255)              cccAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcaccac                                                                                                                                       <  1:21710/149‑1 (MQ=255)               ccAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcaccac                                                                                                                                       >  4:78968/1‑148 (MQ=255)               ccAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcaccac                                                                                                                                       >  4:292781/1‑148 (MQ=255)               ccAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcaccac                                                                                                                                       >  2:209046/1‑148 (MQ=255)                     cTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCTTCGGTCACCACGAACACt                                                                                                                                <  6:132760/149‑1 (MQ=255)                     cTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACt                                                                                                                                >  8:194508/1‑149 (MQ=255)                     cTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACt                                                                                                                                >  2:164154/1‑149 (MQ=255)                         gATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCa                                                                                                                             >  8:1733/1‑148 (MQ=255)                             tCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAAcc                                                                                                                        >  4:97083/1‑149 (MQ=255)                               gTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCtt                                                                                                                      >  8:163629/1‑149 (MQ=255)                               gTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACACGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCtt                                                                                                                      >  4:125221/1‑149 (MQ=255)                                    gcCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGc                                                                                                                >  5:325791/1‑148 (MQ=255)                                            cgGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGa                                                                                                     <  6:160202/149‑1 (MQ=255)                                                           gCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTGCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAgac                                                                                      <  4:92204/149‑3 (MQ=255)                                                            ccACCGAGACCCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAac                                                                                                                                   <  6:199695/107‑1 (MQ=255)                                                            ccACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGcc                                                                                     <  4:67934/149‑1 (MQ=255)                                                                      ccAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACAt                                                                           >  2:82210/1‑149 (MQ=255)                                                                       cAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTATCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGCGACCGGAAGGCCACGCACACATc                                                                          >  2:97042/1‑149 (MQ=255)                                                                        aGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCg                                                                         >  1:91687/1‑149 (MQ=255)                                                                            catcatCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTca                                                                     >  1:21513/1‑149 (MQ=255)                                                                             atcatcTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTcac                                                                    <  3:97083/149‑1 (MQ=255)                                                                                      gcACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGc                                                           <  6:325791/149‑1 (MQ=255)                                                                                        aCCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCga                                                         >  5:102803/1‑149 (MQ=255)                                                                                        aCCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCga                                                         >  6:31093/1‑149 (MQ=255)                                                                                          ccGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CTCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGAACACATCGGTCACCACCCAGGCcat                                                        >  8:302375/1‑148 (MQ=37)                                                                                          ccGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCgat                                                        >  4:188037/1‑148 (MQ=255)                                                                                                 ggTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAg                                                  <  1:164154/147‑1 (MQ=255)                                                                                                 ggTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAg                                                  <  3:188037/147‑1 (MQ=255)                                                                                                 ggTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGtt                                                <  3:125221/149‑1 (MQ=255)                                                                                                  gTCATTGGTACCCAATCCGCTTGCTGCACGGTTTCGTTGAATGGGCCCGAGGCATGGGTCACCAAGTACAATTCAAAAACTTCTTCGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGt                                                 <  1:97042/149‑1 (MQ=255)                                                                                                  gTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCg                                                             >  8:79917/1‑135 (MQ=255)                                                                                                  gTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCg                                                             <  7:79917/135‑1 (MQ=255)                                                                                                  gTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGttg                                               <  8:7629/149‑1 (MQ=255)                                                                                                  gTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGA‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGttg                                               <  2:279899/149‑1 (MQ=255)                                                                                                          tACCCAAGCCGCTTGCTGCAAGGTTTCGTTGAACGTGCCCGATGCATCGGTCAACAAGAACACTTCAAAACCTTC‑‑CGCCATCCAGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCg                                         <  1:260585/149‑1 (MQ=38)                                                                                                             ccAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCg                                    >  2:9925/1‑149 (MQ=255)                                                                                                              cAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGttacttg                                           <  4:77039/141‑1 (MQ=255)                                                                                                              cAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGttacttg                                           >  3:77039/1‑141 (MQ=255)                                                                                                              cAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGt                                   >  2:305521/1‑149 (MQ=255)                                                                                                              cAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGt                                   <  7:1733/149‑1 (MQ=255)                                                                                                               aaGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTg                                  <  3:78968/149‑1 (MQ=255)                                                                                                               aaGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTg                                  >  3:96883/1‑149 (MQ=255)                                                                                                                            cACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGc                       <  6:102803/147‑1 (MQ=255)                                                                                                                            cACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCt                     >  5:117880/1‑149 (MQ=255)                                                                                                                            cACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCt                     <  6:117880/149‑1 (MQ=255)                                                                                                                            cACGGTTTCGTTGAAGCTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATCATCAGTTGCTTGCGGCACGTGGCCTTGATTGCCt                     >  8:291567/1‑149 (MQ=255)                                                                                                                                gTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGAc                 <  7:194508/149‑1 (MQ=255)                                                                                                                                    cGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGaa              <  2:253051/148‑1 (MQ=255)                                                                                                                                      ttGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑AGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCTGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTc           >  8:3036/1‑149 (MQ=255)                                                                                                                                        gAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCtt         >  6:140296/1‑149 (MQ=255)                                                                                                                                        gAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGAGGCCTTGATTGCCTTGACGAAGTCtt         >  5:16452/1‑149 (MQ=255)                                                                                                                                         aaGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCAGACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTc        <  5:31093/149‑1 (MQ=255)                                                                                                                                         aaGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCtt         <  6:16452/148‑1 (MQ=255)                                                                                                                                          tggTGACCGATGCCTCGGGCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCg       <  1:37409/148‑1 (MQ=255)                                                                                                                                          tggTGACCGATGCCTCGGGCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCg       >  2:37409/2‑149 (MQ=255)                                                                                                                                           ggTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTc                                                                      >  4:184639/1‑85 (MQ=255)                                                                                                                                           ggTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTc                                                                      <  3:184639/85‑1 (MQ=255)                                                                                                                                            gTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CTCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGtt     <  7:163629/149‑1 (MQ=255)                                                                                                                                            gTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGtt     <  8:267556/149‑1 (MQ=255)                                                                                                                                             tgaccGATGCCTCGGGCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTg    >  1:193290/4‑149 (MQ=255)                                                                                                                                             tGCCCGAGGCATTGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCTGGAAGGCCACGCACACATCGGTCACCCCACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTg    >  3:48220/1‑149 (MQ=255)                                                                                                                                               cccGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTc  >  2:149175/1‑149 (MQ=255)                                                                                                                                               cccGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTc  <  7:139871/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTC  >  NC_002947/5984553‑5984843

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