Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,557,922	T→G	6.4%	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	6.4%	165.4 / 2.2	59	intergenic (‑293/+38)	PP_4888/purA	methyl‑accepting chemotaxis transducer/adenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base T (26/29);  new base G (2/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

CGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCACGC  >  NC_002947/5557775‑5558068                                                                                                                                                     |                                                                                                                                                   cGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGAACAGACATGCCCACAAATGAAAACGCCGCGGCCCTtc                                                                                                                                                   >  6:4882/1‑149 (MQ=255)  gCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTtct                                                                                                                                                  <  4:329152/149‑1 (MQ=255)   cTCCTACACCGGCCTGGGCCCATCCATCGTGTGGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTtctc                                                                                                                                                 <  2:6197/149‑1 (MQ=38)      cTACAACGGCCTGCGCCCAACCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAgg                                                                                                                                              <  4:169477/149‑1 (MQ=255)      cTACAACGGCCTGCGCCCAACCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAgg                                                                                                                                              >  7:168532/1‑149 (MQ=255)      cTACAACGGCCTGCGCCCAACCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAgg                                                                                                                                              >  6:360341/1‑149 (MQ=255)         cAACGGCCTGCGCCCAACCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAGATGAAAACGCCGCGGCCCTTCTCAGGGcc                                                                                                                                           >  3:50668/1‑149 (MQ=255)           aCGGCCTGCGCCCAACCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCcgc                                                                                                                                         <  4:86054/149‑1 (MQ=255)            cGGCCTGCGCCCAACCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCcgcg                                                                                                                                        <  3:260711/149‑1 (MQ=255)            cGGCCTGCGCCCAACCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCcgcg                                                                                                                                        <  6:195856/149‑1 (MQ=255)              gCCTGCGCCCAACCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCgcggc                                                                                                                                      <  2:234798/149‑1 (MQ=255)               ccTGCGCCCAACCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCgcggc                                                                                                                                      <  2:52587/148‑1 (MQ=255)               ccTGCGCCCAACCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCgcggc                                                                                                                                      >  1:52587/1‑148 (MQ=255)                cctcccccAACCCTTGTTGTGGGCCCCATATTGTCGGGGTTTAACGCGCTAAAAGGGTAGGGGCGGCAACGGTGGATTACCGGTGTATTTTGGCCACAACTCACATCCCAACAAATCAAAACGCCGCGGCCCTTCTCAGGGCCgcggcg                                                                                                                                     <  8:118055/144‑1 (MQ=255)                 tGCGCCCAACCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAGATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTc                                                                                                                                   <  4:50668/149‑1 (MQ=255)                   cgcCCAACCATCGTTGTGTG‑CCCACTTAGGAGCGGGTTTACCCGCGAATAGGGTAGCGGCTGCAACGGTGGATTACGGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTc                                                                                                                                   <  5:194966/147‑1 (MQ=255)                    gcCCAACCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAg                                                                                                                                 >  2:13502/1‑148 (MQ=255)                      ccAACCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGAACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGAt                                                                                                                               >  7:195864/1‑148 (MQ=255)                        aaCCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCCGGGCCGCGGCGTCAGATCAg                                                                                                                            >  6:319229/1‑149 (MQ=255)                           cATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCt                                                                                                                         <  6:159605/149‑1 (MQ=255)                           cATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCt                                                                                                                         >  2:92030/1‑149 (MQ=255)                                ttGGGTG‑ACCACATAGGAGCTGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTTTATTACCGGTGGATATTGGCCAGCACAGACATGCCCAAAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAg                                                                                                                    <  8:168532/149‑1 (MQ=38)                                   tgtg‑CCCAGTTTGGTGGGGGTTTACCCGCTATTACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCGGCGGCGTCAGATCAGCCTGGAAGGc                                                                                                                  <  4:251404/148‑1 (MQ=38)                                        cccACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCa                                                                                                                 <  3:126233/145‑1 (MQ=255)                                        cccACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCa                                                                                                                 >  4:126233/1‑145 (MQ=255)                                          cacaTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGAt                                                                                                              >  8:79374/1‑146 (MQ=255)                                          cacaTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCa                                                                                                            <  7:79374/148‑1 (MQ=255)                                           acaTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCCGGGCCGGGGCGTCAGATCAGCCTGGAAGGCAGATCAgg                                                                                                          <  3:147962/149‑1 (MQ=255)                                           acaTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAgg                                                                                                          >  8:26662/1‑149 (MQ=255)                                            caTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGc                                                                                                         <  5:339392/149‑1 (MQ=255)                                             aTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGc                                                                                                         >  8:45910/1‑148 (MQ=255)                                               aGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGaa                                                                                                      >  7:76674/1‑149 (MQ=255)                                                  aGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGGCCTTGTGAGGGCCGCGGCGTCAGATAAGCCTGGGAGGCAGCTCAGGCGAACgg                                                                                                   >  1:358985/1‑149 (MQ=38)                                                  aGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACgg                                                                                                   <  5:360341/149‑1 (MQ=255)                                                    cGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGt                                                                                                 <  1:13502/149‑1 (MQ=255)                                                    cGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGt                                                                                                 <  7:26662/149‑1 (MQ=255)                                                             ccGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGGGCACAAATGAAAACGCCGCGGCCGTTGTCAGGGGCGCGGCGGCAGATCCGCCTGGAAGGAAGAACAGGCCAACGGGTGACgcagcc                                                                                        >  2:277635/1‑148 (MQ=38)                                                                    tACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCAAAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGt                                                                                 <  7:45910/149‑1 (MQ=255)                                                                      cgttAGGGGCGGCAAGGGGGTATTCCGGTTGTTTATGGCCCACCACTCACATCCCCAAAAATTAAACCGCCGGGGCCCTGCTCTGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGtc                                                                                <  2:358985/145‑1 (MQ=255)                                                                           gcggcgGCAACGGTTTATTACCGGTGGATATTGGCCAGCACAGACATTCCCATAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGtt                                                                           <  1:296436/148‑1 (MQ=255)                                                                              gcggcAACGGTGGATTACCGTTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACTCAGCACGATGGTCTCGTTGCgg                                                                       <  1:277635/149‑1 (MQ=255)                                                                                  cAACGGTTGATTACGGGTGTATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCg                                                                   <  8:343981/149‑1 (MQ=255)                                                                                               ccGGTGGATATTGCCCCCCACACACATGCCCACACATCAACACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGaaa                                                       <  5:4882/148‑1 (MQ=255)                                                                                                  gTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAcatcat                                                   <  8:76674/149‑1 (MQ=255)                                                                                                  gTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTAACAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGTGTGACGCAGCACGATGGGCTCGAGGCGGTCCGCGCCGGGCGAAatg                                                     >  4:97075/1‑147 (MQ=255)                                                                                                       tatTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGACGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCg                                               >  7:162800/1‑148 (MQ=255)                                                                                                            gCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgc                                         >  3:310258/1‑149 (MQ=255)                                                                                                             ccAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgcg                                        >  3:22389/1‑149 (MQ=255)                                                                                                                      aCATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAg                                                                                                                    <  5:20802/64‑1 (MQ=255)                                                                                                                             cacaAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCGGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAatgat                                                   <  1:303716/122‑1 (MQ=255)                                                                                                                             cacaAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCGGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAatgat                                                   >  2:303716/1‑122 (MQ=255)                                                                                                                             cacaAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAg                              >  2:258679/1‑143 (MQ=255)                                                                                                                             cacaAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAg                              <  1:258679/143‑1 (MQ=255)                                                                                                                                     ctgaCGCCGCGGCCCGGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGCCGCAGCAGGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGGGCCAACCAGCTCTTCGATGCGCt                >  3:30690/4‑149 (MQ=255)                                                                                                                                     aaaaCGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATgcgc                 <  1:92030/148‑1 (MQ=255)                                                                                                                                      tgaCGCCGCGGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGAt                     >  7:58618/3‑143 (MQ=255)                                                                                                                                      tgaCGCCGCGGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGAt                     <  8:58618/141‑1 (MQ=255)                                                                                                                                          gccgcGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATg           >  6:185925/1‑149 (MQ=255)                                                                                                                                             gcgGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAg        >  5:160210/1‑149 (MQ=255)                                                                                                                                                 cccTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAc    >  1:328515/1‑149 (MQ=255)                                                                                                                                                  ccTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAcg   >  4:196368/1‑149 (MQ=255)                                                                                                                                                   cTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAcgc  >  4:227910/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTG‑CCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCACGC  >  NC_002947/5557775‑5558068

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