Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,488,275	T→G	8.3%	L535R (CTG→CGG)	PP_4824 →	two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,488,275	0	T	G	8.3%	176.5 / 7.1	60	L535R (CTG→CGG)	PP_4824	two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base G (3/2);  total (31/29)
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.27e-01

GCCAAGGCCGAGCAACAGAACATCGAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGAGTTGCACAGCCACCATTTCCTTTCCAGCAACAA  >  NC_002947/5488127‑5488419                                                                                                                                                       |                                                                                                                                                 gcCAAGGCCGAGAAACAGAACATCGAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCct                                                                                                                                                  <  1:86021/149‑1 (MQ=255)    aaGGCCGAGCAACAGAACATCGAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCcggct                                                                                                                                               <  4:22228/149‑1 (MQ=255)        ccGAGCAACAGAACATCGAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCtggtg                                                                                                                                           <  8:165935/149‑1 (MQ=255)              aaCAGAACATCGAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCGGCTGGTGGTGgcg                                                                                                                                     >  2:98300/1‑149 (MQ=255)                      aTCGAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATc                                                                                                                              >  2:44202/1‑148 (MQ=255)                         gAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAAcg                                                                                                                          >  4:270377/1‑149 (MQ=255)                                tCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAg                                                                                                                   >  4:71188/1‑149 (MQ=255)                                     ttCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGAt                                                                                                                               <  6:42826/132‑1 (MQ=255)                                     ttCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGAt                                                                                                                               >  5:42826/1‑132 (MQ=255)                                        aCCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACg                                                                                                           >  3:208382/1‑149 (MQ=255)                                         cccAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGc                                                                                                          >  4:302020/1‑149 (MQ=255)                                          ccAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGcc                                                                                                         >  6:178694/1‑149 (MQ=255)                                          ccAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGcc                                                                                                         >  6:259336/1‑149 (MQ=255)                                             gCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGc                                                                                                      <  1:329719/149‑1 (MQ=255)                                             gCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCGGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGc                                                                                                      <  5:59697/149‑1 (MQ=255)                                              ccACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCGGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCg                                                                                                     >  4:154689/1‑149 (MQ=255)                                                 caAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTAt                                                                                                  >  8:106111/1‑149 (MQ=255)                                                        cgcgGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAATACCGACCAGTTGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGt                                                                                            <  1:101806/149‑1 (MQ=255)                                                        cgcgGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTATCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCtggtg                                                                                                                                           <  6:212366/101‑1 (MQ=255)                                                                   aGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAgcgg                                                                                <  3:58055/149‑1 (MQ=255)                                                                         gaCCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGcc                                                                          <  3:71188/149‑1 (MQ=255)                                                                         gaCCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGcc                                                                          <  3:270377/149‑1 (MQ=255)                                                                          aCCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTGAAAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGcc                                                                          >  1:84843/1‑149 (MQ=255)                                                                            ccGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGtt                                                                       >  5:154032/1‑149 (MQ=255)                                                                            ccGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGtt                                                                       >  5:120544/1‑149 (MQ=255)                                                                            ccGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGtt                                                                       >  4:80204/1‑149 (MQ=255)                                                                            ccGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGtt                                                                       >  4:75566/1‑149 (MQ=255)                                                                                   gTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCccgcc                                                                >  3:28754/1‑149 (MQ=255)                                                                                         gcCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGAc                                                          >  5:255532/1‑149 (MQ=255)                                                                                           cAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCg                                                        >  6:203152/1‑149 (MQ=255)                                                                                               ccTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGaa                                                    <  5:259336/149‑1 (MQ=255)                                                                                                cTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAg                                                   >  1:140878/1‑149 (MQ=255)                                                                                                   gTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCt                                                >  4:271651/1‑149 (MQ=255)                                                                                                    tCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCtt                                               <  4:208382/149‑1 (MQ=255)                                                                                                      gAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTtgc                                             >  8:40784/1‑149 (MQ=255)                                                                                                      gAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTtgc                                             >  2:194119/1‑149 (MQ=255)                                                                                                      gAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTtgc                                             >  2:240888/1‑149 (MQ=255)                                                                                                      gAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTtgc                                             >  2:85093/1‑149 (MQ=255)                                                                                                      gAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTtg                                              <  6:154032/148‑1 (MQ=255)                                                                                                             ctgGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGt                                                       <  3:244706/132‑1 (MQ=255)                                                                                                                gAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGcc                                     >  4:227639/1‑147 (MQ=255)                                                                                                                gAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGa                                   <  1:194119/149‑1 (MQ=255)                                                                                                                gAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGa                                   <  6:255532/149‑1 (MQ=255)                                                                                                                    aCGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGAGTTg                               <  7:40784/149‑1 (MQ=255)                                                                                                                     cGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGAGTTGc                              <  7:146004/149‑1 (MQ=255)                                                                                                                        ccTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGAGTTGcaca                           <  5:178694/149‑1 (MQ=255)                                                                                                                        ccTG‑AAAAACACCGACCAG‑GGCGAGATCCGGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGAGTTGca                             >  5:293274/1‑147 (MQ=255)                                                                                                                         cTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTtgc                                             >  5:305279/1‑130 (MQ=255)                                                                                                                         cTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTtgc                                             <  6:305279/130‑1 (MQ=255)                                                                                                                          tG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGAGTTGCACAGc                         >  1:19965/1‑149 (MQ=255)                                                                                                                               aaaCACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGAGTTGCACAGccacc                     <  5:203152/149‑1 (MQ=255)                                                                                                                                  cacCGACCAG‑GGCGAGATCCGGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGAGTTGCACAGCCACCAtt                  <  3:154689/149‑1 (MQ=255)                                                                                                                                   acCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCa                                                                                         <  5:314040/77‑1 (MQ=255)                                                                                                                                   acCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCa                                                                                         >  6:314040/1‑77 (MQ=255)                                                                                                                                    ccGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGAGTTGCACAGCCACCAtttc                <  4:28754/149‑1 (MQ=255)                                                                                                                                      gACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGAGTTGCTCAGCCACCAtttcct              <  3:80204/149‑1 (MQ=255)                                                                                                                                      gACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGAGTTGCTCAGCCACCAtttcct              <  3:75566/149‑1 (MQ=255)                                                                                                                                           g‑ggCGAGATCCTGCTGGTGGTGGCGCTGGTTCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGAGTTGCACAGCCACCATTTCCTTTCCa         <  1:44202/149‑1 (MQ=255)                                                                                                                                           g‑ggCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGAGTTGCACAGCCACCATTTCCTTTCCa         <  2:140878/149‑1 (MQ=255)                                                                                                                                           g‑ggCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGAGTTGCACAGCCACCATTTCCTTTCCa         <  3:227639/149‑1 (MQ=255)                                                                                                                                                gagaTCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGAGTTGCACAGCCACCATTTCCTTTCCAGcaa     >  3:244888/1‑149 (MQ=255)                                                                                                                                                   aTCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGAGTTGCACAGCCACCATTTCCTTTCCAGcaacaa  <  3:302020/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                 GCCAAGGCCGAGCAACAGAACATCGAGCTGATCAGCTTCACCCAGCCACAAGTACCGCGGGTCATCAGCGGAGACCCGACGCGTCTGCGCCAGGCCTTGTCGAGCCTGCTGGAAAACGCCCTG‑AAAAACACCGACCAG‑GGCGAGATCCTGCTGGTGGTGGCGCTGGATCAACGCGGCGAGGTGCCACGCCTGCGTATCGCCGTGCAGGACAGCGGCGAGCCGTTGCCCGCCGCCGACCGTGAAGCCTTGCTGCAGGCCGAGTTGCACAGCCACCATTTCCTTTCCAGCAACAA  >  NC_002947/5488127‑5488419

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