Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,542,927	G→C	12.6%	G200G (GGC→GGG)	cysP ←	thiosulfate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,542,927	0	G	C	12.6%	177.7 / 2.3	72	G200G (GGC→GGG)	cysP	thiosulfate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (51/11);  new base C (2/8);  total (53/19)
Fisher's exact test for biased strand distribution p-value = 2.04e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGG  >  NC_000913/2542696‑2543165                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               ttGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGa                                                                                                                                                                                                                                          >  1:17166/1‑238 (MQ=255)    cgcCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTAAGGATGTTGTTCACTTCCGATTCGAAGCTAAGCAGCACATCGCCCAGGCCGCGCTC‑GA‑CAAAAGTGTCGGCCGCGCCACGa                                                                                                                                                                                                                                          >  2:400447/1‑235 (MQ=255)                      ccAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TAACCACTTCAAAGCCCTGCGCTTCATACAGTTTACGTATGTTGTTCACTTCCGTTTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGAGCGCCAGTCTCGACCACTTCa                                                                                                                                                                                                                     >  2:929625/1‑238 (MQ=255)                       cAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAAc                                                                                                                                                                                                                   >  1:554585/1‑239 (MQ=255)                         gTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGACGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCAAGTATCGAACACTTCAACGt                                                                                                                                                                                                                 >  2:1157159/1‑239 (MQ=255)                         gTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGt                                                                                                                                                                                                                 >  1:686400/1‑239 (MQ=255)                           tCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCCCCAGTATCGACCACTTCAACGttt                                                                                                                                                                                                               >  1:216217/1‑239 (MQ=255)                               aTAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCACCGTTTTTCa                                                                                                                                                                                                           >  1:595552/1‑239 (MQ=255)                                    cTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAggc                                                                                                                                                                                                        >  1:158921/1‑236 (MQ=255)                                    cTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCCC‑GG‑CAAAAGTGGTGGTCGCCCCACGACCGCCCGTATCGAACACTCCAACGTTTTTCAgta                                                                                                                                                                                                        >  1:694156/1‑235 (MQ=255)                                     tttGGCGTTTTTTCCCGTTCCGTGGGCCTGCTCGTTTTTACCAACCCACGCCCCCGGGAT‑TCCCGCCAGATTGTTGTTTTT‑GGGAC‑TCCCCACTTCAAGCCCCTGCGCTTCATACTGTTTACGGTTGTTTTTCCCTTCCGATTCGAGCCTAATCACCACACCGCCCAGCCCGCGCTC‑GG‑CAAAGGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACt                                                                                                                                                                                                     <  2:217010/239‑1 (MQ=255)                                             cTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATa                                                                                                                                                                                             >  1:217010/1‑239 (MQ=255)                                             cTTTTCCCGACCCGTTGCCCGGCGTG‑TTTTAACCACCCCCCCCCCCGGGAA‑TTCCCGCAAAATGTTGGTTTTGCGAAA‑CCCCCCCTCCAAACCCCTGCGTTTCATCCTTTTGAGGGGTGTTGTCCCCTTCCGTTCCAAAGCTAACCACCACCCCGCCAGGGCCGCCCTC‑GG‑AAAAGGTGGGGGCCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATa                                                                                                                                                                                             <  2:17166/239‑1 (MQ=255)                                                   cggTACGGTTGGCTGGCTCGTTTTTACCACCCCACG‑CGTAGTGGA‑TTCCCACAAGATGTTGTGTTG‑CGGAACTCCCCACTCCAAAGCCGTGCGCTCCATACTTTTTGCGGATGTTTTACACTTCAGATCCGATGTTAATCAGCCCCTCGCCAAGGCCGGGCTA‑GG‑CAAAGGGGGGGTACGGGCCGCTACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGt                                                                                                                                                                                       <  2:203823/237‑1 (MQ=255)                                                       aCCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CCAAAGGTGTGGTCGCGCCACGACCGCCAGTAACGAACACTTCACCGTTTTTCAGGAACTTGGTCATACACTGCTCgg                                                                                                                                                                                   >  2:450392/1‑239 (MQ=255)                                                         ggtggtctggcgcttttttACAAACCCAC‑CCAACGGGAAGTCCCGCCAGATTTTTGTTTTA‑GGAAA‑TCCCCACTCCAAACCCCTGCGCTTAATCCTGTTTCCGGAGTTTGTCCCCTCCCGATTCGAAGCAACCCACCACTCCCCCCGGGCCGCGCTC‑GG‑CAAAAGGGGGGGCCGCCCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGtt                                                                                                                                                                                 <  2:21099/225‑1 (MQ=255)                                                           ttGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGtttt                                                                                                                                                                               >  1:540823/1‑239 (MQ=255)                                                                   cACGTTTTTATCAACCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGGTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTACTCAGCACATCGCCCAGGCCGCGTTC‑GG‑CCAAGTTGTAGGGCTCCCTAAGACCCCCATTACCCAAAATTTCACGGATTTCAAGGAACTGGCTCATGAACTTTTC‑GTTTTGCCTTTGtt                                                                                                                                                                      >  2:347343/1‑238 (MQ=255)                                                                         tttATCAACCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTcacca                                                                                                                                                                  >  1:624045/1‑238 (MQ=255)                                                                           caaccaCCCACTCCAGCTGAAA‑TTCCGCTGGTATGGTTGTTTC‑GTGAA‑TCCCCACTCCATGCCCCGTCGCTTCATACTGTTTGCGGTTGTGTGTCTGTTCCAGTGCGACGCAACACCGCACACCGCCCAGGCCGCGCTC‑GC‑CAAAAGGGGTGGTCGGCCAATGTCCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCg                                                                                                                                                               <  2:312797/234‑1 (MQ=255)                                                                            accaccccccccccccGGGAA‑TGCCGCGAGACTGTTGGTTGC‑CGGAA‑TCACCCCTTAAACGCCCGCCGCTTCATCTTGTTTCCGGTGGTTTTTCCCTTCCGATGCGAAGCTAACAAGCACACCCCCCGGGCCGCGCTG‑GG‑CAAAAGTGGTGGGCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGt                                                                                                                                                              <  2:872276/225‑1 (MQ=255)                                                                                aCCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCCGCACATCACCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGGAACGAACACTTCAACGGTTTTCAGGAACTGGGTCGTAACCTGGTCGGGTTCGCCTTTGGCACTCCCGACAga                                                                                                                                                          >  2:734788/1‑238 (MQ=255)                                                                                  ccACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGCGGTGGTCGCGCCACGACCGACAGTATCCCACACTTCAACGTTTTTCAGGAAGTGGGTCATAAACTTTTAGGTTTTTCCTTTGTTACCACCGTCACCtt                                                                                                                                                        >  2:874529/1‑239 (MQ=255)                                                                                     cGCCACCGGGAC‑TTCCGCCAGCCTGTTGGTTTT‑CGGCA‑TCACCACTTCAAAGCCCTGCGCTTCCTACTGTTTACGGATGTTGTTCACTTCCGCTTCGAAGCTACTCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACCCCCGCCCGAATCGAACACTTCCACGTTTTTCAGCAACTGGCTCATAAACTGTTCGGCTTTGCCTTTGTCACCACCGTCCGCTTTa                                                                                                                                                      >  2:396820/1‑238 (MQ=255)                                                                                      gCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTAt                                                                                                                                                     >  1:872276/1‑238 (MQ=255)                                                                                      gCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTAt                                                                                                                                                     >  1:783230/1‑238 (MQ=255)                                                                                      gCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTAt                                                                                                                                                     >  1:1051156/1‑238 (MQ=255)                                                                                           cGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCt                                                                                                                                                >  1:137618/1‑238 (MQ=255)                                                                                                   tCCGCAAGAATGTTGTTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGTTGTTGTCCACTCCCGTTTCAAAGCTAACCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCGCGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCa                                                                                                                                        <  1:588952/239‑1 (MQ=255)                                                                                                     cGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCCGCTTTATCCGCTGCGCCCCAt                                                                                                                                       >  1:96643/1‑238 (MQ=255)                                                                                                     cGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTA‑GG‑CAAAAGTGGTGGTCGAGCCACGACCGCCAGTAACGAACACTTCAACGTTTTTCAGGAACTGGGCCATAAACTGTTCGGTTTTGCCTTTGTCACCAACGTAAGCTTTGTCCGCTCCGCCCCAt                                                                                                                                       >  2:247612/1‑238 (MQ=255)                                                                                                           aaTTTTGTTTTT‑GGGAA‑CCCCCCTTACAAAGCCCTGCTTTTCTTACTTTTACCGGTTTTTTTTCCCTCCCAATCGGAAGATAATCCGCCCACCGCCCGGGCCGCGCCG‑GG‑AAAAAGGGGTGGCCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCa                                                                                                                                <  1:580878/239‑1 (MQ=255)                                                                                                                 ggTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGGTCACTTCCGATTCGAAGCTAATAAGCACATCGCCCAGGCCGCGCTC‑GG‑CCAAAGTGGTAGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCCTAAACTGTTGGGATTTCCCTTTGTCCCCACCCGCCGCTTTGTCCCATCCGCCCCCTgccccccactag                                                                                                                           >  2:614595/1‑232 (MQ=255)                                                                                                                        cGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAAc                                                                                                                    >  1:1145598/1‑239 (MQ=255)                                                                                                                                      tCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCg                                                                                                        >  1:303202/1‑238 (MQ=255)                                                                                                                                          aGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCAACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGTTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGATTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACCCGCGttgaccaacg                                                                                                     >  2:918211/1‑229 (MQ=255)                                                                                                                                              gtgcgtcatttccTGTTTACGGTTGTTGTTCTCTTCCGGTAACAACCAACCCATCACACCGCCCAGGCCGCGCTCTGGACAAAGGTGGTGCGCGGCCCCCGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGtttt                                                                                                 <  2:807324/227‑1 (MQ=255)                                                                                                                                               tGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAATGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcgg                                                                                              >  2:837528/1‑239 (MQ=255)                                                                                                                                               tGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGc                                                                                                                                                                             >  1:432480/1‑160 (MQ=255)                                                                                                                                                 cgcTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcgggt                                                                                            >  1:1152486/1‑239 (MQ=255)                                                                                                                                                         ccgttttgcggagtgtgttCTCTCCCGAGTCGAGGCTAACCAGCCCCTCGCCCAGGCCGCGTCC‑GG‑AAAAGGTGGTGGCCGCCCCACCACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGaaa                                                                                    <  2:1145598/224‑1 (MQ=255)                                                                                                                                                           tttttCCGGATGTTTTTCACTCCCATTTCGAAGCTAATCAGCACATCGCCCGGGCCGCCCTC‑GG‑AAAAGGTGGGGTGCCCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAAt                                                                                  <  1:722577/237‑1 (MQ=255)                                                                                                                                                             tttACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CCAAAGTGGTGGTCGCGCCACGACAGCCAGTATCGAACACTTCCACGTTTTTCAGGAACTGGATCATAAACTGTTCGGTTGTGCCTTTCTCACCACCGTGAGCTTTATCCGCTGAGCCCCATTCCGCCAGATAGGTATAACACGCGTTACCCGACGTTTTCGTGTTCTGGAGAAtga                                                                                >  2:234891/1‑237 (MQ=255)                                                                                                                                                              ttACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGaa                                                                                                                                                                                                       <  1:364050/119‑1 (MQ=255)                                                                                                                                                              ttACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGaa                                                                                                                                                                                                       >  2:364050/1‑119 (MQ=255)                                                                                                                                                              ttACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCa                                                                                >  1:558624/1‑238 (MQ=255)                                                                                                                                                              ttACGGATGTTGCTCACTTCCGATTCGCAGCTAATCAGCACATCGCCCAGTCCGCGCTC‑GG‑CAAAAGCAGTGGTCCCGCCACGACCACCAGTATCGAACACTTCAACGTTCTTCTGGAACTGTGTCATACACTCTTCGGTCTTGCCTGTGTCACTCCCGTCCGCTTTATCCGCTGCGACCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGAAGGTTTCGGGTTCGGGAAAGTCa                                                                                >  2:722577/1‑238 (MQ=255)                                                                                                                                                                aCGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCt                                                                             >  1:883211/1‑239 (MQ=255)                                                                                                                                                                     tgttgtTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACg                                                                        >  1:1101882/1‑239 (MQ=255)                                                                                                                                                                      ttttttCATTTCCGATTCAAAGCTAACCAGCACACCGCCCGGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCCCCACAGCCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGt                                                                       <  1:1157159/235‑1 (MQ=255)                                                                                                                                                                       ttgttCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CACAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGTTAGGTATAATGCGCGTTACCCGACGTCTTCGGGTTCGGGACAATCAGCTTAACCTc                                                                      >  2:1048787/1‑239 (MQ=255)                                                                                                                                                                       ttgttCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCTGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGTCTCGGGAAAATAAGCTTCCCGTc                                                                      >  2:826276/1‑239 (MQ=255)                                                                                                                                                                       ttgttCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTc                                                                      >  1:34371/1‑239 (MQ=255)                                                                                                                                                                       ttgttCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGGCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGCCATAAACTGTTCGGTTTTGCGTTTGTCAACACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACACGCGTTCCCCGACTTTTTAGGGTTAGGGaaattactctgcacgtc                                                                      >  2:633937/1‑225 (MQ=255)                                                                                                                                                                         gttCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACACGCGTTACCCGACGTTTTCGGGTTCGGTAAAATCACCTTCACGTCgg                                                                    >  1:269804/1‑239 (MQ=255)                                                                                                                                                                            cACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAgc                                                                 >  2:268668/1‑239 (MQ=255)                                                                                                                                                                               ttCCGGTTCGAAGCTAATCAGCACGTCCGCCAGGCCGCGCCC‑GG‑ACAGAGGGGGGGGCGCGCCACGGCCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCa                                                              <  1:614595/239‑1 (MQ=255)                                                                                                                                                                                 ccGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCAcc                                                            >  1:564395/1‑239 (MQ=255)                                                                                                                                                                                    aTTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCCGGTACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGCAAATCAGCTTCACGTCGCAGCGCCCCAgg                                                         >  1:525469/1‑239 (MQ=255)                                                                                                                                                                                    aTTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAg                                                          >  1:280443/1‑238 (MQ=255)                                                                                                                                                                                            cTAATCCGCCCCCCGCCCAGGCCGCGCTC‑GG‑CAAAAGGGGGGGTCGCCCCACGTCCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCa                                                 <  1:929625/239‑1 (MQ=255)                                                                                                                                                                                                  aGCACACCGCCCGGGCCGCGCTG‑GG‑AAAAAGTGGTGGTCGCCCCACGACCGCCAGTATCGAACACTTCa                                                                                                                                                                                                                     <  1:889319/69‑1 (MQ=255)                                                                                                                                                                                                      caTCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGata                                       >  1:975208/1‑239 (MQ=255)                                                                                                                                                                                                           ccaaGGCGGCGCCC‑GA‑AAAAGGGGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCg                                                                     <  1:627167/201‑1 (MQ=255)                                                                                                                                                                                                            ccAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCa                                                                                                                                                                                                           >  2:889319/1‑69 (MQ=255)                                                                                                                                                                                                                 ccGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCCACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGtt                            >  2:445833/1‑239 (MQ=255)                                                                                                                                                                                                                   tcgcTC‑GG‑AAAAAGTGGTGGACGCGACTCGACCTACAGTATCGCACACTTCAACGTTTTTCATTAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCAGCTGCGCCCGATGCAGCCAGATAGGTATATCGCCCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCCGGTAGTTCCAATACAGGATAGTCTTCGGCTTa                           >  2:125126/2‑238 (MQ=255)                                                                                                                                                                                                                   cccgccgcagcaagtttagtgttCGCCCCACCTCCGCCAGTATCGAACACTTCAGCGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCTCCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGt                             <  2:129590/220‑1 (MQ=255)                                                                                                                                                                                                                          gg‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTAc                    >  1:938329/1‑239 (MQ=255)                                                                                                                                                                                                                                    gttgggCGCGCAACGGCCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAgg            <  1:874529/236‑1 (MQ=255)                                                                                                                                                                                                                                       ggtCTCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGGTCGGTTTTGCCTTTGTCACCACCGCCAGCTTTATCAGCTGAGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGTTTCGGGAAAATCAGCTTCACGTCGGGGCGCACCAGGTCGTTCCACTCGGGGATATTCTTCGGGTTACCCTTACGCCCGAGGatg         >  2:128736/1‑236 (MQ=255)                                                                                                                                                                                                                                        gtCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATAttctt                                  >  1:930559/1‑212 (MQ=255)                                                                                                                                                                                                                                            cgcCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTATTCGGGTTACCCTTACGCACCAGGCTGCCCAtg   >  1:288086/1‑239 (MQ=255)                                                                                                                                                                                                                                             gcCACGACCGCAAGTATCGAACACTTCAACGTTTTGCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCAGCTGCGCCCCATGCCGCCAGATAGGTATACCGCGCGTTACCCGACGTTTGCGGGTTCGGGAAAATCAGATTCACGTCGGATCGCACCCGGTAGTTCCAATCGTAGTTATTATTCGGGTTACCCTTACGCACCAGCAGGCCCAtgg  >  2:964899/1‑239 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               TTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAA‑TTCCGCCAGAATGTTGGTTTT‑CGGAA‑TCACCACTTCAAAGCCCTGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTC‑GG‑CAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGG  >  NC_000913/2542696‑2543165

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