Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,557,922	T→G	8.8%	intergenic (‑293/+38)	PP_4888 ← / ← purA	methyl‑accepting chemotaxis transducer/adenylosuccinate synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,557,922	0	T	G	8.8%	183.4 / 7.2	66	intergenic (‑293/+38)	PP_4888/purA	methyl‑accepting chemotaxis transducer/adenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base G (3/3);  total (31/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.27e-01

GCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAC  >  NC_002947/5557776‑5558066                                                                                                                                                   |                                                                                                                                                 gCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGCAATGCCCACAAATGAAAACGCCGCGGCCCAtc                                                                                                                                                 >  4:6042/1‑148 (MQ=255)     cTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAgg                                                                                                                                            >  3:216046/1‑149 (MQ=255)       aCAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGATGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGCTCAGGGc                                                                                                                                          >  4:116918/1‑149 (MQ=255)         aaCGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCcg                                                                                                                                        >  8:331683/1‑149 (MQ=255)             gCCTGCTCCCAACCATCGTTGGGGGCCCACATCGGAGCGGTTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGAGAAGGGCCgcggc                                                                                                                                    <  6:81266/149‑1 (MQ=37)                ggcgcCAAACCATCTTTGTGTCCCCAAATCGGCGGGTTTTTACCGGCTAATACGGTAGCGGCGGCAACGGTGTATTACCGGTGTATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTc                                                                                                                                 <  4:203501/148‑1 (MQ=255)                ggcgcCAAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGAGAAGGGCCGCGGCGtt                                                                                                                                 <  8:277279/148‑2 (MQ=38)                  cgcCCAACCACCGTTGTGGGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCa                                                                                                                                <  5:188453/148‑1 (MQ=255)                     ccAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGAt                                                                                                                             >  2:102473/1‑148 (MQ=255)                        aCCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGc                                                                                                                         >  4:368450/1‑149 (MQ=255)                              gTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGaa                                                                                                                   <  7:331683/149‑1 (MQ=255)                              gTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGaa                                                                                                                   <  1:34956/149‑1 (MQ=255)                               ttGTTTGCCAAAATAGGTGGGGTTTAACCCGCAAATACGGTAGCGGCGGAAACGGTGGATTACCGGTGTATTTTGGCCAGCACAGACATCCCCACAAATCAAAACGCCGCGGCCCTTCTCTGGGCCGCGGCGTCAGATCAGCCTGGaa                                                                                                                   <  2:200665/148‑1 (MQ=255)                                   gtgCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAgg                                                                                                                 >  4:31232/1‑146 (MQ=255)                                       ccACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGc                                                                                                                >  5:82274/1‑143 (MQ=255)                                       ccACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGc                                                                                                                <  6:82274/143‑1 (MQ=255)                                       ccACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGAt                                                                                                            >  5:1232/1‑147 (MQ=255)                                       ccACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGAt                                                                                                            >  8:113151/1‑147 (MQ=255)                                        cacaTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGACCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAg                                                                                                         >  5:30291/1‑149 (MQ=255)                                         acaTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAgg                                                                                                        <  4:216046/149‑1 (MQ=255)                                                cggggttttACCCGCTAATAGGGTAGCGGCGGCAACGGTGTATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACgg                                                                                                 <  6:322334/143‑1 (MQ=255)                                                 gCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTATCACGGCCGCGGAGTAACATAAACCAGGAAGGCAGCTAAGGCGAACggg                                                                                                >  6:59171/1‑149 (MQ=37)                                                  cGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAg                                                                                                              <  7:106301/134‑1 (MQ=255)                                                  cGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAg                                                                                                              >  8:106301/1‑134 (MQ=255)                                                     gTTTACCCTCTAATACGGTAGCTGCGGCAACGGTTGATTACCGGTGTATATTGGCCAGCACAGACATGCCCACAAATGATAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCTGCCTTCCAGGCTGATct                                                                                                          <  8:292448/135‑2 (MQ=255)                                                     gTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCTGCCTTCCAGGCTGATct                                                                                                          >  7:292448/1‑134 (MQ=38)                                                          cccGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACgcag                                                                                        >  7:78393/1‑148 (MQ=255)                                                                 aTACGGTAGGGGCGGAAACGGTGGATTACCGGTGTATATTGGCCACCACATACATGCCCACAAATCAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGAt                                                                                                                             <  6:10359/104‑1 (MQ=255)                                                                 aTACGGTAGCGGCGGCAACGGTGGATTACCGGTGTATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATg                                                                                 <  2:273163/148‑1 (MQ=255)                                                                 aTACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGAt                                                                                                                             >  5:10359/1‑104 (MQ=255)                                                                     ggTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCACGGCCGCGGCGTCAGATCCCCCAGGAAGGCAGCTCAGGCGAACGGGTGACGCAGCACGAAGGtcac                                                                            >  3:367819/1‑147 (MQ=255)                                                                         gcggcgGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTg                                                                        <  1:102473/149‑1 (MQ=255)                                                                          cggcggCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGc                                                                       >  1:21779/1‑149 (MQ=255)                                                                           ggcggcAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCg                                                                      <  5:127770/149‑1 (MQ=255)                                                                            gcggaAACGGTGTATACCCGGTGTATATTGCCCAGAACTGACATGCCAACAAATGAAAACGCCGCGGCCCTTATCCGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCgg                                                                     <  3:6042/149‑1 (MQ=255)                                                                               gcAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTcc                                                                  <  7:113151/149‑1 (MQ=255)                                                                                cAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCg                                                                 <  6:165305/149‑1 (MQ=255)                                                                                cAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGAGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCg                                                                 <  5:59171/149‑1 (MQ=255)                                                                                        gATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCTGCCTTCCAGGCTGATct                                                                                                          >  3:52905/1‑99 (MQ=255)                                                                                          ttACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCg                                                        >  5:297895/1‑148 (MQ=255)                                                                                             ccGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCGGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAat                                                    >  1:266718/1‑149 (MQ=255)                                                                                                gTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAAGGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAatgat                                                 <  3:31232/149‑1 (MQ=255)                                                                                                gTGGATATGGCCCAGCACAGACATCCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAatgat                                                 <  3:38667/149‑1 (MQ=255)                                                                                                   gATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTc                                              >  1:165778/1‑149 (MQ=255)                                                                                                          gCCAGCACAGACATGCCCACAAATGAAAACGCCCCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGATCGACATCATTTCGACCGgc                                       <  3:230880/149‑1 (MQ=255)                                                                                                           ccAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCAGGGCCGGTCGAAATGATGTCGATCGgcg                                      >  4:106196/1‑149 (MQ=255)                                                                                                            cAGCACAGACATGCCCACAAATGAAAACGCCGCGACCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgcgc                                     <  6:30291/149‑1 (MQ=255)                                                                                                                  aGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAAc                               <  7:352521/149‑1 (MQ=255)                                                                                                                    aCATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCa                             <  2:280991/149‑1 (MQ=255)                                                                                                                    aCATGCCAAAAATTCATAACGCCGCGGCCCTTCTCAGGGCGGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGTGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAAcc                              <  4:367819/148‑1 (MQ=255)                                                                                                                      aTGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGc                           >  1:329934/1‑149 (MQ=255)                                                                                                                         cccACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGGCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCCCGATGGGCTCGTTGCGGGCCGGGGCGGTCGAAATGATGTCGCTCGGCGCGGCAACCAGctct                        >  1:361834/1‑149 (MQ=255)                                                                                                                         cccACAAATGAAAACGCCGCGGCCCTTCCCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGGCTCGTTGCGGTCCGGGCCGGGCGAAATGATGTCGATCGGCGCGCCAACCAGctct                        <  3:116918/149‑1 (MQ=255)                                                                                                                          ccACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCa                                 <  3:333947/139‑1 (MQ=255)                                                                                                                          ccACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCa                                 >  4:333947/1‑139 (MQ=255)                                                                                                                            acaAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTc                      >  4:134817/1‑148 (MQ=255)                                                                                                                              aaaTGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGAt                   <  3:134817/149‑1 (MQ=255)                                                                                                                              aaaTGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGAt                   <  3:106196/149‑1 (MQ=255)                                                                                                                                      aCGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGa           <  3:368450/149‑1 (MQ=255)                                                                                                                                      aCGCCGCGGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGa           >  4:123906/1‑149 (MQ=255)                                                                                                                                        gccgcGGCCCTTCTCAGGGCCGCGGCGTAAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATg         >  8:268145/1‑149 (MQ=255)                                                                                                                                        gccgcGGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCtt                       <  8:349864/135‑1 (MQ=255)                                                                                                                                        gccgcGGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCtt                       >  7:349864/1‑135 (MQ=255)                                                                                                                                         ccgcGGCCCTTCTCAGGACCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGt        >  7:275802/1‑149 (MQ=255)                                                                                                                                             ggCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCGTGactacatc    <  8:78393/149‑8 (MQ=255)                                                                                                                                             ggCCCGTCTCAGGGCCGCGGCGTCAGATCAGCCGGGAAGGCAGATCAGGCGAACGGGTGACGCACCACGATGGTATCGTTGCGGTCCGGGGCGGGCGCAATGATGGCGAGCGGCGCGCCCACCCGCTCGTCGATGCGGGTGATGGAGt     >  8:296326/1‑148 (MQ=255)                                                                                                                                               cccTTCTCAGGGCCGCGGCGTAAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAc  <  7:268145/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAC  >  NC_002947/5557776‑5558066

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑