Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,984,703	A→C	5.8%	D139E (GAT→GAG)	ycaC‑II ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,984,703	0	A	C	5.8%	214.4 / 3.4	69	D139E (GAT→GAG)	ycaC‑II	putative hydrolase
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base C (2/2);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

CAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAG  >  NC_002947/5984559‑5984847                                                                                                                                                 |                                                                                                                                                      cagGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCgg                                                                                                                                                   >  7:267556/1‑149 (MQ=255) cagGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCgg                                                                                                                                                   >  7:7629/1‑149 (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:292781/1‑148 (MQ=255)         ccAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcaccac                                                                                                                                            >  2:209046/1‑148 (MQ=255)         ccAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcaccac                                                                                                                                            >  4:78968/1‑148 (MQ=255)               cTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCTTCGGTCACCACGAACACt                                                                                                                                     <  6:132760/149‑1 (MQ=255)               cTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACt                                                                                                                                     >  2:164154/1‑149 (MQ=255)               cTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACt                                                                                                                                     >  8:194508/1‑149 (MQ=255)                   gATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCa                                                                                                                                 >  8:1733/1‑148 (MQ=255)                       tCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAAcc                                                                                                                            >  4:97083/1‑149 (MQ=255)                         gTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCtt                                                                                                                          >  8:163629/1‑149 (MQ=255)                         gTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACACGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCtt                                                                                                                          >  4:125221/1‑149 (MQ=255)                              gcCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGc                                                                                                                    >  5:325791/1‑148 (MQ=255)                                      cgGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGa                                                                                                         <  6:160202/149‑1 (MQ=255)                                                     gCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTGCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAgac                                                                                          <  4:92204/149‑3 (MQ=255)                                                      ccACCGAGACCCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAac                                                                                                                                        <  6:199695/107‑1 (MQ=255)                                                      ccACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGcc                                                                                         <  4:67934/149‑1 (MQ=255)                                                                ccAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACAt                                                                               >  2:82210/1‑149 (MQ=255)                                                                 cAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTATCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGCGACCGGAAGGCCACGCACACATc                                                                              >  2:97042/1‑149 (MQ=255)                                                                  aGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCg                                                                             >  1:91687/1‑149 (MQ=255)                                                                      catcatCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTca                                                                         >  1:21513/1‑149 (MQ=255)                                                                       atcatcTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTcac                                                                        <  3:97083/149‑1 (MQ=255)                                                                                gcACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGc                                                               <  6:325791/149‑1 (MQ=255)                                                                                  aCCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCga                                                             >  6:31093/1‑149 (MQ=255)                                                                                  aCCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCga                                                             >  5:102803/1‑149 (MQ=255)                                                                                    ccGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CTCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGAACACATCGGTCACCACCCAGGCcat                                                            >  8:302375/1‑148 (MQ=37)                                                                                    ccGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCgat                                                            >  4:188037/1‑148 (MQ=255)                                                                                           ggTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAg                                                      <  3:188037/147‑1 (MQ=255)                                                                                           ggTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAg                                                      <  1:164154/147‑1 (MQ=255)                                                                                           ggTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGtt                                                    <  3:125221/149‑1 (MQ=255)                                                                                            gTCATTGGTACCCAATCCGCTTGCTGCACGGTTTCGTTGAATGGGCCCGAGGCATGGGTCACCAAGTACAAT‑TCAAAAACTTCTTCGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGt                                                     <  1:97042/149‑1 (MQ=255)                                                                                            gTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCg                                                                 <  7:79917/135‑1 (MQ=255)                                                                                            gTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCg                                                                 >  8:79917/1‑135 (MQ=255)                                                                                            gTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGttg                                                   <  8:7629/149‑1 (MQ=255)                                                                                            gTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGA‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGttg                                                   <  2:279899/149‑1 (MQ=255)                                                                                                    tACCCAAGCCGCTTGCTGCAAGGTTTCGTTGAACGTGCCCGATGCATCGGTCAACAAGAACACT‑TCAAAACCTTC‑‑CGCCATCCAGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCg                                             <  1:260585/149‑1 (MQ=38)                                                                                                       ccAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCg                                        >  2:9925/1‑149 (MQ=255)                                                                                                        cAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGttacttg                                               >  3:77039/1‑141 (MQ=255)                                                                                                        cAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGttacttg                                               <  4:77039/141‑1 (MQ=255)                                                                                                        cAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGt                                       >  2:305521/1‑149 (MQ=255)                                                                                                        cAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGt                                       <  7:1733/149‑1 (MQ=255)                                                                                                         aaGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTg                                      >  3:96883/1‑149 (MQ=255)                                                                                                         aaGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTg                                      <  3:78968/149‑1 (MQ=255)                                                                                                                      cACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGc                           <  6:102803/147‑1 (MQ=255)                                                                                                                      cACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCt                         <  6:117880/149‑1 (MQ=255)                                                                                                                      cACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCt                         >  5:117880/1‑149 (MQ=255)                                                                                                                      cACGGTTTCGTTGAAGCTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATCATCAGTTGCTTGCGGCACGTGGCCTTGATTGCCt                         >  8:291567/1‑149 (MQ=255)                                                                                                                          gTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGAc                     <  7:194508/149‑1 (MQ=255)                                                                                                                              cGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGaa                  <  2:253051/148‑1 (MQ=255)                                                                                                                                ttGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑AGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCTGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTc               >  8:3036/1‑149 (MQ=255)                                                                                                                                  gAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCtt             >  6:140296/1‑149 (MQ=255)                                                                                                                                  gAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGAGGCCTTGATTGCCTTGACGAAGTCtt             >  5:16452/1‑149 (MQ=255)                                                                                                                                   aaGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCAGACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTc            <  5:31093/149‑1 (MQ=255)                                                                                                                                   aaGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCtt             <  6:16452/148‑1 (MQ=255)                                                                                                                                    tggTGACCGATGCCTCGGGCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCg           >  2:37409/2‑149 (MQ=255)                                                                                                                                    tggTGACCGATGCCTCGGGCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCg           <  1:37409/148‑1 (MQ=255)                                                                                                                                     ggTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTc                                                                          >  4:184639/1‑85 (MQ=255)                                                                                                                                     ggTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTc                                                                          <  3:184639/85‑1 (MQ=255)                                                                                                                                      gTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CTCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGtt         <  7:163629/149‑1 (MQ=255)                                                                                                                                      gTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGtt         <  8:267556/149‑1 (MQ=255)                                                                                                                                       tgaccGATGCCTCGGGCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTg        >  1:193290/4‑149 (MQ=255)                                                                                                                                       tGCCCGAGGCATTGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCTGGAAGGCCACGCACACATCGGTCACCCCACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTg        >  3:48220/1‑149 (MQ=255)                                                                                                                                         cccGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTc      <  7:139871/149‑1 (MQ=255)                                                                                                                                         cccGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTc      >  2:149175/1‑149 (MQ=255)                                                                                                                                             atgCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAg  <  7:302375/147‑1 (MQ=255)                                                                                                                                                 aTCGGTCACCACGAACACTCTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAg                                                                                                      <  1:174365/46‑1 (MQ=255)                                                                                                                                                 aTCGGTCACCACGAACACTCTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAg                                                                                                      <  1:169193/46‑1 (MQ=255)                                                                                                                                                 aTCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAg                                                                                                      >  2:169193/1‑45 (MQ=255)                                                                                                                                                 aTCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAg                                                                                                      >  2:174365/1‑45 (MQ=255)                                                                                                                                                 |                                                                                                                                                      CAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAG  >  NC_002947/5984559‑5984847

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