Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	954,395	G→A	21.9%	V24V (GTC→GTT)	focA ←	formate channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	954,395	0	G	A	21.9%	50.2 / 2.7	44	V24V (GTC→GTT)	focA	formate channel
Reads supporting (aligned to +/- strand):  ref base G (15/18);  new base A (5/5);  total (21/23)
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.93e-01

AACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCA‑GAGAGAAGCAAATGCCGCCAACCAGTTTT‑GCCATGCCGA‑AGGGCATTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑TAA‑TCGCCAG‑ATAGAAAGTCTTA‑AGCGG‑ATGTTTCG‑TTGCTTTATAGACA‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCT‑CTTTATTAAGTCGGCGACG‑AGATACTAACAAAGCAT‑TATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTA  >  NC_000913/954173‑954629                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  agaaaaaaccAACAGGGTGGAAGAAAAGAGATCGGCTCCGCAGAAAACAAAAAGAACCAGCCCAA‑GAGAGAAGAAAATGCCCCCAACCAGTTTT‑GCCATGCCGA‑AGGGCTTTGTGCCTGTGCCAGTGGT‑TGCTGTAATATAGAAGACGAATGCGATTGAGATGAAAACCCC‑GGCGG‑TAA‑TCGCCAG‑ATAGAAAGTCTTA‑AGCGG‑ATGTTTCG‑TTGTTTTATAGACA‑CCCG‑CCTCT‑TCgg                                                                                                                                                                                                                                <  1:829098/230‑1 (MQ=255)             cgggtggaaagaaaagacgccGGCCCCGCAGACAACAAAAAGAACCACCCCGA‑GAGAGAAGCAAGCGCCGCCAA‑CAGTTTTGGCCATGCCGA‑AGGTCATTGTGCCTGTGCCAGTGGT‑TGCTGTGATAAAGAAGACGAATGCGATTGAGATGAAAACCCC‑GGGGG‑TAA‑CCGCCAG‑ATAGAACGTCTAA‑AGCGG‑ATGTTTCG‑TTGCTTTATAGACC‑CCCG‑CCTCT‑TCgg                                                                                                                                                                                                                                <  2:330254/207‑1 (MQ=255)                  tGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCA‑GAGAGAAGCAAATGCCGCCAACCAGTTTT‑GCCATGCCGA‑AGGGCATTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑GAA‑TCGCCAG‑ATAGAAAGTCTTA‑AACGG‑AGGTTTCG‑TTGCTTTTTAGAAA‑CCCC‑CCCCT‑TTGG‑CCCCCTTTGccctgtct                                                                                                                                                                                                              >  2:829098/1‑233 (MQ=255)                     aaataaaggtagtcgCTCCGAACAAAACAAAAAGACTCCGCCCCA‑GAAGGAGACAAATGCCGCCAACCATTTTG‑GCCATGCAGA‑AGGGAATTGTGCGTGTGCAAGGTGT‑TGCTTTAATATAAAAGACTAAGGCGTTGGAGATAAAACCCCC‑GGCGG‑TAA‑CCGCCAG‑ATAGAAAGTCTAA‑AGGGG‑ATGTTTCG‑TTGTTTTATAGCCC‑ACCG‑CCTCT‑TCGG‑CCACTTTGGCCATTgcagca                                                                                                                                                                                                           <  2:191430/225‑1 (MQ=255)                       gTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCA‑GAGAGAAGCAAATGCCGCCAACCAGTTTT‑GCCATGCCGA‑AGGGCATTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑TAA‑TCGCCCG‑ATAGAAAAGCTTA‑AGCGG‑GTGTTTCG‑TTGCTTTATAGACA‑CCCG‑CCCCT‑TTGG‑CCCCCTTGGcccgttctcttat                                                                                                                                                                                                         >  1:330254/1‑228 (MQ=255)                       gTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCA‑GAGAGAAGCAAATGCCGCCAACCAGTTTT‑GCCATGCCGA‑AGGGCATTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑TAA‑TCACCCG‑CTAGAAAGTCTTA‑AGCGG‑ATGGTTAC‑TTGCTTTTTAAACA‑CCCG‑CCTCT‑TCGG‑CCCCCTTTGTCATTGGAGCAgg                                                                                                                                                                                                         >  1:902807/1‑239 (MQ=255)                           agagATCGGCTCCGCAGACAACACAAAGAATCAGCCCCA‑GAGAGAAGCAAATGCCGCCAACCAGTTTT‑GCCATGCCGA‑AGGGCATTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑TAA‑TCGCCAG‑ATAGAAAGTCTTA‑AACGG‑GTGTTTCG‑TTGCTTTTTAGAAA‑CCCC‑CCCCT‑TCGG‑CCACTTTTGCCCTTGCCGCAGGAAga                                                                                                                                                                                                     >  1:562052/1‑238 (MQ=255)                              gaTCGGCTCCGCAGACAACACAAAGAATCAGCCCCA‑GAGAGAAGCAAATGCCGCCAACCCGTTTT‑GCCATGCCGA‑AGGGCATTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑GAA‑TCCCCGG‑ATCGAAAGTCTTA‑AGCGG‑ATGGTTCG‑TTGCTTTATAGACA‑CACG‑CCCCT‑TCGG‑GCACTTTTGCCCTTTACGCAGGAAGTtaa                                                                                                                                                                                                  >  1:427710/1‑239 (MQ=255)                              gaTCGGCTCCGCAGAAAACAAAAAGAATC‑GCCCCAGGAAAGAAGCAAATGCCGCCAACCAGTTTG‑GCCAGGCCGA‑AGGGATTTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGGGATTGAGATAAAACCACC‑GGCGG‑TAA‑TCGCAAG‑ATAGAAAGTCTAA‑GGCGG‑ATGTTTCG‑TTGTTTTATAGACA‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTaa                                                                                                                                                                                                   <  1:88115/238‑1 (MQ=255)                               aTCGGCTCCGCAGACAACACAAAGAATCAGCCCCA‑GAGAGAAGCAAATGCCGCCAACCAGTTTT‑GCCATGCCGA‑AGGGCCTTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑GAA‑TCGCCAG‑CAAGAAAAACTTA‑AGCGG‑AGTGTTCG‑TTGGTTTATAAACA‑ACCGCCCCCT‑GCGC‑CCACTTT‑GCCCTTTCAGGAGGAAATaaaa                                                                                                                                                                                                 >  2:888775/1‑239 (MQ=255)                               aTCGGCTCCGCAGACAACACAAAGAATCAGCCCCA‑GAGAGAAGCAAATGCCGCCAACCAGTTTT‑GCCATGCCGA‑AGGGCATTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑TAA‑TCGCCAG‑ATAGAAAGTCTTA‑AACGG‑ATGTTTCG‑TTGCTTTATAGAAA‑CCCC‑CCCCT‑TCGG‑CCCCCTTTGCCCTTTCAACAGGAAGGaaaa                                                                                                                                                                                                 >  1:291429/1‑239 (MQ=255)                               aTCGGCTCCGCAGACAACACAAAGAATCAGCCCCA‑GAGAGAAGCAAATGCCGCCAACCAGTTTT‑GCCATGCCGA‑AGGGCATTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑GAA‑TCGCCAG‑ATAGAAAAGCTTA‑AACGG‑GTGTTTCC‑TTGCTTTTTTAACA‑CCCC‑CCTCT‑TTGGTCCACTTT‑GCCATTGTCGCAGGGAGTaaaa                                                                                                                                                                                                 >  1:385794/1‑239 (MQ=255)                                       cGCAGACAACACAAAGAATCAGCCCCA‑GAGAGAAGCAAATGCCGCCAACCAGTTTT‑GCCATGCCGA‑AGGGCATTGTGCCTGTGCCAGAGGT‑TGCTGTGATATAGGAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑TAA‑TAGCCCG‑ATAGAAAAGCTTA‑AACTG‑ATTGTTCG‑TTGCTTTTTAGACA‑CCCT‑CCCCTCTCGG‑‑CACATTGGCCCTTGC‑CAAGGGAGTAAAAGATCtaaaa                                                                                                                                                                                        >  1:86830/1‑238 (MQ=255)                                            acaacaCAAAGAATCAGCCCCA‑GAGAGAAGCAAATGCCGCCAACCAGTTTT‑GCCATGCCGA‑AGGGCATTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGGGATGAAAACACC‑GGCGG‑TAA‑TTGCCCG‑ATAAAAAGTCTTA‑AACGG‑GTGTTTCC‑GTGCTTTTTTGACA‑CCCC‑CCCCC‑TCGG‑CCCCTTTGGCCCATTCAGCAGGAAAGAAAAAAACAAAAagggt                                                                                                                                                                                    >  2:729803/1‑237 (MQ=255)                                              aacaCAAAGAATCAGCCCCA‑GAGAGAAGCAAATGCCGCCAACCAGTTTT‑GCCATGCCGA‑AGGGCATTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑GAA‑TCGCCAG‑ATAGAAAGTCTTA‑AACGG‑ATGTTTCG‑TTGCTTTTTAAACA‑CCCG‑CCTCT‑TCGG‑CCCCTTTGGCCATTGCACCAGGCAGTAAAAAATCAAAAGGGttt                                                                                                                                                                                   >  1:529095/1‑237 (MQ=255)                                                       aacccgCCCCA‑GGAAAAAGCAAATGCCCACAACCAGTTTT‑GCCATGCCGA‑GGGGAATTGTGCCTGTGCCAGTGGT‑TGCGGGGATAAAAAAGACGAATGCGATTGAAATGAAAACCCC‑GGCGG‑TAA‑GCGCCAG‑ATAGAAAGTCTTA‑AGCGG‑ATGTTTCG‑TTGCTTTATAGACA‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTca                                                                                                                                                                         <  1:729803/234‑1 (MQ=255)                                                          cccccccA‑GAGAGAAGAAAACTCCGACAACCATTTGT‑CCGATGCCGA‑AGTGTATTGTGCCTGTGCCGGTGTT‑TGCGGTAAAAAAGAAGACGAATGTAATGGAGATGAAAACCCC‑GGCGT‑AAA‑GCGCAAG‑ATAGAAAGTCTTA‑AGCGG‑TTGTTTTG‑TTGTTTTAAAGACA‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCacac                                                                                                                                                                      <  2:371535/236‑1 (MQ=255)                                                             ccccA‑GAGAGAAGCAAATGCCGCCAACCAGTTTT‑GCCATGCCGA‑AGGGCATTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑TAA‑TCCCCAG‑ATAGAAAGTCTTA‑AGCGG‑ATGTTTCG‑TTGCTTTATAGACA‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGGAAAAGATCAAAAAGGTTGTCAGCTTTTACCCTa                                                                                                                                                                    >  1:781482/1‑238 (MQ=255)                                                                     ggaaagaaacggccaacaactgttttt‑GCCGCGCGGA‑GGGGCTTTGGGCGGGTCCCGGTGGT‑GGCTGTAATAAAAAAGAGGAAGGCGTTGGGAATGAAAACCCC‑GGGGG‑AAA‑CCGCAAG‑A‑AAAAAATTTTA‑AGGGGTGTGTTTGT‑TTGTTTTATAAACC‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActctct                                                                                                                                                             <  2:1034226/215‑1 (MQ=255)                                                                     gaagaacaaatgccccccaaagttttt‑GGCATGCCGAGGGGGCATTGTGCGGGTG‑CAGTGGTGGGGTTAGATATAAAAGATAAAGGCGTTGGGAATAAAAACCCC‑GG‑GG‑TAACTCGCAGG‑AGAAAAATCCTAA‑AGGGT‑TTTTTTCG‑TTGTTTTATAGACC‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActc                                                                                                                                                                <  2:781482/213‑1 (MQ=255)                                                                        aaGCAAATGCCGCCAACCAGTTTT‑GCCATGCCGA‑AGGGCATTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑TAA‑TTGCCAC‑GTAGAAAAACTTA‑AGTGG‑GTGTTTCT‑GTTCTTTTTATAAC‑CCCG‑CCCCT‑TCAG‑CCACTTTGGCCATTTCTACTGGAAAGAAAAAAACAAAAAGGGTGTCaccattacaaataaaactcccttt                                                                                                                                                         >  2:185357/1‑216 (MQ=255)                                                                                  cgcaAACCGGTTTT‑GCCTGGGCGA‑AGGGATTTGTGCCTGGGCCTTGGGT‑TGGTGGGAAAAAGAAGCAGAATGCGA‑TGAGAAGGAAACCCC‑GGCGGTTAA‑TC‑CCAGAATAGAAAGTCTAA‑AGGGG‑ATTTTTCG‑TTGCTTTAAAGACA‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCT‑CTTTATTAAGTCgg                                                                                                                                              <  2:241807/239‑1 (MQ=255)                                                                                            tttt‑GCCATGCCGA‑AGGGCATTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑TAA‑TCGCCCG‑ATAGAAAGTCTTA‑AGCGG‑ATGTTTCG‑TTGCTTTATAGACA‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAAATCAAAAAGGGTGTCAGCTTTCACACTCACTCTAT‑CTTTATTAAATCGTCGACg‑agga                                                                                                                                    >  1:885825/1‑236 (MQ=255)                                                                                                 gCCATGCCGA‑AGGGCATTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑TAA‑TCGCCAG‑ATAGAAAGTCTTA‑AGCGG‑ATGTTTCG‑TTGCTTTATAGACC‑CCCG‑CCTCT‑TCGG‑CCACCTTTGCCATTGCAGCAGGAAGTAAAAGATCAAAAAGGTTTTCAGCTTTCACACTAACCCTCT‑CTTTATTAAGTCGGCGACG‑GGATACTaa                                                                                                                               >  1:167228/1‑239 (MQ=255)                                                                                                  ccATGCCGA‑AGGGCATTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑TAA‑TCGCCAG‑AAAGAAAGACTTA‑AGCGG‑AGGTTTTG‑TTGCTTTATAAACA‑CCCG‑CCCCT‑TCGG‑CCCCTTTTGCCATTGCGACAGGACGGAAAAGAACAAAAAGGTTTTCAGCTTTTACACCAAC‑CTCTCCTTTTTTaattagggggagaaaaaaaaaa                                                                                                                               >  2:572044/1‑219 (MQ=255)                                                                                                      tgcGA‑AGGGTATTGGGCCTGTGCAGGTGGT‑TGCGGTGATAAAAAAGAAGAATGCGATTGAGATGAAAACACC‑GGCGG‑TAA‑TCGCCAG‑ATAGAAAGTCTTA‑AGCGG‑ATTTTTCG‑TTGTTTTATAGACC‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTaaaa                                                                                                                                                                                                 <  1:706024/168‑1 (MQ=255)                                                                                                       ccGA‑GGGGCTTTGTGCCTGGGCCAGTGTT‑TGCTGTAAAAAAAAAGACGAATGCGATTGAAATGAAACCCCC‑GGCGG‑TAA‑CCCCAAG‑ATAAAAAGTCTTA‑AGCGG‑TTGTTTCG‑TTGCTTTATAGACA‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCT‑CTTTATTAAGTCGGCGACG‑AGATACTAACaa                                                                                                                            <  1:810344/236‑1 (MQ=255)                                                                                                        ccgaggggatttggcccggggccAGTGTT‑TGCTGGGATATAGAATACGAATGCGATTGAGATAAAAACCCC‑GGGGT‑TAA‑TCGCAAG‑AAAAAAAGTTTAA‑AGGGT‑ATTTTTCT‑TTGCTTTATAGACA‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCT‑CTTTATTAAGTCGGCGACG‑AGATACTAACaa                                                                                                                            <  1:805843/216‑1 (MQ=255)                                                                                                                      gcctgtgcaAGTGGT‑TGGTGTGATATAGAAGCCGAATGCGATGGAGATGAAACCCCC‑GGCGG‑TAC‑CCGCCAG‑ATAGAAAGTCTTA‑AGCGG‑ATGTTCCG‑TTGCTTTATAAACC‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCT‑CTTTATTAAGTCGGCGACG‑AGATACTAACAAAGCAT‑TATAGATGAGaa                                                                                                          <  1:258114/239‑1 (MQ=255)                                                                                                                        ttggcccggggtt‑tGGTGTAATAAAGAAGAAGAATGCGTTTGAGATGAAACCCCG‑GGGGG‑TAA‑CCGACAG‑ATAGAAAGTCTTA‑GGGGG‑TTGTTTCG‑TTGCTTTATAGACC‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCT‑CTTTATTAAGTCGGCGACG‑AGATACTAACAAAGCAT‑TATAGATGAGAAAt                                                                                                        <  1:700175/227‑1 (MQ=255)                                                                                                                         gtggtcggttggttgtgttaatatAAAAGACGAAGGTGTTTGAAATAAAACCACC‑GGCGG‑AAA‑TCCCAAA‑ATAGAAATTTTAA‑AGCGG‑TTTTTTCT‑GTGTTTAAAAGACC‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCT‑CTTTATTAAGTCGGCGACG‑AGATACTAACAAAGCAT‑TATAGATGAGaaa                                                                                                         <  2:562052/218‑1 (MQ=255)                                                                                                                          gtgccAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑TAA‑TCGCCAG‑ATAGAAAGTCTTA‑AGCGG‑ATGTTTCG‑TTGCTTTATAGACA‑CCCG‑CCTCT‑TCGG‑CCACTTTTGCCCTTGCAGCAGGAAGTAAAAAATCAAAAAGGGTGTCAGCTTTCACCCTAACTCCTT‑CTTTATTAAATAGGGGAAG‑AGATACTAACaaaagcatataaagaggaaaatt                                                                                                       >  1:967075/1‑219 (MQ=255)                                                                                                                          gtgccAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑TAA‑TCGCCAG‑ATAGAAAGTCTTA‑AGCGG‑ATGTTTCG‑TTGCTTTATAGACA‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCCACCTAACTCTCT‑CTTTATTAAGTCGGCGACG‑CGATACAAACCAAGAAt‑ttaaacttttcctttt                                                                                                      >  1:466646/1‑224 (MQ=255)                                                                                                                              gaggggT‑TGTGGTGATATAGAAGACGAATGCGATGGAGATAAAACCACC‑GGGGG‑TAA‑CCGCAAG‑ATAGAAAGTCTGA‑AGCGG‑ATTTTTCG‑TTTCTTTATAGACA‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCT‑CTTTATTAAGTCGGCGACG‑AGATACTAACAAAGCAT‑TATAGATGAGAAATTGatat                                                                                                  <  1:888775/235‑1 (MQ=255)                                                                                                                                ggtgtgtgctggaagagagaagtCAAATGCGAGAGAGAGGAAACCACC‑GGCGG‑TAA‑TCGACAG‑ATAAAAAGTCTAA‑AGCGG‑TTGTGGAG‑GTTGTTTAAAACCC‑CACG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCT‑CTTTATTAAGTCGGCGACG‑AGATACTAACAAAGCAT‑TATAGAGGAGAAATTGatat                                                                                                  <  1:842093/221‑1 (MQ=255)                                                                                                                                gTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑TAA‑TTGCCAG‑ATAGAAAATCTTT‑AGAGG‑ATGTTTCG‑TTGCTTTTTAGACA‑CCCG‑CCTCT‑TCGG‑CCAACTTGGCCCTTGCTGCCGGAAGGAAACAATCCAAAAGGGTGTCCGCTTTTACACCACACCTCC‑TTTTATTTAAGCGGCGGCGAAGACAC‑AACAAAACATATAAACATAAAAAATTGataaa                                                                                                 >  2:915126/1‑237 (MQ=255)                                                                                                                                   gttgttgtatataaAGAAGAGTAAGGCATTTGAAATAAAAACCCC‑GGAGG‑AAA‑CCGCGAA‑GAAGAAA‑TGTTATAGGGG‑ATGTTTTG‑TTGGTTAATAATCA‑CCAG‑CCTCT‑TCGG‑CCTATTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCT‑CTTTATTAAGTCGGCGACG‑AGATACTAACAAAGCAT‑TATAGATGAGAAATTGATATAGAt                                                                                              <  2:385794/229‑1 (MQ=255)                                                                                                                                                               aTTGAGATGAAAACACC‑GGCGG‑TAA‑TCGCCCG‑ATAGAAAGGCTTA‑AGCGG‑ATGTTTCG‑TTGCTTTATAGACA‑CCCG‑CCTCT‑TCGG‑CCACTTTGGc                                                                                                                                                                                                                     >  1:617411/1‑95 (MQ=255)                                                                                                                                                                      taaaaccccc‑GGGGG‑TAA‑CCGCAAA‑AAAAAAAGCTTTA‑AGCGG‑‑TTTCTCGTTTGTTTAAAAGCCA‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCT‑CTTTATTAAGTCGGCGACG‑AGATACTAACAAAGCAT‑TATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCa                                                           <  1:105377/231‑1 (MQ=255)                                                                                                                                                                       taaaaCCCCGGGCGG‑AAA‑TCGCAAG‑AAAGAAATTCTTA‑AGCGG‑TGTTTTCG‑TTGCTTAATAGCCC‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCT‑CTTTATTAAGTCGGCGACG‑AGATACTAACAAAGCAT‑TATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCa                                                           <  1:572044/238‑1 (MQ=255)                                                                                                                                                                        aaaaCGCG‑GGGGG‑AAA‑TCGCAAG‑ATAGAAA‑TTTTA‑AGGCG‑GTTTTTCGTTTTTTTTAAAGCCC‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCT‑CTTTATTAAGTCGGCGACG‑AGATACTAACAAAGCAT‑TATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCa                                                           <  1:130892/237‑1 (MQ=255)                                                                                                                                                                                         a‑TCGCCAG‑AGAAAAAGTTTTA‑AGCGG‑TTGTTTCG‑TTGTTTTA‑AAACACCCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCT‑CTTTATTAAGTCGGCGACG‑AGATACTAACAAAGCAT‑TATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATa                                          <  1:185357/239‑1 (MQ=255)                                                                                                                                                                                           accccaaattaaaaaTTCTTA‑GGCGG‑TTTTTTCG‑TTGCTTTAAAGCCC‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCT‑CTTTATTAAGTCGGCGACG‑AGATACTAACAAAGCAT‑TATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATa                                          <  1:247895/227‑1 (MQ=255)                                                                                                                                                                                                                                      tataGACA‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCT‑CTTTATTAAGTCGGCGACG‑AGATACTAACAAAGCAT‑TATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAAAATTTGATTATCCAAATAAAAATAAATTTTa  >  2:690961/1‑239 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  AACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCA‑GAGAGAAGCAAATGCCGCCAACCAGTTTT‑GCCATGCCGA‑AGGGCATTGTGCCTGTGCCAGTGGT‑TGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACC‑GGCGG‑TAA‑TCGCCAG‑ATAGAAAGTCTTA‑AGCGG‑ATGTTTCG‑TTGCTTTATAGACA‑CCCG‑CCTCT‑TCGG‑CCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCT‑CTTTATTAAGTCGGCGACG‑AGATACTAACAAAGCAT‑TATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATAAATTTTA  >  NC_000913/954173‑954629

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