Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	162,281	A→T	5.9%	intergenic (‑315/+59)	PP_0153 ← / ← scpC	hypothetical protein/propionyl‑CoA:succinate CoA transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	162,281	0	A	T	5.9%	203.9 / 2.8	68	intergenic (‑315/+59)	PP_0153/scpC	hypothetical protein/propionyl‑CoA:succinate CoA transferase
Reads supporting (aligned to +/- strand):  ref base A (35/29);  new base T (2/2);  total (37/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 = 5.51e-01

ACTTATTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGGCC  >  NC_002947/162136‑162420                                                                                                                                                   |                                                                                                                                            aCTTATTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGc                                                                                                                                          >  4:1027/1‑149 (MQ=255)  cTTATTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGcc                                                                                                                                         >  1:172904/1‑149 (MQ=255)   ttATTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGccc                                                                                                                                        >  1:20957/1‑149 (MQ=255)   ttATTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGccc                                                                                                                                        >  3:143840/1‑149 (MQ=255)      tttGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAGAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCcgcg                                                                                                                                     >  1:131425/1‑149 (MQ=255)          gTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCTTGCCCGCGaat                                                                                                                                  <  7:184471/148‑2 (MQ=255)             cGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCTCTTCGCGGGCTTGCCCGCGaatgggc                                                                                                                              <  5:315087/149‑6 (MQ=255)              gTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGcc                                                                                                                             >  4:255176/1‑149 (MQ=255)              gTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGcc                                                                                                                             <  3:90145/149‑1 (MQ=255)                      tCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGc                                                                                                                     >  7:107467/1‑149 (MQ=255)                      tCGAAAAAGGGGATTTCTCAGACGGGAAGGACTGAAAGAACAAGCACAGGAGTACAGATTATGCTGCTAGAAGAGGGGTGGCGAAGTGATAC‑AAGGGGTTGACTGTTATGGCCCATTCGCGGGCAAGACCGGTAAGAGCCCGGTACGGc                                                                                                                     >  6:364105/1‑149 (MQ=255)                       cGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGc                                                                                                                     <  3:1027/148‑1 (MQ=255)                          agaagGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGaa                                                                                                                 <  3:200080/149‑1 (MQ=255)                              ggggATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAg                                                                                                             <  5:297912/149‑1 (MQ=255)                               gggATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAgg                                                                                                            <  7:334265/149‑1 (MQ=255)                               gggATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAgg                                                                                                            >  7:315812/1‑149 (MQ=255)                                     tGTCAGTCGGGAAGGTCTGAGAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGgccgc                                                                                                       <  2:131425/148‑1 (MQ=255)                                         aGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGtt                                                                                                   <  8:315812/148‑1 (MQ=255)                                         aGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTc                                                                                                  <  2:193809/149‑1 (MQ=255)                                          gTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTcc                                                                                                 <  4:31909/149‑1 (MQ=255)                                           tCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCg                                                                                                >  3:307868/1‑149 (MQ=255)                                                        aaaGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAgg                                                                                                            >  3:316171/1‑124 (MQ=255)                                                        aaaGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAgg                                                                                                            <  4:316171/124‑1 (MQ=255)                                                        aaaGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGa                                                                                   >  3:231203/1‑149 (MQ=255)                                                         aaGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGAt                                                                                  <  2:20957/149‑1 (MQ=255)                                                              ccAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAg                                                                                                             >  7:250416/1‑117 (MQ=255)                                                              ccAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAg                                                                                                             <  8:250416/117‑1 (MQ=255)                                                                aGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGc                                                                                                                     >  4:319950/1‑107 (MQ=255)                                                                aGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGc                                                                                                                     <  3:319950/107‑1 (MQ=255)                                                                 gCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGc                                                                          >  1:167836/1‑149 (MQ=255)                                                                          taCAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGaa                                                                                                                         <  1:280401/93‑1 (MQ=255)                                                                          taCAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGaa                                                                                                                         >  2:280401/1‑93 (MQ=255)                                                                          taCAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCAGATCCGGATCa                                                                                >  6:25492/1‑134 (MQ=255)                                                                          taCAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCAGATCCGGATCa                                                                                <  5:25492/134‑1 (MQ=255)                                                                          taCAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATgcg                                                                 >  1:22667/1‑149 (MQ=255)                                                                          taCAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATgcg                                                                 >  2:213472/1‑149 (MQ=255)                                                                                ttATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGt                                                           <  2:22667/149‑1 (MQ=255)                                                                                 tATGCTGCTAGTAGAGGGGTGGGGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGt                                                           <  4:190408/148‑1 (MQ=255)                                                                                     ctgctAGTAGAGGGGTGGCGTAGTGATAGAAAGGGTTTGACTGTTATGGCCCATACGCGGGAAACCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTttc                                                        >  1:28033/1‑148 (MQ=255)                                                                                     ctgctAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTttctt                                                      <  8:107467/149‑1 (MQ=255)                                                                                     ctgctAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTttctt                                                      >  4:161341/1‑149 (MQ=255)                                                                                      tgctAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTttcttc                                                     >  7:24538/1‑149 (MQ=255)                                                                                           gtagAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAgg                                                 <  4:307868/148‑1 (MQ=255)                                                                                                 gggTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGt                                          >  6:264585/1‑149 (MQ=255)                                                                                                      gCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAg                                     <  4:231203/149‑1 (MQ=255)                                                                                                      gCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAg                                     >  6:148314/1‑149 (MQ=255)                                                                                                           gTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCa                                 <  8:24538/148‑1 (MQ=255)                                                                                                           gTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAg                                >  8:264485/1‑149 (MQ=255)                                                                                                                    aGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTc                          <  8:4745/147‑1 (MQ=255)                                                                                                                    aGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTc                          >  7:4745/1‑147 (MQ=255)                                                                                                                      ggTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCgcgc                      <  5:264585/149‑1 (MQ=255)                                                                                                                      ggTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCgcgc                      <  3:255176/149‑1 (MQ=255)                                                                                                                        tttGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAg                    >  5:344600/1‑149 (MQ=255)                                                                                                                        tttGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAg                    >  1:220141/1‑149 (MQ=255)                                                                                                                          tGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAgg                   <  1:213472/148‑1 (MQ=255)                                                                                                                            aCTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGAtg                >  3:171327/1‑149 (MQ=255)                                                                                                                              tGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGAtgtg              >  2:166062/1‑149 (MQ=255)                                                                                                                               gTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGc             >  5:135589/1‑149 (MQ=255)                                                                                                                                ttATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCg            >  8:162033/1‑149 (MQ=255)                                                                                                                                   tGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCggg          <  5:148314/148‑1 (MQ=255)                                                                                                                                   tGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCggg          <  4:143840/148‑1 (MQ=255)                                                                                                                                    ggCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGt        >  8:315142/1‑149 (MQ=255)                                                                                                                                      cccATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCgggg         <  6:8444/146‑1 (MQ=255)                                                                                                                                      cccATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGtt       >  5:8444/1‑147 (MQ=255)                                                                                                                                       ctcttCGCGGGCTTGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATg     >  4:318619/4‑149 (MQ=255)                                                                                                                                       ccATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATg     >  5:377099/1‑149 (MQ=255)                                                                                                                                       ccATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATg     <  1:166062/149‑1 (MQ=255)                                                                                                                                          ttCGCGGGCTTGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGgcc  >  8:123690/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                            ACTTATTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAG‑AAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGGCC  >  NC_002947/162136‑162420

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