Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,256,865	T→A	8.6%	intergenic (+17/+28)	PP_5703 → / ← PP_4634	type III effector protein HopJ/alpha/beta family hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,256,865	0	T	A	8.6%	153.8 / 5.4	58	intergenic (+17/+28)	PP_5703/PP_4634	type III effector protein HopJ/alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base T (28/25);  new base A (3/2);  total (31/27)
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.77e-01

CCTTCGGCGAGCACTACCGTGATGTGTTGGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAA  >  NC_002947/5256718‑5257008                                                                                                                                                    |                                                                                                                                                 ccTTCGGCGAGCACTACCGTGATGTGTTGGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCttt                                                                                                                                                 <  3:338790/149‑1 (MQ=255)  gttCGGCGAGCACTAGCGTGATGTGTTGGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGGACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTg                                                                                                                                                <  4:235611/148‑1 (MQ=255)       gCGAGCACTACCGTGATGTGTTGGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGcc                                                                                                                                           >  7:89459/1‑149 (MQ=255)       gCGAGCACTACCGTGATGTGTTGGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCAAAGCGGcc                                                                                                                                           >  5:66705/1‑149 (MQ=255)           gCACTACCGTGATGTGTTGGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCAAAGCGGcc‑ccgg                                                                                                                                      <  8:256816/149‑3 (MQ=255)            cACTACCGTGATGTGTTCGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGCCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTCGAGATTGCGGGGGCCGCTTTGCGGCC‑CCACg                                                                                                                                     <  6:236344/149‑1 (MQ=255)                 ccGTGATGTGTTGGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCAAAGCGGcc‑ccggcaat                                                                                                                                  <  6:66705/147‑7 (MQ=255)                           ttGGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGc                                                                                                                                               >  5:188192/1‑125 (MQ=255)                           ttGGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGc                                                                                                                                               <  6:188192/125‑1 (MQ=255)                              gCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGCTCCTCCGCCCACCCTc                                                                                                                    >  1:284452/1‑148 (MQ=255)                              gCCACGCCTGAGGGTAGCGACCATGGCAACATCAGTGCACTGAGCAAACACGGTCTGGCGGGCGTGAAGTATGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTTGGGGC‑GCACGTTGCTGAGccccggcga                                                                                                                    >  8:283964/1‑142 (MQ=255)                               ccACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCa                                                                                                                   >  1:341250/1‑148 (MQ=255)                                cACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACCCTCCAGCGCCCCCCCTCaa                                                                                                                  >  1:359125/1‑148 (MQ=255)                                   gCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCa                                                                                                              <  5:85999/149‑1 (MQ=255)                                    ccTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCat                                                                                                             <  7:276540/149‑1 (MQ=255)                                          ggTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGt                                                                                                       <  2:341250/149‑1 (MQ=255)                                              gCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCa                                                                                                   <  1:159377/149‑1 (MQ=255)                                               cGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATcc                                                                                                               <  8:107311/136‑1 (MQ=255)                                               cGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATcc                                                                                                               >  7:107311/1‑136 (MQ=255)                                                   cATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGc                                                                                                                                               >  4:336683/1‑101 (MQ=255)                                                   cATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGc                                                                                                                                               <  3:336683/101‑1 (MQ=255)                                                      ggCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGa                                                                                                                                    <  4:336930/108‑1 (MQ=255)                                                      ggCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGa                                                                                                                                    >  3:336930/1‑108 (MQ=255)                                                        cAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGc                                                                                         >  7:56532/1‑149 (MQ=255)                                                         aaCATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGc                                                                                                                                               >  8:6409/1‑95 (MQ=255)                                                         aaCATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGc                                                                                                                                               <  7:6409/95‑1 (MQ=255)                                                              ccGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCTCTTTGCGGCC‑ACACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTc                                                                                     >  1:238625/1‑147 (MQ=255)                                                               cGTGCACTCATCACACACGGGCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCa                                                                                   <  2:286814/148‑1 (MQ=255)                                                                gTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGc                                                                                 <  8:56532/149‑1 (MQ=255)                                                                    aCTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTcc                                                                             >  6:153046/1‑149 (MQ=255)                                                                      tGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTcccc                                                                           >  3:296156/1‑149 (MQ=255)                                                                        aTCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTg                                                                         >  2:288544/1‑149 (MQ=255)                                                                         tCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGt                                                                                                       >  1:155436/1‑118 (MQ=255)                                                                           aaaCACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCg                                                                      >  6:109186/1‑149 (MQ=255)                                                                           aaaCACGGTCTGGCGCGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGTGGCCGCGTTGAGGAC‑CCACGATCCTACGCCCCCCCTCAATCAATGATGGTGCCACGCCCCCCGCCTGACAGACTCCCCGGCCg                                                                      >  4:250967/1‑149 (MQ=255)                                                                            aaCACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCg                                                                      >  5:112112/1‑148 (MQ=255)                                                                             acacGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATACAGCGCCCCACCGCAATCCAAGATTGTGCCAGGCCCCCCGCCATCCAGCCAACCCTGCCGc                                                                     >  8:150541/1‑148 (MQ=255)                                                                               acGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCCCAAGGCTTGAGATTGCCGGGGCCGCATTGCGGCC‑CCCCGAGCCAGCCCCCACCCACAAACCAAGTACGTGGCAGCCCCCCCCCCATGCACCCGACCCAGCCGcccc                                                                  >  4:79733/1‑146 (MQ=255)                                                                                     tGGCGGGAGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGTCCCCTTAGCGGCC‑CAACGATACTGCGCCCAACCTAAATACAAGCTCTTGCCAGGCCCCCCCCCTTCCACCCTCCCATCACGCAACCCCtc                                                             >  4:313559/1‑147 (MQ=255)                                                                                         gggcgTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAgg                                                        <  1:133701/149‑1 (MQ=255)                                                                                         gggcgTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAgg                                                        <  6:48060/149‑1 (MQ=255)                                                                                                 aGTTTGCTGCGCAGCGGCTTGGGTGCAAGGCTTCATCTGGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGc                                                <  2:284452/149‑1 (MQ=255)                                                                                                       ctgcGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAg                                           <  1:288544/148‑1 (MQ=255)                                                                                                       ctgcGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAg                                           <  8:15688/148‑1 (MQ=255)                                                                                                        tgcGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCCCCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGc                                         >  7:69003/1‑149 (MQ=255)                                                                                                              gCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTcgc                                   <  7:150541/149‑1 (MQ=255)                                                                                                                cGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCGTTGGGGCC‑CCACGATCCTGCGCCCCCCCTCCATCCATGATGGTGCCAGGCCCCCCGCCGTGCCGCCGCCCCGGCCGCCACCCCGGCAAGGCCCCCGGCCACCGGCGTTggcg                                  >  4:146473/1‑148 (MQ=255)                                                                                                                     gcgcgcAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCCGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCg                            >  3:51697/1‑149 (MQ=255)                                                                                                                          cAAGGCTTGAGATTGCCGTGGCCGCTTTGCGGCC‑CCACGATCCTGCTCCCACCCTCAATCCATTATGGTGACAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGTCCACAGGCTTTCGCGAAACCGGTCAt                       <  2:238625/149‑1 (MQ=255)                                                                                                                              gCTTGAGATTGCCGGGGCCGCTTTGCGGCCGCCACGATCATGCGCCCCCCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTGCCAGCCGCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACCGGCTTTCGGGAAACCGGTCATgg                     >  1:35469/1‑148 (MQ=255)                                                                                                                              gCTTGAGATTGCCGGGCCCGCGTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGaa                   <  7:283964/149‑1 (MQ=255)                                                                                                                                  cagaTTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAAc                                                                  <  8:337095/97‑1 (MQ=255)                                                                                                                                       tcgtggggCCGCAAAGCGGCC‑CCGCGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCACGGCCACCGGCCACAGGCGTTCGCGAAACCGGTCAt                       >  5:264054/5‑136 (MQ=255)                                                                                                                                        cgtggggCCGCAAAGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGaaa         >  5:389176/4‑149 (MQ=255)                                                                                                                                         ccGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGt                                                                                                       >  3:117041/1‑54 (MQ=255)                                                                                                                                         ccGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGt                                                                                                       <  4:117041/54‑1 (MQ=255)                                                                                                                                         ccGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGaaa         <  2:359125/148‑1 (MQ=255)                                                                                                                                              gCCGCTTTGCGTCC‑CCACGATCCTGCGCCCACCCTCAAAACATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGAGAAATGGCCa   <  3:313559/149‑1 (MQ=255)                                                                                                                                               ccGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACAGCAAGCCCCCCGGCCACACGCTTTCGCGAAACCCGGCATGGAAAAACGCCAAATGGGCaa  >  8:111137/1‑149 (MQ=255)                                                                                                                                                   tttGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGa                    >  6:111203/1‑127 (MQ=255)                                                                                                                                                   tttGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGa                    <  5:111203/127‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CCTTCGGCGAGCACTACCGTGATGTGTTGGCCACGCCTGAGGGTAGCGACCATGGCAACATCCGTGCACTGATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCC‑CCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAA  >  NC_002947/5256718‑5257008

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