Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,539,112	G→T	7.0%	T625K (ACG→AAG) ‡	PP_4872 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,539,112	0	G	T	7.0%	224.4 / 5.4	71	T625K (ACG→AAG) ‡	PP_4872	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/32);  new base T (3/2);  total (37/34)
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.91e-01

AGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTT  >  NC_002947/5538971‑5539260                                                                                                                                              |                                                                                                                                                     agcaTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggccgtggcc                                                                                                                                                 >  3:35432/1‑147 (MQ=255) agcaTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGa                                                                                                                                               <  4:178241/149‑1 (MQ=255) agcaTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGa                                                                                                                                               <  4:182145/149‑1 (MQ=255)     tCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGc                                                                                                                                            >  4:254110/1‑148 (MQ=255)        gCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGACCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCAtt                                                                                                                                        >  4:130776/1‑149 (MQ=255)         ccAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTg                                                                                                                                       <  3:254110/149‑1 (MQ=255)         ccAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTg                                                                                                                                       <  2:180874/149‑1 (MQ=255)           aGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCTCGGCCACGGCCAtg                                                                                                                                        <  8:228304/146‑2 (MQ=255)           aGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCTCGGCCACGGCCAtg                                                                                                                                        >  7:228304/1‑145 (MQ=255)             cccTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGa                                                                                                                                   <  2:271282/149‑1 (MQ=255)                  cGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGACCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGActg                                                                                                                                <  3:130776/147‑3 (MQ=255)                        gTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACg                                                                                                                        >  1:290550/1‑149 (MQ=255)                           gTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCt                                                                                                                     >  6:112711/1‑149 (MQ=255)                                gcAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGAGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTAcc                                                                                                                >  7:167675/1‑149 (MQ=255)                                 cAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAg                                                                                                                                              >  1:68533/1‑118 (MQ=255)                                 cAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAg                                                                                                                                              <  2:68533/118‑1 (MQ=255)                                  aaGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTg                                                                                                              >  2:198707/1‑149 (MQ=255)                                    gAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAgg                                                                                                                                             <  4:74512/116‑1 (MQ=255)                                    gAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAgg                                                                                                                                             >  3:74512/1‑116 (MQ=255)                                    gAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGt                                                                                                            <  1:285533/149‑1 (MQ=255)                                          cTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGgcgc                                                                                                      >  7:293819/1‑149 (MQ=255)                                           ttGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCTTGGCCGAGGCCATTgg                                                                                                                                      >  7:234434/1‑116 (MQ=255)                                           ttGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTgg                                                                                                                                      <  8:87103/116‑1 (MQ=255)                                           ttGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTgg                                                                                                                                      <  8:234434/116‑1 (MQ=255)                                           ttGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGcgcg                                                                                                     >  2:103068/1‑149 (MQ=255)                                           ttGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGcgcg                                                                                                     >  2:305075/1‑149 (MQ=255)                                           ttGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACACCTTCGCTCTTGTAGGTCACCTTGAGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTgg                                                                                                                                      >  7:87103/1‑116 (MQ=255)                                              tgCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGaa                                                                                                  <  5:282270/149‑1 (MQ=255)                                                cAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGc                                                                                                <  2:290550/149‑1 (MQ=255)                                                     ggCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTg                                                                                           <  8:293819/149‑1 (MQ=255)                                                      gCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGc                                                                                          <  3:77431/149‑1 (MQ=255)                                                                    tctTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAgg                                                                            >  8:154428/1‑149 (MQ=255)                                                                    tctTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGAGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAgg                                                                            >  3:73642/1‑149 (MQ=255)                                                                        gatgATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGc                                                                          >  1:51382/4‑147 (MQ=255)                                                                        gACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGcc                                                                         <  5:112711/148‑1 (MQ=255)                                                                         atgATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCa                                                                        <  2:51382/146‑1 (MQ=255)                                                                          cGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAgg                                                                      <  7:154428/149‑1 (MQ=255)                                                                           gATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGccc                                                                                        <  6:204337/130‑1 (MQ=255)                                                                           gATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGccc                                                                                        >  5:204337/1‑130 (MQ=255)                                                                             tctcCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGa                                                                                                                                               <  6:238836/73‑1 (MQ=255)                                                                             tctcCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGa                                                                                                                                               >  5:238836/1‑73 (MQ=255)                                                                               tcCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTtcat                                                                 >  1:316995/1‑149 (MQ=255)                                                                               tcCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTtcat                                                                 >  6:321169/1‑149 (MQ=255)                                                                                ccACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAAc                                                                                                                         <  7:242321/92‑1 (MQ=255)                                                                                ccACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAAc                                                                                                                         >  8:242321/1‑92 (MQ=255)                                                                                 cACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAgg                                                                      <  8:167675/142‑1 (MQ=255)                                                                                 cACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTtcatca                                                               <  5:148985/149‑1 (MQ=255)                                                                                 cACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTtcatca                                                               >  6:148985/1‑149 (MQ=255)                                                                                 cACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTtcatc                                                                <  1:198707/148‑1 (MQ=255)                                                                                  aCCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTcatcat                                                              <  5:321169/149‑1 (MQ=255)                                                                                    cccTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGt                                                            >  8:10618/1‑149 (MQ=255)                                                                                       ttCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCtt                                                         >  2:92023/1‑149 (MQ=255)                                                                                       ttCGCTCTTGTAGGTCACCTTGCGCAGCGGGTCCTTGCCTGAGTGGTACATGGCCGTGGCCTAGGCCATTGGCGACGGGTAGACCGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCtt                                                         >  6:71262/1‑149 (MQ=255)                                                                                                gTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCggtgg                                                >  3:271188/1‑149 (MQ=255)                                                                                                gTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCggt                                                  <  5:286022/147‑1 (MQ=255)                                                                                                    gTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTAcc                                            >  6:13316/1‑149 (MQ=255)                                                                                                    gTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTAcc                                            >  7:15510/1‑149 (MQ=255)                                                                                                    gTCACCTTGCGCAGCGGGTTCTTGCCTGACTGCAACATGGCAGTGGCCGAGGCCAAAGGCGACGGGTAGAACGCCTGTACCTGGTCGCCGCGGAAGCCGTTGCCCATCAGCCACAGGGCCAGGTACATCATGTCATCGTCGGTGGTAcc                                            >  3:232366/1‑149 (MQ=255)                                                                                                    gTCACCTTGAGCACCCGGTACTTTCCTGAGTGGAAAATGGCCCTGGCCGAGGCCATTGGCGACGCGTAGAACGCCTGTACCTCGTCGGCGCGGAAGCCCTTCCCCTTCAGCCACAGGGCCAGGATCATCATTTCTTCGTCTGTGGTAca                                            >  8:278048/1‑148 (MQ=255)                                                                                                     tCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCg                                           >  6:71716/1‑149 (MQ=255)                                                                                                              gcAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCTACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGcgcg                                  >  3:164311/1‑149 (MQ=255)                                                                                                                           gCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTa                       <  4:35432/147‑1 (MQ=255)                                                                                                                           gCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGACAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACg                     >  6:323824/1‑149 (MQ=255)                                                                                                                            ccTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACgg                    <  8:15510/149‑1 (MQ=255)                                                                                                                            ccTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACg                     <  1:92023/148‑1 (MQ=255)                                                                                                                            ccTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACg                     <  5:71262/148‑1 (MQ=255)                                                                                                                              tGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGTGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCGTCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGaa                          <  1:26675/141‑1 (MQ=255)                                                                                                                              tGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCGTCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGaa                          >  2:26675/1‑141 (MQ=255)                                                                                                                               gcgTGGGAAATGTCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGaa                          <  6:221420/138‑1 (MQ=255)                                                                                                                               gAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGaa                          >  5:221420/1‑140 (MQ=255)                                                                                                                                     ccaaTGGCCTCGGCCACGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTtcgt                                                      <  2:228428/103‑1 (MQ=255)                                                                                                                                     ccaaTGGCCTCGGCCACGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTtcgt                                                      >  1:228428/4‑106 (MQ=255)                                                                                                                                         tggccgtggccgAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGTGGCAATGGAGTACGGGATCAGGTACTGc       <  2:316995/149‑1 (MQ=255)                                                                                                                                           gccgtggccgAGGCCTTTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCTCAGGGCCAGGTTGATGATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCt      <  4:232366/148‑1 (MQ=255)                                                                                                                                              gtggccgAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACTGGATCAGGTACTGCTCCtt  <  4:73642/149‑1 (MQ=255)                                                                                                                                              gtggccgAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCtt  <  4:271188/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                     AGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTT  >  NC_002947/5538971‑5539260

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