Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,557,924	T→G	11.4%	intergenic (‑295/+36)	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,924	0	T	G	11.4%	135.1 / 12.0	62	intergenic (‑295/+36)	PP_4888/purA	methyl‑accepting chemotaxis transducer/adenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base T (26/27);  new base G (5/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.80e-01

CTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAC  >  NC_002947/5557780‑5558066                                                                                                                                                 |                                                                                                                                               cTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCa                                                                                                                                              <  2:239845/147‑1 (MQ=255) cTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAgg                                                                                                                                            >  1:23419/1‑149 (MQ=255) cTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCGTCTGAgg                                                                                                                                            >  5:208956/1‑149 (MQ=255)                 ccAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCCGGGCCGCGGCGTCAGATc                                                                                                                            <  3:161236/149‑1 (MQ=255)                   aaCCAGCGTGGTGCGCCCACAACGGAGCGGGTTAACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGTATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAg                                                                                                                          <  3:238690/149‑1 (MQ=38)                          gTTGTGGCCCCACATAGCACCGGTTTAACCCGCGAATACGTTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGaa                                                                                                                   <  1:105008/149‑1 (MQ=255)                                 gCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCTGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAgg                                                                                                                 >  5:74361/1‑144 (MQ=255)                                  cccACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGAt                                                                                                            >  2:209422/1‑148 (MQ=255)                                  cccACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATc                                                                                                           <  7:214241/149‑1 (MQ=255)                                      caTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGc                                                                                                       <  3:57375/149‑1 (MQ=255)                                      caTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGc                                                                                                       >  2:273706/1‑149 (MQ=255)                                      aaTAGGAGGGGTTTTACCCGCTAATACGGTAGCGGCGGAAACGGTGGATTACGGGTGTATATTGGCCCGCACAGACATCCCCACAAATCAAACCGCCGCGGCCCTTCTCCGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGc                                                                                                       <  6:208956/148‑1 (MQ=255)                                       atcggtgCGGTTTTACCCGGTAATACTGTAGCGGCGGCAACGTTGGATTACCGTTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCGGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCg                                                                                                      <  2:257026/143‑1 (MQ=255)                                         aGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCacgaga                                                                                                     >  2:115501/1‑143 (MQ=255)                                          ggAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAAc                                                                                                   <  4:128562/149‑1 (MQ=255)                                            aGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACgg                                                                                                 <  1:209422/149‑1 (MQ=255)                                              ctatcttACCCGCTAATAGGGTTGGGGGGCCAACGGTGTATTACGGGTGTATTTTGGCCAGAACATACATGCCCAAAAATCAAAACGCCGGGGCCCTTCTCAGGGCGGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGt                                                                                               <  6:35789/144‑1 (MQ=255)                                               gggTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCGTATGACGGCCCCCGCGTCAAAAAAACCAGGAAAGGAGAACACCCGAACCGGt                                                                                               >  4:7377/1‑148 (MQ=21)                                                ggTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGa                                                                                             >  4:233906/1‑149 (MQ=255)                                                 gTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGAc                                                                                            >  2:250488/1‑149 (MQ=255)                                                     aCCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCGTAAGAGGGCCCCCGCGTCAGATAAGCCTGGAAGGCAGATCATGCGAACGGGTGagatag                                                                                        >  5:35789/1‑144 (MQ=38)                                                      cccGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCTGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACgcag                                                                                        <  6:74361/148‑1 (MQ=255)                                                         gcgAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCAc                                                                                     >  2:172689/1‑148 (MQ=255)                                                         gcgAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCCGCCCGTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAg                                                                                                                  >  1:305818/1‑119 (MQ=255)                                                              tACGGTAGCGGCTGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCTCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCATCACGATGGt                                                                               <  6:127328/149‑1 (MQ=255)                                                                 ggTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGtctc                                                                            >  3:242600/1‑149 (MQ=255)                                                                       ggcggcAACGGTGGATTACCGGCGGATATTGGCCTGCACAGACATGCCCACAATTGAAAACGGCGCTGCGCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCg                                                                      <  5:290065/149‑1 (MQ=255)                                                                             aaCGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGa                                                                                             <  7:177884/120‑1 (MQ=255)                                                                             aaCGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCgg                                                                >  1:12963/1‑149 (MQ=255)                                                                                ggTGGATTACCGGTGGATATTGGCCAGCACAGCCATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCCGGcc                                                             >  7:222256/1‑149 (MQ=255)                                                                                  tGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCgg                                                           >  3:312391/1‑149 (MQ=255)                                                                                     ataccgggTGTATATGGCCCACAACACACACCCCCACAAAAAAACCCGCCGCGGCCCTTCTCTGGCCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCg                                                        <  3:7377/143‑1 (MQ=255)                                                                                          cGGTGTATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAatg                                                   <  4:242600/149‑1 (MQ=255)                                                                                          cGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAat                                                    >  4:101851/1‑148 (MQ=255)                                                                                          cGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAat                                                    >  1:82247/1‑148 (MQ=255)                                                                                          cGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAat                                                    >  1:86590/1‑148 (MQ=255)                                                                                                      gCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgc                                       <  2:111402/149‑1 (MQ=255)                                                                                                      gCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTAAAAGGGCCGCGGCGTCAGATCAAACAGGAAGGCAGAACAGGCGAACCGGTTACGCACAACGATGCTCTCGTTGCGGTACGGGCCGGTCGAAATGATGTCGATCCgc                                       >  6:1741/1‑149 (MQ=255)                                                                                                       ccAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgcg                                      >  3:284212/1‑149 (MQ=255)                                                                                                          gCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgcgcgc                                   <  1:172689/149‑1 (MQ=255)                                                                                                           cacaGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACgcagc                                                                                       >  7:9287/1‑96 (MQ=255)                                                                                                           cacaGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCg                                                                 <  8:184162/118‑1 (MQ=255)                                                                                                                       cacaAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTc                      <  2:82247/149‑1 (MQ=255)                                                                                                                       cacaAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTc                      <  2:86590/149‑1 (MQ=255)                                                                                                                       cacaAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTc                      >  4:14193/1‑149 (MQ=255)                                                                                                                        acaAATGAAACCGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCg                     <  1:273706/149‑1 (MQ=255)                                                                                                                           aaTGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGACTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCTCGATGGGCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATg                  >  8:230314/1‑149 (MQ=255)                                                                                                                                tgaCGCCGCGGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCtt             >  2:86166/3‑149 (MQ=255)                                                                                                                                ccccgccggGGCCCTTCTCTGGCCGGCGGCGTCAGATCAGCCGGGACGGCAGAGGAGGCGCACGGGTGACGCTGCACGATGGTCTCGTGGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCt              <  6:178108/145‑1 (MQ=255)                                                                                                                                aaaCGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGATCGACATCATTTCGACCgg                                        <  6:283047/122‑1 (MQ=255)                                                                                                                                aaaCGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGATCGACATCATTTCGACCgg                                        >  5:283047/1‑122 (MQ=255)                                                                                                                                  aCGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATAAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCCGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCCACCAGCACTTAGATGCGCTTGa           >  4:89806/1‑149 (MQ=255)                                                                                                                                  aCGCCGCGGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGa           <  1:86166/149‑1 (MQ=255)                                                                                                                                  aCGCCGCGGCCCTGAGAAGGGCCGCGGCGTCAGATAAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGa           >  6:16743/1‑149 (MQ=255)                                                                                                                                    gccgcGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATg         >  3:35730/1‑149 (MQ=255)                                                                                                                                    gccgcGGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCg                                                            >  5:94827/1‑98 (MQ=38)                                                                                                                                    gccgcGGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCg                                                            <  6:94827/98‑1 (MQ=38)                                                                                                                                    gccgcGGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGCGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATg         >  8:183613/1‑149 (MQ=255)                                                                                                                                      cgcgGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATg         >  1:221080/1‑147 (MQ=255)                                                                                                                                      cgcgGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGaaa                                                     >  7:112655/1‑103 (MQ=39)                                                                                                                                         ggCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTc    <  3:89806/149‑1 (MQ=255)                                                                                                                                           cccTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCGCGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAc  <  3:233906/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                               CTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAC  >  NC_002947/5557780‑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: ‑