Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,984,708	T→G	5.6%	T138P (ACC→CCC)	ycaC‑II ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,984,708	0	T	G	5.6%	219.7 / 3.0	71	T138P (ACC→CCC)	ycaC‑II	putative hydrolase
Reads supporting (aligned to +/- strand):  ref base T (34/33);  new base G (2/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

GGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTT  >  NC_002947/5984561‑5984852                                                                                                                                                    |                                                                                                                                                      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                                                                                                                                                 >  4:78968/1‑148 (MQ=255)       ccAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcaccac                                                                                                                                                 >  2:209046/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                                                    <  4:77039/141‑1 (MQ=255)                                                                                                      cAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGttacttg                                                    >  3:77039/1‑141 (MQ=255)                                                                                                      cAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGt                                            <  7:1733/149‑1 (MQ=255)                                                                                                      cAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGt                                            >  2:305521/1‑149 (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                              >  5:117880/1‑149 (MQ=255)                                                                                                                    cACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCt                              <  6:117880/149‑1 (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:169193/46‑1 (MQ=255)                                                                                                                                               aTCGGTCACCACGAACACTCTCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAg                                                                                                           <  1:174365/46‑1 (MQ=255)                                                                                                                                               aTCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAg                                                                                                           >  2:169193/1‑45 (MQ=255)                                                                                                                                               aTCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAg                                                                                                           >  2:174365/1‑45 (MQ=255)                                                                                                                                                tCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGtt  >  7:46011/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACT‑TCAAAACCTTC‑‑CGCCAGC‑‑GCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTT  >  NC_002947/5984561‑5984852

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