Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,502,992	C→T	9.7%	T96M (ACG→ATG) ‡	speB →	agmatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,502,992	0	C	T	9.7%	202.0 / 8.4	72	T96M (ACG→ATG) ‡	speB	agmatinase
Reads supporting (aligned to +/- strand):  ref base C (38/27);  new base T (4/3);  total (42/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-01

CCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGCGCTGA  >  NC_002947/2502852‑2503138                                                                                                                                             |                                                                                                                                                   ccccGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGt                                                                                                                                            <  3:249135/149‑1 (MQ=255)   ccGGCGCGCGCCTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCg                                                                                                                                          >  4:409921/1‑149 (MQ=255)    cGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCgg                                                                                                                                         <  1:427534/149‑1 (MQ=255)         cgcgcTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACTGCGGCGTGCCGCTGTCgaaggga                                                                                                                                    <  5:39140/149‑7 (MQ=255)            gcTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCACTTCGACAGCCGCACGCCGCACTCTGTACCTGa                                                                                                                                  >  5:117994/1‑148 (MQ=255)              ttCGGCCCGCGGGCGATCCGTACTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAg                                                                                                                               >  5:184583/1‑149 (MQ=255)                cGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGc                                                                                                                              >  2:37625/1‑148 (MQ=255)                     cgcgGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCCGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAg                                                                                                                        <  4:245581/149‑1 (MQ=255)                      gcgGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAgg                                                                                                                       <  7:5454/149‑1 (MQ=255)                       cgGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGc                                                                                                                      <  6:79833/149‑1 (MQ=255)                                  tgctgcCTCGGTCCAGCAGGCCTGGGGCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAg                                                                                                            >  7:117281/1‑148 (MQ=255)                                    ctgcATCGGTACAGCAGGCCTGTGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGca                                                                                                          <  6:117994/148‑1 (MQ=255)                                          cGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTg                                                                                                   >  5:17463/1‑149 (MQ=255)                                           ggTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGc                                                                                                  <  6:184583/149‑1 (MQ=255)                                              ccAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAgg                                                                                               >  8:309041/1‑149 (MQ=255)                                               cagcagGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGc                                                                                              >  2:123178/1‑149 (MQ=255)                                                agcagGCCTGGGCCCGCCAATGGCCGTGGGCGTTCGACCCGTTCGACAACCTGGCGGTCATCGACTATGGAGACTGTCCCTTCGAAAGCGGCACGCCGCAGTCGGTACCGGACAGAATCGAGGCACAAGCCTAGCACATTCTTAAGGCg                                                                                             >  6:227390/1‑149 (MQ=255)                                                 gcagGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCgg                                                                                            <  1:123178/149‑1 (MQ=255)                                                 gcagGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCgg                                                                                            <  4:74327/149‑1 (MQ=255)                                                       cTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTgcg                                                                                      >  6:250427/1‑149 (MQ=255)                                                       cTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTgcg                                                                                      >  6:314297/1‑149 (MQ=255)                                                          ggCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCTGGTTGCGCCa                                                                                   <  5:66028/149‑1 (MQ=255)                                                               gCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCgg                                                                                            >  6:330666/1‑135 (MQ=255)                                                               gCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCgg                                                                                            <  5:330666/135‑1 (MQ=255)                                                               gCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCt                                                                               >  8:410728/1‑148 (MQ=255)                                                                  aCTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGAATGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGAcc                                                                           <  4:347053/149‑1 (MQ=255)                                                                   cTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGAccc                                                                          >  7:155447/1‑149 (MQ=255)                                                                    tGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCt                                                                         >  5:403008/1‑149 (MQ=255)                                                                    tGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCt                                                                         >  8:100073/1‑149 (MQ=255)                                                                      gCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtg                                                                       <  7:289118/149‑1 (MQ=255)                                                                       ccGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtgg                                                                      <  8:117281/149‑1 (MQ=255)                                                                           gggCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtgg                                                                      >  3:153521/1‑145 (MQ=255)                                                                           gggCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtgg                                                                      <  4:153521/145‑1 (MQ=255)                                                                           gggCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGc                                                                  <  6:403008/149‑1 (MQ=255)                                                                               gttcgacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGAc                                                               >  6:52642/1‑148 (MQ=255)                                                                               gttcgacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGAc                                                               >  1:249342/1‑148 (MQ=255)                                                                                ttcgacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCa                                                             >  3:27950/1‑149 (MQ=255)                                                                                    acccgttcgaccACCTGGCGGTCATCGACTATGGAGACTGTCCCTTCGACAGCTGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGTCGGGTTTCGCCATGCTGACCCTTGGTGGCGACCACTtc                                                         <  5:227390/149‑1 (MQ=255)                                                                                        gttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGTACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGa                                                                <  5:290688/138‑1 (MQ=255)                                                                                        gttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGTACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGa                                                                >  6:290688/1‑138 (MQ=255)                                                                                         gatgaccGCCAGGTGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAg                                                    >  7:383090/4‑149 (MQ=255)                                                                                            gaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTa                                                 >  3:299323/1‑149 (MQ=255)                                                                                               caccTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGAGCCATGCTGACCCTAGGTGGCGACAACTTCATCAGCTAcac                                              >  6:403928/1‑147 (MQ=255)                                                                                                              gACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTc                                <  7:410728/148‑1 (MQ=255)                                                                                                              gACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCa                               >  7:130107/1‑149 (MQ=255)                                                                                                               aCTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAc                              >  8:179384/1‑149 (MQ=255)                                                                                                                     gCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAcgcccgc                        <  5:413237/149‑1 (MQ=255)                                                                                                                        aCTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCcgccgc                     <  4:299323/149‑1 (MQ=255)                                                                                                                         cTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCAACCCGCTGCTCAAGGCTCACGCCCgccgcc                    <  1:37625/149‑1 (MQ=255)                                                                                                                          tGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCa                   >  4:271128/1‑149 (MQ=255)                                                                                                                          tGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCa                   >  5:209467/1‑149 (MQ=255)                                                                                                                            tCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATg                 >  7:280469/1‑149 (MQ=255)                                                                                                                              ccTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGc               <  5:52642/149‑1 (MQ=255)                                                                                                                               cTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCACCGCTCACGCCCGCCGCCATGGcc              >  3:104873/1‑149 (MQ=255)                                                                                                                                ttCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGAc                                                                            >  4:59648/1‑86 (MQ=255)                                                                                                                                ttCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGAc                                                                            <  3:59648/86‑1 (MQ=255)                                                                                                                                  ccacacCGGCACGCCGCCGCCGGCACCGCACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGccc             <  3:319065/145‑1 (MQ=255)                                                                                                                                  cGACTGCGGCGTGCCGCTGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATg                 <  3:267124/143‑1 (MQ=255)                                                                                                                                  cGACTGCGGCGTGCCGCTGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATg                 >  4:267124/1‑143 (MQ=255)                                                                                                                                  cGACTGCGGCGTGCCGCTGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGc           >  4:115320/1‑149 (MQ=255)                                                                                                                                  cGACAGCGGCACGCCGCAGTCGGTACCGGACAGGATCGAGGCCCATGCCGAGCAGAGTCTGGAGGCGCGTTGCGGCGGGCTGGCCCGGGGTGGCGAGCGCGTCGTGGGCGGGCCGCTGCGGCCGGCTCAGGCCGGCCGCCAGGGcgcgg           >  8:313415/1‑148 (MQ=255)                                                                                                                                   gACTGCGGCGTGCCGCTGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCt          >  1:352174/1‑149 (MQ=255)                                                                                                                                    actgcggcgTGCCGCTGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGc               >  3:326193/4‑143 (MQ=255)                                                                                                                                    actgcggcgTGCCGCTGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGc               <  4:326193/140‑1 (MQ=255)                                                                                                                                    aCAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTg         >  3:20587/1‑149 (MQ=255)                                                                                                                                      aGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGgc       >  1:426811/1‑149 (MQ=255)                                                                                                                                        cGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGgcgc     >  8:117124/1‑149 (MQ=255)                                                                                                                                        cGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGgcgc     >  7:325660/1‑149 (MQ=255)                                                                                                                                         ggCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGCGCt    <  4:20587/149‑1 (MQ=255)                                                                                                                                           cACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGCGCTGa  >  5:133669/1‑149 (MQ=255)                                                                                                                                           cACGCCGCAGTCGGTACCGGACAGCATCGAGG‑CCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTc                                >  2:423454/1‑118 (MQ=255)                                                                                                                                           cACGCCGCAGTCGGTACCGGACAGCATCGAGG‑CCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTc                                <  1:423454/118‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                   CCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGCGCTGA  >  NC_002947/2502852‑2503138

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