Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	2,372,052	G→A	13.2%	Q104Q (CAG→CAA)	arnC →	undecaprenyl‑phosphate 4‑deoxy‑4‑formamido‑L‑arabinose transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	2,372,052	0	G	A	13.2%	177.9 / 2.3	76	Q104Q (CAG→CAA)	arnC	undecaprenyl‑phosphate 4‑deoxy‑4‑formamido‑L‑arabinose transferase
Reads supporting (aligned to +/- strand):  ref base G (26/38);  new base A (2/8);  total (30/46)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCACCACCACAGCCTGTGAATCGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGAACGCAGCACC  >  NC_007779/2371824‑2372286                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           gcgcACCACCACAGCCTGTGAATCGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTTAGTCAAGTTACATGCGACTTGATTTTTA‑ACCCT‑TA‑TTTCGAT‑CTTCAGAAACcccccg                                                                                                                                                                                                                                  >  1:849734/1‑237 (MQ=255)           gcaGCCTGTGAAGCGTGGGGGAAGAAGTTTGAGTTCCTGTGTTTTAATGACGACTGTAGCGATAATTCCTCACATATGCTGGTCGAAGCCACACAGGCGGAGAACCGCCATATTGTGTCTATTTTGCTTAACCGCAATTCCGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCCCGTTACTGGCGACTTAATTATTC‑CCCTT‑GA‑TGCCGAT‑CTCAAGAATCCGCCAGAAGAAATccc                                                                                                                                                                                                                        <  2:49612/238‑1 (MQ=255)               ccTGTGAATCGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGGGTTAGTCTCGTTATTGGCGGATTTATTATTT‑CCCTT‑TT‑TTCCGAT‑CTCCCCAATCCCCCCCAAGAAAAcccccc                                                                                                                                                                                                                     >  1:132995/1‑237 (MQ=255)                 tgtgAATCGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGGGTCAGTCACGGTTCTGGCGCGCTTTTTATTA‑ACCTT‑TA‑TGTCGAG‑CTCCCCAAACCGCCCGAAGAAATaccccc                                                                                                                                                                                                                     >  1:1058144/1‑235 (MQ=255)                       tCGTTGGGGAAAGAGTATGTGTTCCTTCTTGTTGATGACGACAGTGGGGATAATTCCGCGCATATGCTGGTCGAAGCCTCAAAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGc                                                                                                                                                                                                            <  2:550303/239‑1 (MQ=255)                        cGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGGGTCTGTCACGATAACGGCGAGTTAATTTTTA‑TCCTT‑GA‑TGCCGAG‑CTGCAGAAACCGCCAGAAGAAAACCCCCGCCcggggggc                                                                                                                                                                                                           >  1:479249/1‑237 (MQ=255)                        cGTTGGGAAAAGAGTTTGAGATCCTGTTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGACCAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCa                                                                                                                                                                                                           <  1:1656005/239‑1 (MQ=255)                            ggggAAAGAGTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAg                                                                                                                                                                                                       <  1:1551295/239‑1 (MQ=255)                                 aaGAGTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGAt                                                                                                                                                                                                  <  1:1841697/239‑1 (MQ=255)                                 aaGAGTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGAt                                                                                                                                                                                                  <  1:948596/239‑1 (MQ=255)                                       aTGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGCGGGGATTTA‑TTACAGTACTTGCGACATAATTATTA‑CCCTT‑CATTGCCGTT‑TACCGAAATCCCCCAAAAAAAAACACCCCCCGGGGGGCaaaaccgaggatagt                                                                                                                                                                                             >  1:1660649/1‑228 (MQ=255)                                       aTGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGAGTTAGTCACGTTACTGGCGACTTTATTATTA‑TCCTT‑TT‑TGCCGAT‑CTCCAGAATCCGCCCGAAAAAATACCCCGCCTTGTGGCCACAACCGACTAAGGt                                                                                                                                                                                            >  2:992764/1‑239 (MQ=255)                                          agaTCCTGCTGATGGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTa                                                                                                                                                                                          <  1:824467/238‑1 (MQ=255)                                             tCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGTTGGTCGAAGCCTCACAAGCGGAAAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTacgac                                                                                                                                                                                      <  1:1448800/239‑1 (MQ=255)                                             tCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTacg                                                                                                                                                                                        <  1:992764/237‑1 (MQ=255)                                             actggttgattggtGACGACAGTAGCGATATTCCCGCGCTTTGGCTGTTCGAAGCCTCACAACCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTC‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTacgac                                                                                                                                                                                      <  2:1431033/233‑1 (MQ=255)                                               ctgctgATTGAGGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCATGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGt                                                                                                                                                                                    <  1:1671188/239‑1 (MQ=255)                                                tgctgATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACCTTATTGGCGACCTCATTTTTA‑CCCCT‑TA‑TTCCGAT‑CTCCCAAATTCCCCCGAAGAAAACCCCCCCCCCGTGGCCAAAACCCCTGAAGGGTacgacggg                                                                                                                                                                                   >  1:1015386/1‑237 (MQ=255)                                                   tgattgatGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGt                                                                                                                                                                                 <  1:204760/238‑1 (MQ=255)                                                   tgattgatGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGt                                                                                                                                                                                 <  1:530489/238‑1 (MQ=255)                                                   tgattgagGAGGGCAGTAGCGATATTTCCGCGATTATGCTGGTCGAGGCCTCAAAAGCGGAAAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGCGGCAAAAGCCGATGAAGGTTACGACGTGGt                                                                                                                                                                                 <  2:747934/238‑1 (MQ=255)                                                     attgatGACGGCAGTAGCGATAATCCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGCGGGGATTCAG‑CACCTTTAAGGGGGCTTTATTTTTATCCCTT‑AT‑TGCCG‑G‑CTCCCGAACCCCCCCCAAGAAAACCCCCCCCCCGCGGCAAAAAACCAAGGAGGGTAAGACGTgggggg                                                                                                                                                                              >  1:1431033/1‑235 (MQ=255)                                                      ttgatGACGGCATTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAgg                                                                                                                                                                              <  1:1589094/238‑1 (MQ=255)                                                      ttgatGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAgg                                                                                                                                                                              <  1:1867907/238‑1 (MQ=255)                                                       tgatgaCGGCGGTGGCAATAATTCCGCGCTTATGCTGGTGGAAGCCTCACAGGCGGAGAACAGCCTTATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTCCTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTacgac                                                                                                                                                                                      <  1:842699/229‑1 (MQ=255)                                                          tgaCGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTg                                                                                                                                                                         <  1:1328721/239‑1 (MQ=255)                                                              ggCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAgg                                                                                                                                                                                             <  1:1449430/215‑1 (MQ=255)                                                                                                   cccaaaaggaggAGACCGGCGATTTTGTTTTTTTTTTGCATAACCGCAATTCGGGGCAATGTT‑A‑GCAGATTCGGGG‑GGGTTTAGGTCGCG‑TATCGGGGATTTTATTAATT‑CCCTTGGA‑CGCCGAC‑CCCAGGAATCCGCAAAAAGAAACCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCt                                                                                                                                <  2:1060851/229‑1 (MQ=255)                                                                                                                 aaCACCCATATTGTTTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCATTCACGTTACGGGCGACTTAATTATTC‑CCCTT‑GA‑TCCCGAT‑CTCCAGAACCCGCCAGAAGAAACCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTaa                                                                                                                  <  1:1858817/239‑1 (MQ=255)                                                                                                                    aGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCtg                                                                                                                                                                                                                 <  1:294823/141‑1 (MQ=255)                                                                                                                      ccttttggtgtTTATTTTGCTTACCCGAAATTCGGGACAACATTCA‑GC‑TATTAGGGG‑GGGACACATTCACGTTACTGGCGATTTATTTATTC‑CCCTG‑GA‑‑GCCCATGCCCCAAACCCCCCCAGAAAAAACCCCCCGCCTGGTGGCAAAAGCCGATGAAGTTTACTGCGTGGTAGGGACTGTACGCAAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGc                                                                                                             <  2:1353815/232‑1 (MQ=255)                                                                                                                         tatTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTT‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCaaa                                                                                                                                                                                                         <  1:998679/144‑1 (MQ=255)                                                                                                                               ttgtATTTTGCTGAACC‑AAATTTCCGGCAACTTTGACGC‑TATTGGGGG‑GGGTTCCACACGCGTTCGTGGCGTTTTATTTATTC‑CCCTT‑GT‑GGCTGTT‑CCCAAAAACCCGCCAAAAAAATCCCCCCGCCTGGTGGCAAAGGCCGATAAGGTTTACGACGTGTTAGGGACTGTACGCCAGACCCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAgc                                                                                                    <  2:963324/236‑1 (MQ=255)                                                                                                                                                gCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCGCGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCg                                                                                   <  1:45139/239‑1 (MQ=255)                                                                                                                                                  aaTTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTC‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGAt                                                                                 <  1:284739/239‑1 (MQ=255)                                                                                                                                                      aCGGGCAACATTCA‑GG‑GTTTATGGC‑GGGATTCAGTCACGTTCCTGGGAACTTAATTATTC‑CCCTT‑GA‑TGCCGAT‑CTCCAAAATCCGCCAAAAGAAACCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGt                                                                             <  1:73459/239‑1 (MQ=255)                                                                                                                                                           cAACATTCA‑GC‑GATTATGGG‑GGGATTCAGTCCCGTTACTGGCGACTTAATTATTC‑CCCTT‑GA‑TGCCGAT‑CCCAAAAACCCGCCAGAAAAAACCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTa                                                                        <  1:1168750/239‑1 (MQ=255)                                                                                                                                                            aaCATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTC‑CCCTT‑GA‑TGCCGAT‑CTCCAGAACCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTa                                                                        <  1:1129149/238‑1 (MQ=255)                                                                                                                                                                                     ttccattCACTTTACTGGCGACTTAATTATTC‑CCCTG‑GA‑GGCCGAT‑CCCCAGAACCCGCCAGAAGAAACCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTAt                                                  <  1:1892224/232‑1 (MQ=255)                                                                                                                                                                                      ttcgttcccGTTCCTGGCGATTTAATTATTC‑CCCTT‑GA‑TGCCGAT‑CTCCAGAACCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCg                                                <  1:2709/231‑1 (MQ=255)                                                                                                                                                                                              cGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCata                                        >  1:188663/1‑239 (MQ=255)                                                                                                                                                                                                ttACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATAtt                                      >  2:1351176/1‑239 (MQ=255)                                                                                                                                                                                                  aCTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTg                                     >  1:923282/1‑238 (MQ=255)                                                                                                                                                                                                        gatttAATTATTC‑CCTTT‑GA‑TGCCGTT‑CTCCAGATTCCGCCAGAAGAAACCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTCCGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTcgatg                               <  1:239815/235‑1 (MQ=255)                                                                                                                                                                                                         aCTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATg                                                                                >  2:985108/1‑188 (MQ=255)                                                                                                                                                                                                         aCTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATg                                                                                <  1:985108/188‑1 (MQ=255)                                                                                                                                                                                                         aCTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTcgatgcg                             >  1:1814201/1‑239 (MQ=255)                                                                                                                                                                                                         aCTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTcgatgcg                             >  1:676806/1‑239 (MQ=255)                                                                                                                                                                                                         aCTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTcgatgcg                             >  1:1774763/1‑239 (MQ=255)                                                                                                                                                                                                          cTTAATTATTC‑CCCTT‑GA‑TGCCGAT‑CCCCAAAACCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTcgatgcga                            <  1:1861650/239‑1 (MQ=255)                                                                                                                                                                                                                 atta‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAACCCCCCGCCTGGGGGCAAAGGCCGATGAAGGTTCCGACGTTGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCa                     <  1:218196/239‑1 (MQ=255)                                                                                                                                                                                                                  tta‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAg                                                                                                                          >  1:295021/1‑137 (MQ=255)                                                                                                                                                                                                                    a‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGcc                                              >  2:384964/1‑211 (MQ=255)                                                                                                                                                                                                                    a‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGcc                                              <  1:384964/211‑1 (MQ=255)                                                                                                                                                                                                                    a‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACt                   >  1:921898/1‑238 (MQ=255)                                                                                                                                                                                                                    a‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACt                   >  1:55916/1‑238 (MQ=255)                                                                                                                                                                                                                    a‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACt                   >  1:207690/1‑238 (MQ=255)                                                                                                                                                                                                                      tcccctgtttcccGAT‑CCCCAAAACCCCCCAAAAAAATCCCCCCGCCTGGTGGCAAACACCGATGAAGGTTACGCCGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACt                   <  1:1036206/228‑1 (MQ=255)                                                                                                                                                                                                                         tt‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCAGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCAt              >  1:1911321/1‑239 (MQ=255)                                                                                                                                                                                                                         tt‑GA‑TGCCGAT‑CCCCAGAATCCGCCAGAAGAAACCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCAt              <  1:382081/239‑1 (MQ=255)                                                                                                                                                                                                                          t‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATg             >  1:431251/1‑239 (MQ=255)                                                                                                                                                                                                                            gA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGa                                                                                                                                                                            <  1:675142/79‑1 (MQ=255)                                                                                                                                                                                                                            gA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGa                                                                                                                                                                            >  2:675142/1‑79 (MQ=255)                                                                                                                                                                                                                            gA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTAtt                                                                                                        >  2:1325038/1‑147 (MQ=255)                                                                                                                                                                                                                            gA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGa            >  1:657903/1‑239 (MQ=255)                                                                                                                                                                                                                            gA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGa            >  1:555212/1‑239 (MQ=255)                                                                                                                                                                                                                            gA‑TGCCGAT‑CTCCAGAATCCGCAAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTAtt                                                                                                        <  1:1325038/147‑1 (MQ=255)                                                                                                                                                                                                                            gA‑TGCCGAT‑CCCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGa            >  2:1794623/1‑239 (MQ=255)                                                                                                                                                                                                                                 ccGAT‑CCCCAGAACCCGCCGAAAGAAACCCCCCGCCTGGTGGAAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGAACg         <  2:1814201/238‑1 (MQ=255)                                                                                                                                                                                                                                   gAT‑CTCCAGACCCCGCAAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGAACgcag      <  1:1487149/239‑1 (MQ=255)                                                                                                                                                                                                                                   gAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGAACgca       >  1:949519/1‑238 (MQ=255)                                                                                                                                                                                                                                   gAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGAACgca       >  1:734711/1‑238 (MQ=255)                                                                                                                                                                                                                                   gAT‑CCCCAGAACCCGCAAGAAGAAACCCCCCGCCGGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGGACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGAACgcag      <  1:1864854/239‑1 (MQ=255)                                                                                                                                                                                                                                    aT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCTTATTGTCGCTGCGATGTTGCACTGCCATGAACgcag      >  2:104232/1‑238 (MQ=255)                                                                                                                                                                                                                                    aT‑CCCCAAAATCCGCCAAAAGAAACCCCCCGCCTGGGGGCAAAAGCCGATGAAGGTTACGACGGGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGAACgcag      <  1:405959/238‑1 (MQ=255)                                                                                                                                                                                                                                        cccAAAAGCCGGCAAAAGAAACCCCCCGCCGGGGGGAAAAGGCGGAGGAGTGTTACGACGTGGTAGGGCCTGTCCGCCGGAACCGCCAGGACAGCTTGTTTCGTAAAACCGCTTCGAAGATGATTATCCGGCTTATTAAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTACCGCTGTCATATTGTCGATGCGATGTTGCACTGCCATGAACGCAGCAcc  <  2:280598/238‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           GCGCACCACCACAGCCTGTGAATCGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCA‑GC‑GATTATGGC‑GGGATTCAGTCACGTTACTGGCGACTTAATTATTA‑CCCTT‑GA‑TGCCGAT‑CTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGAACGCAGCACC  >  NC_007779/2371824‑2372286

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