Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,567,948	T→A	9.9%	intergenic (‑21/‑82)	PP_2250 ← / → PP_2251	LysR family transcriptional regulator/membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,567,948	0	T	A	9.9%	161.7 / 5.7	60	intergenic (‑21/‑82)	PP_2250/PP_2251	LysR family transcriptional regulator/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (23/31);  new base A (2/4);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAAGGTTGCCAATGGCCATGCTCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCAAGTTGGTCTACGTGACGTTCTATGAACTGATCGGGCTGTGCATGT  >  NC_002947/2567802‑2568093                                                                                                                                                    |                                                                                                                                                  ttCAAGGTTGCCAATGGCCATGCTCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGATAAGTTCATGGCAGCCTCCTGTCAGTTGATCt                                                                                                                                                 >  2:34368/1‑149 (MQ=255) ttCAAGGTTGCCAATGGCCATGCTCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCt                                                                                                                                                 >  4:204260/1‑149 (MQ=255) ttCAAGGTTGCCAATGGCCATGCTCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCt                                                                                                                                                 >  4:208305/1‑149 (MQ=255)   cAAGGTTGCCAATGGCCATGCTCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCt                                                                                                                                                 <  6:268466/147‑1 (MQ=255)    aaGGTTGCCAATGGCCATGCTCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGa                                                                                                                                               <  3:208305/148‑1 (MQ=255)    aaGGTTGCCAATGGCCATGCTCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGa                                                                                                                                               <  3:204260/148‑1 (MQ=255)          gCCACTTGCCATGCTCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGATCAACTGACAggag                                                                                                                                         <  7:269616/148‑4 (MQ=255)           ccAATGGCCATGCTCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGATCAACTGACAggagg                                                                                                                                        <  8:102411/148‑5 (MQ=255)                  ccATGCTCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACtttt                                                                                                                                >  6:47776/1‑149 (MQ=255)                  ccATGCTCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACtttt                                                                                                                                >  6:43071/1‑149 (MQ=255)                   cATGCTCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTg                                                                                                                               <  8:5337/149‑1 (MQ=255)                     tGCTCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCt                                                                                                                             <  8:239740/149‑1 (MQ=255)                       cTCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGt                                                                                                                           >  2:32706/1‑149 (MQ=255)                        tCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGt                                                                                                                           >  1:66215/1‑148 (MQ=255)                        tCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGATCAACTGAc                                                                                                                                              >  5:289181/1‑129 (MQ=255)                        tCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGATCAACTGAc                                                                                                                                              <  6:289181/129‑1 (MQ=255)                          aCGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTc                                                                                                                          <  2:83321/147‑1 (MQ=255)                          aCGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGATCAACTGACAggagg                                                                                                                                        <  5:81148/133‑5 (MQ=255)                          aCGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGATCAACTGACAggagg                                                                                                                                        >  6:81148/1‑129 (MQ=255)                          aCGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGAGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGAAGTTGTCGCTGGAAAAGTTCATGACAGCCTCCTGTCAGTTTCTCTGACAGCTGCTGACTTATGCTGTCaa                                                                                                                        >  6:54425/1‑149 (MQ=255)                           cGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAg                                                                                                                       <  5:187586/149‑1 (MQ=255)                               tGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCg                                                                                                                   >  5:31091/1‑149 (MQ=255)                                 aaGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGtt                                                                                                                 >  6:270237/1‑149 (MQ=255)                                 aaGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCCGACTTTTGCTGTCAAGTGCGtt                                                                                                                 <  4:102419/149‑1 (MQ=255)                                     cACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTg                                                                                                                >  6:115062/1‑146 (MQ=255)                                     cACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGaa                                                                                                              <  5:115062/148‑1 (MQ=255)                                     cACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAg                                                                                                             >  6:228509/1‑149 (MQ=255)                                         cgcTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATg                                                                                                         <  2:258111/149‑1 (MQ=255)                                               aGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGcc                                                                                                   <  2:66215/149‑1 (MQ=255)                                                  gcgcgcgcTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGcccc                                                                                                 >  1:123365/1‑148 (MQ=255)                                                    gcgcgcTGCGGCGGAAAACGAACCGTGGTGCGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCa                                                                                               >  1:179501/1‑149 (MQ=255)                                                        gcTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCt                                                                                          <  6:31091/149‑1 (MQ=255)                                                         cTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCtt                                                                                         <  5:228509/149‑1 (MQ=255)                                                                       gAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGATCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTgg                                                                              <  1:224710/146‑1 (MQ=255)                                                                       gAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGATCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTgg                                                                              >  2:224710/1‑146 (MQ=255)                                                                                           ccAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGAt                                                        <  6:201736/149‑1 (MQ=255)                                                                                                     cTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCaa                                              >  8:55435/1‑149 (MQ=255)                                                                                                     cTGGATGTGGTCGCCGGAAAAGTTCATGGCAGCCTCCTGGCAGTGGATCTGAGAGCGGCGGACTTTTGCGGTCAAGTGCGTTGAAGCATCATCGCCCCCATCGCTTTTTGTCGTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCaa                                              <  8:283060/149‑1 (MQ=255)                                                                                                      tGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCG‑CCCCATCGCTTTTTGTCTTTgg                                                                              <  3:60168/115‑1 (MQ=255)                                                                                                      tGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCG‑CCCCATCGCTTTTTGTCTTTgg                                                                              >  4:60168/1‑115 (MQ=255)                                                                                                        gATGTTTTCGCTGGATAAGTTCATTGGAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCTCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCTAGt                                            <  3:137578/148‑1 (MQ=255)                                                                                                            ttgtCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCAAGTTGGTc                                       >  4:162928/1‑149 (MQ=255)                                                                                                               tCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCAAGTTGGTCTAc                                    <  5:43071/149‑1 (MQ=255)                                                                                                               tCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCAAGTTGGTCTAc                                    <  5:47776/149‑1 (MQ=255)                                                                                                                     gAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCAAGTTGGTCTACGTGa                                <  5:270237/147‑1 (MQ=255)                                                                                                                        aaGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCAAGTTGGTCTACGTGACGTTc                           <  2:179501/149‑1 (MQ=255)                                                                                                                        aaGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTACTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCAAGTTGGTCTACGTGACGTTc                           >  4:7541/1‑149 (MQ=255)                                                                                                                              aTGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCAAGTTGGTCTACGTGACGTTCTATGaa                     <  1:32706/149‑1 (MQ=255)                                                                                                                              aTGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCAAGTTGGTCTACGTGACGTTCTATGaa                     <  1:34368/149‑1 (MQ=255)                                                                                                                                ggCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCAAGTTGGTCTACGTGACGTTCTATGAACt                   >  1:180837/1‑149 (MQ=255)                                                                                                                                     cctcctGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCAAGTTGGTCTACGTGACGTTCTATGAACTGATCg              <  2:180837/149‑1 (MQ=255)                                                                                                                                        cctGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGaaa                                                   <  1:90971/109‑1 (MQ=255)                                                                                                                                        cctGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGaaa                                                   >  2:90971/1‑109 (MQ=255)                                                                                                                                        cctGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCAAGTTGGTCTACGTGACGTTCTATGAACTGATCGGGc           >  8:255796/1‑149 (MQ=255)                                                                                                                                         ctGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCgg                                                                  <  4:73812/93‑1 (MQ=255)                                                                                                                                         ctGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCgg                                                                  >  3:73812/1‑93 (MQ=255)                                                                                                                                         ctGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCAAGTTGGTCTACGTGACGTTCTATGAACTGATCggg            <  7:55435/147‑1 (MQ=255)                                                                                                                                            tCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCAAGTTGGTCTACGTGACGTTCTATGAACTGATCGGGCtgtg       <  3:162928/149‑1 (MQ=255)                                                                                                                                            acAGTTGATCTGACAGCTGCTGACTGTTGCTGTCAAGTTCGTGCAAGCATCATCGCCCCCATCGCATTTTGTCTTTGGATGTAACCCCGGTTGCAGGGATAGAAACGCAAGTTGGTCTACGTGACGTTCTATGAACTGATCGGGCtgt        <  8:108182/147‑1 (MQ=255)                                                                                                                                                  gATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCGATCGCTTTTTGTCTTTGGAGGTAACGCCGGGTGCAGGGATTGAAACGCAAGTTGGTCTACGTGACGTTCTATGAACTGATCGGGCTGTGCATGt  <  5:54425/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TTCAAGGTTGCCAATGGCCATGCTCACGGCTGAAGGCACGCGCTTGAGGGCGCGCGCTGCGGCGGAAAACGAACCGTGGT‑CGAGCACGGCCAGAAACAGCTGGATGTTGTCGCTGGAAAAGTTCATGGCAGCCTCCTGTCAGTTGATCTGACAGCTGCTGACTTTTGCTGTCAAGTGCGTTGAAGCATGATCGCCCCCATCGCTTTTTGTCTTTGGAGGTAACCCCGGGTGCAGGGATTGAAACGCAAGTTGGTCTACGTGACGTTCTATGAACTGATCGGGCTGTGCATGT  >  NC_002947/2567802‑2568093

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