Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	817,815	G→C	14.7%	A258P (GCC→CCC)	moaA →	molybdopterin biosynthesis protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	817,815	0	G	C	14.7%	113.1 / 3.5	54	A258P (GCC→CCC)	moaA	molybdopterin biosynthesis protein A
Reads supporting (aligned to +/- strand):  ref base G (36/9);  new base C (2/7);  total (38/16)
Fisher's exact test for biased strand distribution p-value = 1.61e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCATCGCCCTATCCAGCTGCGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTC‑TGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGAGAAGAAACAGACCCATTTCCTGCATCAAAACAACACCGGTATTACGCAAAACTTATCGTACATTGGCGGCTAAAACGTC  >  NC_000913/817588‑818039                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   cgcaccccctcttCCAGCTGGCTTTCATCGAACTGATGGAATCGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCCCATCTC‑TGGTCAGGTTCTGCGGGACGAGCTACTGCGTCGCTGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTgcg                                                                                                                                                                                                                                    <  2:360633/216‑1 (MQ=255) aGCATCCCCCTTTCCCGCTGCGTTGTACCGACCTGATGGAAACGGGGGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTC‑TGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTgcg                                                                                                                                                                                                                                    <  1:509825/228‑1 (MQ=255)        cccTATCCAGCTGCGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTC‑TGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGCCGGT‑CCCGCGCAAG‑GCTTTTGCCATCCAG‑C‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAACAACTTCTGCGCCCCttt                                                                                                                                                                                                                             >  1:540243/1‑227 (MQ=255)             tCCAGCTGCGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTC‑TGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑GCTTTTGCCATCAAG‑A‑TTACACC‑GGAGCGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCAGCCCCACATGCA‑Acc                                                                                                                                                                                                                       >  1:354296/1‑228 (MQ=255)                  cTGCGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTC‑TGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑TCTTTTTCCATCCAG‑A‑TTAAGCC‑GGAGAGATTGGCCTTCTCA‑T‑GCCGTATGAAAAAGATTTCTGCGCCACTCGCA‑ACCGCCt                                                                                                                                                                                                                   >  1:468485/1‑227 (MQ=255)                  cTGCGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTC‑TGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGCAAAAGACTTCTGCGCCACTTGCA‑ACCGCCt                                                                                                                                                                                                                   >  2:465635/1‑227 (MQ=255)                  cTGCGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTC‑TGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCACCGACGGT‑CCCGCGCAAG‑GCTTTTGCCCTCCAG‑A‑TTACCCC‑CGGGGAATGGGCGTTATCA‑T‑GCCGTATGAAAAAACATTTTGCGCCACTTGCA‑ACccccca                                                                                                                                                                                                                  >  1:489110/1‑226 (MQ=255)                  cTGCGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTC‑AGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCCCCAATTACGTCAACGCAGCGACAGT‑CCCGCGCAAG‑TCTTTTGCCCTCCAG‑C‑TTACGCC‑GGAGAGAGTGGCCTTATCA‑T‑GCCGCATGAACAAGCATTCGCCGGCACTTGCA‑ACCACCTg                                                                                                                                                                                                                  >  1:165106/1‑228 (MQ=255)                     cGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTC‑TGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCCATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAGGACTTTTGCCCAACTTGCA‑ACCCCCTgccc                                                                                                                                                                                                               >  2:536373/1‑226 (MQ=255)                     cGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTC‑TGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GCAGAGATTGGCCTTACCA‑T‑CCCGTATGAAAAAAACTTCTCCGCCACTTGCA‑ACccccacccc                                                                                                                                                                                                               >  1:43575/1‑223 (MQ=255)                               aaCTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTC‑TGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑ACTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCa                                                                                                                                                                                                     >  2:189806/1‑228 (MQ=255)                                                  gggCAGCGAGCTCTTCCGTAAGCATCACATCTC‑TGGTCAGGTTCTGCGTGACGAGATACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑GCCTTTGCCATCCAC‑A‑TTACGCC‑GGAGAGATTGGCCTTGTCA‑T‑GCCGTATGAAAAAGAATTCTGCGCCCCTTGCA‑ACCGCCTCCGCGTTTCCACAAGTAGCAAACCCCAActctc                                                                                                                                                                                  >  1:414753/1‑227 (MQ=255)                                                        gggctttttCCGTAAGCATCAC‑CCTCTTGGCGAGGTTCGGGGGGGCGGGGTAGTGCGTGGGGGCTG‑GATCCCCCAATTACGTCAACGCGGCGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTCGGCC‑GGAAAGATGGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGAA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGt                                                                                                                                                                             <  1:512316/220‑1 (MQ=255)                                                                     aaGCATCACATCTC‑TGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAggcgg                                                                                                                                                                >  2:554865/1‑227 (MQ=255)                                                                                      gTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGCAAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGAGAAGACGGCGTTAAACTGCGCGATc                                                                                                                                               >  2:476883/1‑228 (MQ=255)                                                                                      gTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATc                                                                                                                                               >  2:369730/1‑228 (MQ=255)                                                                                      gTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑GCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCCTCTCTCCCCGGTTGGTGAAGGCGGCGTTAAcacccgcaacc                                                                                                                                               >  1:394541/1‑218 (MQ=255)                                                                                           gctCTGCTTGACGAGCTACTGCGTCGCGGCTG‑GATTCACCAATTACGTCACCGCAGCGACGGC‑CCCGCGCAAG‑TCTCTTGCCATCCAG‑A‑TTACGCC‑GGCGAGATTGGCCTTATCA‑T‑CCCGCATGCAACAGACTTCTACGCCACTCGCC‑ACCGCCAGCGAGCGTCCTCCCTCGGTAAACTCCCACCCTGCCTGCTGGGTGCAGCCGGCGTCAACCTgcaccacctcccg                                                                                                                                          >  1:695554/3‑217 (MQ=255)                                                                                           gTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGCAGCCGGCGTTAACCTGCGCGATctgct                                                                                                                                           >  2:400053/1‑227 (MQ=255)                                                                                           gTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GCAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCGCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGAGCCGTTAACCGGCCGGATctcct                                                                                                                                           >  2:191651/1‑227 (MQ=255)                                                                                            ttCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCCCTGCCGGTTTGGTGAAGCCGCCGTTAACCTGCGCGCTctgctg                                                                                                                                          >  2:415723/1‑227 (MQ=255)                                                                                              cTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGGGAAGGCGGCGTTAACCTGCGCGATCTCCTGGaa                                                                                                                                       >  1:351489/1‑228 (MQ=255)                                                                                                gCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATCAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGGGACGGCGGCGTTAACCCCCCCGATctctcgaacga                                                                                                                                     >  1:227361/1‑219 (MQ=255)                                                                                                        cgcctcccgggtccgGCTGTGATCCCCCAATGACACAAACGCAGCG‑CGGC‑GCCGGCAAAGTTCTGTTGCCA‑CCAGAA‑TTACGCC‑GTATAGATTGCCCTTATCA‑T‑CCCGTATGAAAAAGCCTTCTGCCCCATTTGCA‑CCCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATAcc                                                                                                                              <  1:476883/215‑1 (MQ=255)                                                                                                              cgggTCGGGGGTG‑GAGCCCCCAATTTCCTAAAAGCGGAGGCGGTCCCCGCGAAAG‑TTTTTTGCCCTCAAG‑T‑TTCCGCC‑GGAGAGATGGGCTTTTCCA‑T‑GCCGTATAAAAAAGCCTTCTGCGCCACTGGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCAGCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCaa                                                                                                                        <  1:554865/225‑1 (MQ=255)                                                                                                                  tCGGGGCTG‑GATCCACCAATTCCGTCAACGCAGCGCCGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTCCGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑gg                                                                                                                   <  2:118598/227‑1 (MQ=255)                                                                                                                  ggcggggctctaTCCCCAAATTCCTTCACCGCAGCGACGGT‑CCCGCACAAG‑TCTTTTGCCACACTG‑T‑TTCCGCC‑GGGGAGATGGGCCTTATCA‑T‑GCCGTATGAAAAAGCCTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑gg                                                                                                                   <  1:400053/217‑1 (MQ=255)                                                                                                                          g‑gATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGTTGAAGGCGCCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAg                                                                                                          >  2:303400/1‑228 (MQ=255)                                                                                                                                          gTCAACGCAGCGACGGT‑CCCGCGCATG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCACGTTTCCTCCATTGGTAAACTCCATACCTGCCTGTTTGGTTCAGTCGGCGTTAACCTGCTCTATCTCACGTAAACCGATCCACAGCAACA‑CGCCCTgcgcaggcggatttagacgg                                                                                           >  1:560345/1‑209 (MQ=255)                                                                                                                                                  caggcgcggtgcacggagaG‑TTTTGTGCCATCCAG‑A‑TT‑CGCCTGGAGAG‑TTGGGCTTTTCC‑TGGCCGTATGAAAAAGGCTTCTGCCCCTCTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCg                                                                                    <  2:165106/210‑1 (MQ=255)                                                                                                                                                    cGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCCGGAATTTCACCGGCGCTGCggg                                                                                  >  2:73480/1‑227 (MQ=255)                                                                                                                                                    cGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGACCTGCTGGAAGACGATACCCAGCAACA‑GCGGCTGGAGCCGCGAATTTCAGCGGCGCTGCggg                                                                                  >  1:234547/1‑227 (MQ=255)                                                                                                                                                         gT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGagatga                                                                            >  2:267289/1‑228 (MQ=255)                                                                                                                                                                ggcaaattttttgGCCATCAGA‑A‑TTCCGCC‑GGAAAGTTTGGCCTTACAA‑G‑GCGGTAGGAAAAAGCCTTCTGGCCCCCTGGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGAGAAGAAACa                                                                        <  2:354296/218‑1 (MQ=255)                                                                                                                                                                 gcAAG‑TTTTTGGCCATCCAG‑A‑TTCCGCC‑GGAGAGATTGGCCTTATCA‑G‑GCCGTATGAAAAAGCTTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGAGAAGAAACAGa                                                                      <  2:43575/227‑1 (MQ=255)                                                                                                                                                                   gtaagtttttGGCAACCAG‑AGTTCCGCC‑GGAAAGTTTTGCCTTTCAA‑T‑GCCGTATGAAAAAGCCTTCGGCCCCCCTTCCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGAGAAGAAACAGAc                                                                     <  2:361573/222‑1 (MQ=255)                                                                                                                                                                     g‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGGGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATCCCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTCCGGGAGAAGAAACaaaacaat                                                                 >  2:699891/1‑221 (MQ=255)                                                                                                                                                                               cacccAG‑T‑TTAGGCC‑GGAGTGGTTGGCCTTCAGC‑T‑GCCTTAGAAAAAAGTATT‑TGCGCGCCTTGGACACCGGCGGCGCGTTTCCTTCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGGTCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGAGAAGAAACAGACCCATTTCCTGCAt                                                        <  1:536373/225‑1 (MQ=255)                                                                                                                                                                                   cAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATATCTGCCTGTTTGGTGAGGGCGGCGTTAACCTGCGCGCTCTGCTGGAAGACGTTACCCAG‑AACACGGAGCTGGAAGCGCGTGTTTCAACGGCGCTGCGGGCGAAGAAACAGAACCATTTCCTGCATCaca                                                    >  2:688594/1‑226 (MQ=255)                                                                                                                                                                                             gcccgcgaagttggcccTTCTCC‑T‑GCCTTAAAAAAACGACTGCCGCCCCCCTGGCA‑GCCGCCTGCGCAGTTCCTCCATTGGTAAACTCCAGCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGAGAAGAAACAGACCCATTTCCTGCATCAAAACAacac                                             <  1:191651/213‑1 (MQ=255)                                                                                                                                                                                                  gagagaTTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGCAGCGCGTATTTCAGCAGCGCTGCGGGAGCAGAAAAAGACCCAATTCCTGAATCAAAACAACACCggga                                        >  1:554222/1‑226 (MQ=255)                                                                                                                                                                                                   gagggatGGCCCTATCAGT‑GCCTTATCAAAAAGCCTCCGCCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGAGAAGAAACAGACCCATTTCCTGCATCAAAACAACACCGGTa                                        <  2:227361/222‑1 (MQ=255)                                                                                                                                                                                                   agagaTTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGa                                                                                                            >  2:76149/1‑159 (MQ=255)                                                                                                                                                                                                   agagaTTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GCCGCTGGAAGCGCGGATTTCAGCGGCGCTGCGGGCGAAGAAACAGACCCATTTCCTTCATCAAAACAACACCggact                                       >  2:188487/1‑225 (MQ=255)                                                                                                                                                                                                          ggCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACAGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCCGGGGGAGAAGAAACAGACACATTTCCCGCATCAAAACAACACCGGCATTACGCa                                 >  2:501815/1‑227 (MQ=255)                                                                                                                                                                                                          ggCCTTACCA‑T‑GCCGTATGAAAAAGATTTCGGCCCCCATAGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGAGAAGAAACAGACCCATTTCCTGCATCAAAACAACACCGGTATTACGCaa                                <  2:489110/228‑1 (MQ=255)                                                                                                                                                                                                            gctttaccct‑GCCGTATGAAAAAGCTTTCTGCCCCCCTGGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGAGAAGAAACAGACCCATTTCCTGCATCAAAACAACACCGGTATTACGCaaa                               <  2:540243/219‑1 (MQ=255)                                                                                                                                                                                                              ttcacA‑T‑GCCGTATGAAAAAGCCTTCTGCGCCGTTGGCA‑CCCGGCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGAGAAGAAACAGACCCATTTCCTGCATCAAAACAACACCGGTATTACGCAAAAc                             <  2:394541/223‑1 (MQ=255)                                                                                                                                                                                                                        ccGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAACCTCCATCTCTGCCTGTTTGGGGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGCGCCGAAACGGACCCATTTCCTCACTCCAAACAGAACCGGTATTACGCacgtcccttatac                     >  1:40037/1‑215 (MQ=255)                                                                                                                                                                                                                        ccGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGAGAAGAAACAGACCCATTTCCTGCATCaaaa                                                   >  2:452252/1‑197 (MQ=255)                                                                                                                                                                                                                        ccGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGAGAAGAAAAAGACCCATTTCCTCAATCAAAACAACACCGGTATTAGGCAAAAATAATCGTaa                    >  1:194734/1‑227 (MQ=255)                                                                                                                                                                                                                            aTGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTaa                                                                                                                                                          >  2:470793/1‑91 (MQ=255)                                                                                                                                                                                                                             tGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAAc                                                                                                                                                         <  1:470793/91‑1 (MQ=255)                                                                                                                                                                                                                             tGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGAGACGAAACGAACCCATTTCCTGAATCAAAACAACACCGGATTTACGCAAAACCTAACGGACAATgg               >  1:375502/1‑228 (MQ=255)                                                                                                                                                                                                                                     cctttcgcgcCACGTGCA‑ACCGCCTGCGCGGTTCCTCCATTGCTAAACTCCAACTCTGCCTGTTTGGGGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATCCCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGAGAAGAAACAGACCCATTTCCTGCATCAAAACAACACCGATATTACGCAAAACTTATCGTACATTGTCGGCTaaa       <  2:560345/222‑1 (MQ=255)                                                                                                                                                                                                                                     aCTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTTCGGGCGAAGAAACAGCCCCATTTCCGGCATCCAAACACAACCGGAATTACGCAAAACTCATCGAAAATTGGCGGCCaaa       >  1:278270/1‑228 (MQ=255)                                                                                                                                                                                                                                          tGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGAGAAGAAACAGACCCATTTCTTGCATCAAAACAACACCGGTATTACGCAAAACTTATCGTACATTGGCGGCCAAAAcacc  >  2:651758/1‑225 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   AGCATCGCCCTATCCAGCTGCGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTC‑TGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTG‑GATCCACCAATTACGTCAACGCAGCGACGGT‑CCCGCGCAAG‑TCTTTTGCCATCCAG‑A‑TTACGCC‑GGAGAGATTGGCCTTATCA‑T‑GCCGTATGAAAAAGACTTCTGCGCCACTTGCA‑ACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTTAACCTGCGCGATCTGCTGGAAGACGATACCCAGCAACA‑GGCGCTGGAAGCGCGTATTTCAGCGGCGCTGCGGGAGAAGAAACAGACCCATTTCCTGCATCAAAACAACACCGGTATTACGCAAAACTTATCGTACATTGGCGGCTAAAACGTC  >  NC_000913/817588‑818039

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