Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,218,828	T→A	9.9%	E158V (GAA→GTA)	nfxB ←	HTH‑type transcriptional regulator NfxB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,218,828	0	T	A	9.9%	254.6 / 14.2	90	E158V (GAA→GTA)	nfxB	HTH‑type transcriptional regulator NfxB
Reads supporting (aligned to +/- strand):  ref base T (41/40);  new base A (5/4);  total (46/44)
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.98e-01

TTGTCCAGCCACTTGGCCCCAGGG‑ATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGCATCCAGGGCGTCGGCATAGGCCCGCCACGGGCGGTTCGTCT  >  NC_002947/3218681‑3218975                                                                                                                                                      |                                                                                                                                                      ttGTCCAGGCACTTGGCCCCAGGG‑AGCTATGCCCTGGGGTTTTGTTGGCTTGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCGCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCg                                                                                                                                                     <  2:574921/150‑1 (MQ=255)  tGTCCAGCCACTTGGCCCCAGGG‑ATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCgg                                                                                                                                                    <  1:411031/150‑1 (MQ=255)      cAGCCACCTGGCCCCAGGG‑ATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCTGCATc                                                                                                                                                <  1:974199/150‑1 (MQ=255)       aGCCACTTGGCCCCAGGG‑ATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATc‑g                                                                                                                                              >  2:1252016/1‑150 (MQ=255)        gCCACTTGGCCCCAGGG‑ATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TACGGCATc‑gt                                                                                                                                             >  2:636312/1‑149 (MQ=255)         ccACTTGGCCCCAGGG‑ATCGATGCCCTGGGGTTTTTTTGGCCGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATc‑gac                                                                                                                                            >  2:1830375/1‑150 (MQ=255)          cACTTGGCCCCAGGG‑ATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATc‑gacg                                                                                                                                           >  2:258796/1‑150 (MQ=255)           aCTTGGCCCCAGGG‑ATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGTCGCCCGCCCGCGTCG‑TTCGGCATc‑gacga                                                                                                                                          >  2:2123771/1‑150 (MQ=255)           aCTTGGCCCCAGGG‑ATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATc‑gacga                                                                                                                                          <  2:434013/150‑1 (MQ=255)           aCTTGGCCCCAGGG‑ATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTAACGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATc‑gacga                                                                                                                                          >  1:2262945/1‑150 (MQ=255)             ttGGCCCCAGGG‑ATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCTGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATg                                                                                                                                        >  1:1725273/1‑150 (MQ=255)                  cccAGGG‑ATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCAt                                                                                                                                                 <  2:430563/137‑1 (MQ=255)                  cccAGGG‑ATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCAt                                                                                                                                                 >  1:430563/1‑137 (MQ=255)                  cccAGGG‑ATCGATGCACGGAGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTa                                                                                                                                   >  1:2302259/1‑150 (MQ=255)                   ccAGGGCATCGAT‑CCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCg                                                                                                                                                     >  2:799725/1‑132 (MQ=255)                   ccAGGGCATCGAT‑CCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCg                                                                                                                                                     <  1:799725/132‑1 (MQ=255)                   ccAGGG‑ATCGATGCCCTGGGGTTTTTTTTTCTTTTTAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGTCGTCCCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑CACGATGCAGtcc                                                                                                                                  >  2:838944/1‑148 (MQ=255)                   ccAGGG‑ATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAACGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TACGGCATC‑GTCGATGCCGAac                                                                                                                                  <  1:636312/150‑1 (MQ=255)                   ccAGGG‑ATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TACGGCATC‑GTCGATGCCgaa                                                                                                                                   >  1:860105/1‑147 (MQ=255)                   ccAGGG‑ATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TACGGCATC‑GTCGATGCCgaa                                                                                                                                   <  2:860105/149‑3 (MQ=255)                      ggg‑ATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACAcc                                                                                                                               >  2:985683/1‑150 (MQ=255)                             gATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCAAACGGGGTCCTCATCAACGTCGAGGAAAAGCAGCTCGTTCTTTATCTGTGAGTTCGCGCTTGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GAAGATGCAGTACACCATGC‑Tc                                                                                                                        >  1:1112961/1‑150 (MQ=255)                              aTGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCa                                                                                                                       <  1:963212/150‑1 (MQ=255)                                 cccTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGa                                                                                                                    <  1:1904234/150‑1 (MQ=255)                                  ccTGGGGTTTTGTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGaa                                                                                                                   <  1:717788/150‑1 (MQ=255)                                         tttttttGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCTCGGACGTCGCGCTGGCCGCGCGGACG‑GTAGAAATGCCGGC‑TACGATGCCGTACACCATGC‑TCAGGAACATTTCg                                                                                                            >  1:1686091/1‑150 (MQ=255)                                           tttttGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTGCgtt                                                                                                          <  2:2262945/150‑3 (MQ=255)                                             tttGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGa                                                                                                        <  1:1890053/150‑1 (MQ=255)                                              tgggCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGaa                                                                                                       <  1:1699510/148‑1 (MQ=255)                                                 gCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGAt                                                                                                    >  2:1916421/1‑150 (MQ=255)                                                  cTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATc                                                                                                   <  2:2183789/150‑1 (MQ=255)                                                   tGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATcg                                                                                                  >  1:1408276/1‑150 (MQ=255)                                                    ggTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATcgc                                                                                                 <  2:691800/150‑1 (MQ=255)                                                     gTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATcgcc                                                                                                >  2:1945073/1‑150 (MQ=255)                                                      tGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGACA‑AT‑GGCGCCTGACGGTGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATcgccg                                                                                               >  2:653063/1‑150 (MQ=255)                                                        aaGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TACGGCATC‑GTCGATGCCgaac                                                                                                                                  >  2:710681/1‑111 (MQ=255)                                                        aaGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TACGGCATC‑GTCGATGCCGAac                                                                                                                                  <  1:710681/114‑1 (MQ=255)                                                            cAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGCACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGa                                                                                         <  1:1830375/150‑1 (MQ=255)                                                             aGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGat                                                                                        >  2:1577866/1‑150 (MQ=255)                                                                    cgagcgGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATa                                                                                 <  2:1408276/150‑1 (MQ=255)                                                                    cgagcgGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TACGGCATC‑GTCGATGCCga                                                                                                                                    <  1:92511/100‑2 (MQ=255)                                                                    cgagcgGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TACGGCATC‑GTCGATGCCga                                                                                                                                    >  2:92511/1‑99 (MQ=255)                                                                     tcgacggTGCTGCGGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGTCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATAc                                                                                >  2:1948820/6‑150 (MQ=255)                                                                     gagcgGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATAc                                                                                <  1:838267/150‑1 (MQ=255)                                                                       gcgGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACgg                                                                              <  1:838944/150‑1 (MQ=255)                                                                           gTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAg                                                                          <  1:258796/150‑1 (MQ=255)                                                                                gcagcaCCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACAttt                                                                                                              <  2:1837268/109‑1 (MQ=255)                                                                                gcagcaCCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACAttt                                                                                                              >  1:1837268/1‑109 (MQ=255)                                                                                     aCCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCctg                                                                <  1:1252016/150‑1 (MQ=255)                                                                                      ccGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCAGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGCAGATCGCCGCGCTGATATCCATACGGAAAGCCCCCGCctgc                                                               >  2:1007154/1‑150 (MQ=255)                                                                                       cggcGAGGAAAAGCAGCGCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCctgct                                                              <  2:1112961/147‑1 (MQ=255)                                                                                       cGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGa                                                                                                     <  1:720726/111‑1 (MQ=255)                                                                                       cGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGa                                                                                                     >  2:720726/1‑111 (MQ=255)                                                                                             ggAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTCCACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCcgc                                                        >  1:1102594/1‑150 (MQ=255)                                                                                              gAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCcgcg                                                       >  2:2105077/1‑150 (MQ=255)                                                                                                  agcagcTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAgg                                                   <  1:2010875/150‑1 (MQ=255)                                                                                                   gcagcTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATcgcc                                                                                                <  1:1951998/104‑1 (MQ=255)                                                                                                   gcagcTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATcgcc                                                                                                >  2:1951998/1‑104 (MQ=255)                                                                                                    cagcTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGCGTCA‑ACAGAGGGCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGgaa                                                 >  1:735025/1‑150 (MQ=255)                                                                                                        tCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGaagaag                                             <  1:985683/150‑1 (MQ=255)                                                                                                        tCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGaagaag                                             >  1:675121/1‑150 (MQ=255)                                                                                                        tCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGaagaag                                             >  1:2254043/1‑150 (MQ=255)                                                                                                          cAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACAcc                                                                                                                               >  1:1095183/1‑67 (MQ=255)                                                                                                          cAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACAcc                                                                                                                               <  2:1095183/67‑1 (MQ=255)                                                                                                          cAGCGTTTGCACGGACGTCGCGCCGGCCGCCCGCCCGCGTCG‑TTCGGCCTC‑GACGATGCCGTCCACCATGC‑TCAGGACAATTGCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGc                                           <  1:1916421/150‑1 (MQ=255)                                                                                                             cGTTTGCACGGACGTCACGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGCATc                                        >  1:2295301/1‑150 (MQ=255)                                                                                                                ttGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGCATCCAg                                     <  2:1686091/150‑1 (MQ=255)                                                                                                                  gCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGc                                           >  1:1802161/1‑142 (MQ=255)                                                                                                                  gCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGc                                           <  2:1802161/142‑1 (MQ=255)                                                                                                                    aCGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCa                                                     <  2:341850/130‑1 (MQ=255)                                                                                                                    aCGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCa                                                     >  1:341850/1‑130 (MQ=255)                                                                                                                    aCGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGc                                           <  1:1966463/140‑1 (MQ=255)                                                                                                                    aCGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGc                                           >  2:1966463/1‑140 (MQ=255)                                                                                                                      ggACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGCATCCAGGGCGTc                               >  1:2256686/1‑150 (MQ=255)                                                                                                                             acgcTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGCATCCAGGGCGTCGGCATAg                        <  2:2295301/149‑1 (MQ=255)                                                                                                                              cgcTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGCATCCAGGGCGTCGGCATAgg                       >  2:297334/1‑150 (MQ=255)                                                                                                                               gcTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGACGCCCAGGAAGAAGGCATCCAGGGCGTCGGCATAGGc                      >  2:1937246/1‑150 (MQ=255)                                                                                                                                      gcccgcccgcGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGCATCCAGGGCGTCGGCATAGGCCCGCCAc               >  2:1655705/1‑150 (MQ=255)                                                                                                                                        ccgcccgcGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGaagaag                                             >  2:1273241/1‑118 (MQ=255)                                                                                                                                        ccgcccgcGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGaagaag                                             <  1:1273241/118‑1 (MQ=255)                                                                                                                                           cccgcGTCG‑TTCGGCATC‑GACGATGCCGTACACCCTGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGCATCCAGGGCGTCGGCATAGGCCCGCCACGGGCg          <  2:1102594/150‑1 (MQ=255)                                                                                                                                             cgcgTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGCATCCAGGGCGTCGGCATAGGCCCGCCACGGGCggg        >  1:1809183/1‑149 (MQ=255)                                                                                                                                              gcgTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCGCCGCCTGCTGGCCGCGCAGGAAGAAGGCATCCAGGGCGTCGGCATAGGCCCGCCACGGGCGGtt       <  2:780525/150‑1 (MQ=255)                                                                                                                                               cgtcg‑tTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGCATCCAGGGCGTCGGCATAGGCCCGCCACGGGCGGTTc      <  1:1577866/150‑1 (MQ=255)                                                                                                                                               cgtcg‑tTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGCATCCAGGGCGTCGGCATAGGCCCGCCACGGGCGGTTc      <  1:1781557/150‑1 (MQ=255)                                                                                                                                               cgtcg‑tTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGCATCCAGGACGCCGGCATAGGCCCGCCACGGGCGGTTc      <  1:653063/150‑1 (MQ=255)                                                                                                                                                 tcg‑tTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGCATCCAGGGCGTCGGCATAGGCCCGCCACGGGCGGTTCGt    >  2:2033504/1‑150 (MQ=255)                                                                                                                                                   g‑tTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGc                                                                         >  2:2035587/1‑79 (MQ=255)                                                                                                                                                   g‑tTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGc                                                                         <  1:2035587/79‑1 (MQ=255)                                                                                                                                                   g‑tTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGCATCCAGGGCGTCGGCATAGGCCCGCCACGGGCGGTTCGtct  >  1:84878/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TTGTCCAGCCACTTGGCCCCAGGG‑ATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCG‑TTCGGCATC‑GACGATGCCGTACACCATGC‑TCAGGAACATTTCGGTGAAGATCGCCGCGCTGATATCGATACGGAAAGCCCCCGCCTGCTGGCCGCGCAGGAAGAAGGCATCCAGGGCGTCGGCATAGGCCCGCCACGGGCGGTTCGTCT  >  NC_002947/3218681‑3218975

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑