Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,502,985	A→T	6.7%	S94C (AGC→TGC)	speB →	agmatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,502,985	0	A	T	6.7%	174.3 / 3.0	60	S94C (AGC→TGC)	speB	agmatinase
Reads supporting (aligned to +/- strand):  ref base A (33/23);  new base T (2/2);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

CCAGTAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGC  >  NC_002947/2502842‑2503133                                                                                                                                                 |                                                                                                                                                     ccaGTAACCACCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCgg                                                                                                                                                  >  5:177278/1‑148 (MQ=255)   aGTAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCAc                                                                                                                                               >  2:339785/1‑149 (MQ=255)     tAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACTGCGGcgtgc                                                                                                                                             >  8:68657/1‑145 (MQ=255)       aCCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCAcgccg                                                                                                                                           >  2:271564/1‑149 (MQ=255)          gCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAg                                                                                                                                        >  8:2499/1‑149 (MQ=255)             ccGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACTGCGGCGTGCCGCTGTCg                                                                                                                                     <  4:46344/149‑1 (MQ=255)             ccGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCg                                                                                                                                     <  4:23666/149‑1 (MQ=255)                gcgcgcgcTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACTGCGGCGTGCCGCTGTCgaag                                                                                                                                  <  7:68657/149‑4 (MQ=255)                     cgcTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGa                                                                                                                             <  8:294197/149‑1 (MQ=255)                       cTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACa                                                                                                                           <  4:181075/149‑1 (MQ=255)                           ggCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCa                                                                                                                        <  1:271564/148‑1 (MQ=255)                            gCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATc                                                                                                                      >  5:87895/1‑149 (MQ=255)                              ccGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTCGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCAcgccg                                                                                                                                           <  1:82325/127‑1 (MQ=255)                              ccGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCAcgccg                                                                                                                                           >  2:82325/1‑126 (MQ=255)                               cgcgGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGc                                                                                                                         >  8:216739/1‑143 (MQ=255)                               cgcgGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGc                                                                                                                         <  7:216739/143‑1 (MQ=255)                                 cgGGCGATCCGTGCTGCCTCGGTCCTGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGc                                                                                                                 <  7:227335/149‑1 (MQ=255)                                 cgGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGc                                                                                                                 <  1:365476/149‑1 (MQ=255)                                        tCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGcc                                                                                                          >  1:357116/1‑149 (MQ=255)                                           gTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAg                                                                                                       <  5:53763/149‑1 (MQ=255)                                                  cTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTc                                                                                                >  8:258826/1‑149 (MQ=255)                                                      gTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTg                                                                                              >  6:271746/1‑147 (MQ=255)                                                        ccAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGCCCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACTGCGGCGTGCCGCTGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAgg                                                                                          >  2:76894/1‑149 (MQ=255)                                                             agGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGt                                                                                     >  6:186252/1‑149 (MQ=255)                                                               gCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTg                                                                                   >  2:362828/1‑149 (MQ=255)                                                                ccTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTgc                                                                                  >  4:111191/1‑149 (MQ=255)                                                                 cTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCg                                                                                                         >  5:113462/1‑125 (MQ=255)                                                                 cTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCg                                                                                                         <  6:113462/125‑1 (MQ=255)                                                                      cccGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATg                                                                            <  1:339785/149‑1 (MQ=255)                                                                        cGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCt                                                                          >  7:27566/1‑149 (MQ=255)                                                                          ccACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGa                                                                        <  2:357116/149‑1 (MQ=255)                                                                           cACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGAc                                                                       <  5:376837/149‑1 (MQ=255)                                                                                   gTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtggtg                                                               <  8:346059/149‑1 (MQ=255)                                                                                              acccgttcgaccACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCAGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACAATATGCAGGCGGGTTGAGCCATTCTGACCCTAGGTTGCGACCACtt                                                     >  7:179864/1‑148 (MQ=255)                                                                                              acccgttcgaccACCT‑GGCGGACATCGACTATGGCGACTGTCCCTTCGACAGCGGAACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATACCGAGCACATTCTGCATGCGGGTTGCGCCATGCTGACCCTTGGTCGCGACCAcat                                                     >  6:174282/1‑146 (MQ=255)                                                                                               cccgttcgaccACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtca                                                   >  4:36007/1‑149 (MQ=255)                                                                                                ccgttcgaccACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtcat                                                  >  2:277056/1‑149 (MQ=255)                                                                                                  gttcgaccACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACNGCATCGAGNNCCNTGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtcatca                                                <  7:2499/149‑1 (MQ=255)                                                                                                  gttcgaccACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtcatca                                                <  3:111191/149‑1 (MQ=255)                                                                                                  gttcgaccACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtcat                                                  >  6:11157/1‑147 (MQ=255)                                                                                                  gttcgaccACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtcat                                                  >  6:16453/1‑147 (MQ=255)                                                                                                  gttcgaccACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtcat                                                  >  3:46294/1‑147 (MQ=255)                                                                                                   ttcgaccACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAg                                               <  6:177278/149‑1 (MQ=255)                                                                                                                    tCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAg                               >  5:87203/1‑149 (MQ=255)                                                                                                                     cATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAgg                              >  8:214761/1‑149 (MQ=255)                                                                                                                      aTCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGc                             >  8:233666/1‑149 (MQ=255)                                                                                                                       tCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGc                             >  8:233899/1‑148 (MQ=255)                                                                                                                          aCTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAc                         <  6:87895/149‑1 (MQ=255)                                                                                                                               ggCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAcgcccg                    >  3:104203/1‑149 (MQ=255)                                                                                                                                 cGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAcgcccgcc                  <  6:87203/149‑1 (MQ=255)                                                                                                                                 cGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAcgcccgcc                  <  5:186252/149‑1 (MQ=255)                                                                                                                                    cTGTCCCTTCGACAGCGGCACGGAGCAGTCGGTACAGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCcgccgc                >  7:55827/1‑148 (MQ=255)                                                                                                                                         ccTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATgg           <  3:347704/148‑1 (MQ=255)                                                                                                                                         ccTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGc          <  5:11157/149‑1 (MQ=255)                                                                                                                                         ccTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGc          >  2:305372/1‑149 (MQ=255)                                                                                                                                         ccTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGc          >  5:33920/1‑149 (MQ=255)                                                                                                                                         ccTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGc          <  5:16453/149‑1 (MQ=255)                                                                                                                                          cTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGcc         >  6:245211/1‑149 (MQ=255)                                                                                                                                               aCAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTg    >  6:325861/1‑149 (MQ=255)                                                                                                                                               aCAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTg    >  6:118380/1‑149 (MQ=255)                                                                                                                                                 aGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTATCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGgc  <  1:277056/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CCAGTAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCT‑GGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGC  >  NC_002947/2502842‑2503133

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