Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	954,388	C→T	28.8%	A27T (GCA→ACA)	focA ←	formate channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	954,388	0	C	T	28.8%	51.5 / 10.5	52	A27T (GCA→ACA)	focA	formate channel
Reads supporting (aligned to +/- strand):  ref base C (23/13);  new base T (2/13);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.62e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

GATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑G‑CGGATGTTT‑CGTTGCTTTATAGA‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACG  >  NC_000913/954152‑954591                                                                                                                                                                                                                                                        |                                                                                                                                                                                                               gaagcccccccctgacatAGCAACAACCATCAACCCGGGGGAAGTAAAGGAAACGGCTCCGGCAACAACAC‑AAAGAAGCAGCCCCAGCGTGAAGCAAATGCCG‑CCAACCAGGTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑G‑CGGATGTTT‑CGTTGCt                                                                                                                                                                                                               <  2:387548/221‑1 (MQ=255)   tGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATTGCC‑CGATAGAAAAGCTTAACG‑GGAATGTTT‑‑GTTGCtttt                                                                                                                                                                                                            >  1:1135772/1‑238 (MQ=255)   tGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATAGTC‑AGCTAGAAAAGCTTAA‑A‑CGGGAGTTT‑CGTGGCTTta                                                                                                                                                                                                            >  1:490041/1‑239 (MQ=255)      cccccCTCCCCTTTACAACCAACATCAACCCGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACAC‑AAAGACACAGCCCCAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑G‑CGGATGTTT‑CGTTGCTTTATAg                                                                                                                                                                                                         <  2:1094599/238‑1 (MQ=255)                    gCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TATTCGCC‑AGGTAGAAAAGCTTTA‑G‑CGGATGTTT‑CTTTGCTTTATAGA‑CACCCCCCTCTT‑c                                                                                                                                                                                         >  1:937656/1‑239 (MQ=255)                                       tGGAAGTAAAGAGATCGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAAACGCC‑ACATAGAAAGTCTTAA‑A‑CGGATGTTT‑CGTTGCTTTTTAGA‑CACCCCCCCCTT‑CGGCACCTTTGTCCCTTcca                                                                                                                                                                      >  1:818736/1‑239 (MQ=255)                                           ataaaaGAGCCCGGCCCCCCACACCACAC‑AAAAAAACACCCCAAGAGAGAAGCAAAGGCCC‑CCAACCATCTTTG‑CCAGGACGAAGGTCTTTGT‑GCCGGTGCCAGTGGTTGCTTTGATAGAGAAGACG‑AATGGTATTGAG‑ATAAAAACCCGGGGGG‑TAATCGCC‑AGAAAGAAAGTCTTAA‑G‑GGGATGTTT‑CGTTGCTTTATAGA‑CCCCCGCCTCTT‑CGGCCACTTTGGCCATTgcagca                                                                                                                                                                   <  2:1134153/235‑1 (MQ=255)                                           aaaaaaggagtcggccccaaaaacaacaA‑AAAGAACCCGCCCAAGGGAAAACAAAATGCCC‑CCAACCATTTTGC‑CCAGGGCGAAGGTTTTTGG‑CCCTGGGCCATGGGTTGCTGATATAAAGAAGAAG‑AATGCGATTGAG‑ATGAAACCCCCGGCGG‑TAACCGCC‑AGAAAGAAAGTCTTAA‑G‑GGGATGTTT‑CGTTGTTTTAAAGA‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGcagcag                                                                                                                                                                  <  2:46164/217‑1 (MQ=255)                                            agacaagagttggCTCCTCAGAAAAAATAAAAGAATCAGCCCGTGAGAAAACCAAATGCCG‑CCAACTATTTTGC‑CCAGGCCGAAGGGCGTGGG‑CCATGTCGATGTGGTGGCTGGTATAAAGAAGAGG‑AGTGGAATTTAA‑ATGAAATCCCCGGGGG‑TACTCCCA‑AGATAAAAATTCGTAT‑G‑CGGAGTTTT‑CGTTGTTTTAGAAA‑CCCCCGCCTCTT‑CGGCCACTTTGGCCATTgcagca                                                                                                                                                                   <  2:185993/227‑1 (MQ=255)                                               aaGAGATCGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAAAACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAC‑G‑CGGATGGTT‑CGTTGTCTTCTAGA‑CACACGCCTCTT‑CGGCCACCTTTGCCAtttctcccggtag                                                                                                                                                              >  1:548219/1‑228 (MQ=255)                                                    aTCGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGCG‑ATGAAAACACCGGCGG‑CAAATGCC‑CCGTAGAAAATCTTAA‑A‑CGGAAGTTT‑AGTTGTTTTATAAA‑AACCCGCCCCTT‑CGTCCCCTTTTGACACTtcagcaagggaaaaaa                                                                                                                                                         >  2:1028746/1‑229 (MQ=255)                                                    aTCGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGAC‑AGATAGAAAGGCTTAA‑A‑CGGGTGTTT‑CGGTGCTTTTTAGA‑AACCCCCCCCTT‑CCGCCCCTTTGGCCCtttccaaaggaaagaaaa                                                                                                                                                         >  2:982834/1‑223 (MQ=255)                                                      cGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAAGCTTCA‑G‑CGGCTGGTT‑CGTTGCTTTATAGA‑CACCCGCCCCTT‑CGGCCACTTTTGCCCTTTCAGCAGGAAGTaaaaaa                                                                                                                                                       >  1:657322/1‑237 (MQ=255)                                                         cTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATAGCC‑AGCTTGAAAGGCTTAA‑G‑ACGGTGTTT‑CGTTGCTTTTTAGA‑CACCCGCCTCTT‑CGGCCCCTTTTGGCATTGCAGCAGGCAGTAAAAGAtaa                                                                                                                                                    >  1:35702/1‑237 (MQ=255)                                                            cGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑GAATCGCC‑AGATAGAAAGGCTTAA‑G‑CGGATGTTT‑CGTTGCTTTATAGA‑CACCCGCCCCTT‑CGGCCCCTTTGGCCATTGCAGCAGGAAGTCAAAGATAaaaa                                                                                                                                                 >  1:269483/1‑239 (MQ=255)                                                                   aacaC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑G‑CGGATGGTT‑CGTTGCTTTATAGA‑AACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGGAAAAAATCAAAAGGGGTg                                                                                                                                           >  1:234073/1‑238 (MQ=255)                                                                   aacaC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGGCTTAA‑A‑TGGAGGGTT‑CGTTGCTTTATAGG‑CACCCGCCCCTT‑TGGTCACCTTGGCCATTGCAGCAGGAAGTAAAAGACGAAAAGGGGTgc                                                                                                                                          >  1:692570/1‑238 (MQ=255)                                                                   aacaC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑GAAACGCC‑AGATAGAAAGACTTAA‑A‑CGGGAGGTT‑TGTTGATTTTTTGA‑ACCCCCCCCCTTAGGGCCACTTT‑GCCATTTCGGCAGGGAGTaaaaaacaaaaaagggagt                                                                                                                                          >  2:351088/1‑224 (MQ=255)                                                                     cacaaaaaaaaccgCCCAAGAAAGAA‑CAAAGGCCCCCCAACCTGTTTTG‑CCATGGCGAGGGGATTTGG‑GCGGGTGCCGGTGTTGTGTGTGATATA‑AAAACGTAGTGCGATTGAA‑ATGAACACCCCGGCGG‑TAATCCGC‑AGAAAGAAGTTCTTAG‑G‑CGGATTTTT‑CGTGTTTTTAAAAA‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTc                                                                                                                                         <  2:1048863/226‑1 (MQ=255)                                                                     acaaaaaaatcgcccccaAGAGAGAACAAAAGCCCC‑CCAACCTGTTTTC‑CCAGGCCAAAGGGCTTTGG‑CTCTGGCCCTGTGGTACTGGTTATAAAGAAGACG‑AATGCTATAGAA‑ATGAAACCCCGGCCGG‑AAACCACC‑AAAAAAAATTTCTAAA‑G‑CGATTTTTT‑CGTTGTTTTATAAC‑CCCCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTc                                                                                                                                         <  2:251614/225‑1 (MQ=255)                                                                               tcgccccaAGAGGAAAGCAAGTGCCG‑CCAACTGTTTTTG‑CCAGGCCGAGGGGTTTTGG‑CCCTGTGCCAGGGGTTGCTGTAAAAAAGAAGACG‑AATGCGA‑TGAGAATAAAACCCCCGGCGG‑TAACCGCC‑AAATAGAACGTCTTAA‑G‑CGGATGTTT‑CGTTGTTTTATAGA‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTcaca                                                                                                                               <  2:269483/235‑1 (MQ=255)                                                                                cAGCCCCAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATGGCC‑AGATAGAAAGTCTTAA‑G‑CGGAAGTTT‑CGTTGCTTTATAGA‑AACCCCCCCCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAAATCAAAAGGGTTGTAACCTTTCacac                                                                                                                              >  1:796566/1‑239 (MQ=255)                                                                                   ccccAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGGTAGAAAGTCTTAA‑A‑CGGATGTTT‑CGTTGCTTTATAGA‑AACCCGCCCCTT‑CGGCCACCTTTGCCCTTGCAGCCGGAAGGAAAAAATTCAAAGGGGTTTCAGCCTTTCCACTaa                                                                                                                           >  1:1268126/1‑239 (MQ=255)                                                                                             aaGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAAGCTTAA‑G‑CGGATGTTT‑CGTTGCTTTTTAAA‑TCCCCCCCCCTT‑CGGCCCCTTTTGCCATTGCCGCCGGAAGAAAAAAATCAAAAAGGTGGTCACCCTTTCCCCTAAAT‑CCCTCttt                                                                                                                >  2:776548/1‑239 (MQ=255)                                                                                             aaGAAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGCAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGCTATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TCACTACC‑CAATAGACAAATTTTA‑G‑CGAGAGTTT‑CCTTGCATTATTCG‑CACCACGCTCTT‑TGTCCACTTTTGCCCTTTCAACCAGAGATAAAAAATTCAAAGGGTTATCAGCCTTCAACAAAACT‑CTCCCttt                                                                                                                >  2:290818/1‑239 (MQ=255)                                                                                                   aTGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑G‑CGGATGTTT‑CGTTGCTTTATAGA‑CACCCGCCTCTT‑CGGCCCCTTTGGCCCTTTCAGCCGGAAGTAAAAAATTAAAAGGGGTGTTAGCTTTTACCCTAACT‑CTCCCTTTATTAAg                                                                                                          >  1:612697/1‑239 (MQ=255)                                                                                                            aacagtttttG‑CCATG‑CGGAGGGGATTGTGGCCTGTGCCGTTGGTGGCTGTGATATAGAAGCCG‑AATGCGATTGAA‑ATGAAAACACGGGCGG‑TAACCGCC‑AGATAGAAAGTTTTAA‑G‑CGGATGTTT‑CGTTGCTTTATAGC‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGcgac                                                                                                  <  1:477529/233‑1 (MQ=255)                                                                                                             acatgTTTTGCCCATGCGAAAGGGTTGTGT‑CCTTGGGACAGTGTTTGCTGAGATAAAGAAGAAA‑ATTGGGATGAAA‑ATCAAAACCCCGGGGT‑TAATCCCA‑AGATAGAAACTTTTAA‑G‑GGGTTGTTT‑CGTTGCTTTATAGA‑CCCCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGcgac                                                                                                  <  2:1268126/235‑1 (MQ=255)                                                                                                                            gCCAAAGGGCATTGT‑GCCTGTGCGAGTGGTTGCTGTGAAATAGAAGACA‑AATGCGATTGAG‑ATGAAACCCCGGGCGG‑TACTCGCC‑AGATAGAATGTCTAAA‑G‑CGGATGTTT‑CGTTGCTTTATAGA‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAg                                                                                   <  1:351088/239‑1 (MQ=255)                                                                                                                                  gggCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAAACTTAA‑G‑CGGGTGTTT‑CGTTGCTTTTTAGA‑CACCCGCCCCTT‑CGGCCCCTTTTGCCATTGCAGCAGGAAGTAAAAAATCAAAAGGGGTGTCAGCTTTTCCCCTAACT‑CCCTCTTTTTTTAATGGGcgcaggggaaacaaaaaaacattat                                                                             >  2:1073578/1‑216 (MQ=255)                                                                                                                                      aTTGG‑GCCGGCGCAAGGGGTTTCTGTGATAAAGAAGGCG‑AAGGGGATTGAA‑ATGAAAACACGGGGGT‑TAATCGCA‑AGATAGAAAGTCTTAG‑G‑GGGATGTTT‑CGTTGTTTTATAGA‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                                                                      <  1:334579/226‑1 (MQ=255)                                                                                                                                         tgggcctggaccAGTGGTTGTTGTAATATAGAAGCAG‑AAGGCGATTGAA‑ATGAAACCCCCGGCGG‑TAACCGCA‑AGAAAGAAAGTCTTAA‑G‑GGGTTGTTT‑CGTTGCTTTATAGA‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATga                                                                       <  1:776548/236‑1 (MQ=255)                                                                                                                                              tgtggctggtggttgtggttATAAAGAAGACG‑AATGAGATAGAA‑ATAAAAGCCGCGGAGG‑TAAAGGCG‑AAGAAAAAAGTCTAAA‑G‑CGGTTGTTT‑TGTGGTTTTAACAA‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGGCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGaaa                                                                   <  2:548219/230‑1 (MQ=255)                                                                                                                                               ggggccggggtgtgggggaaaataGAAGACG‑AATGGGATTGGA‑ATGAAAACCCCGGGGG‑TAACGGCA‑AGAAAAAAAGTCTTAG‑G‑CGTATTTTT‑CGTTGTTTAATAGA‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAAt                                                                  <  1:1073578/217‑1 (MQ=255)                                                                                                                                                gtgccAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑CGATAGAAAGTCTTAA‑G‑CGGATGTTT‑CGTTGCTTTATAGA‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGa                                                                                                                                                       >  1:121102/1‑153 (MQ=255)                                                                                                                                                 tgccAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCCCC‑CGATAGAAAGTCTTAA‑G‑CGGATGTTT‑TGTTGCTTTATAAA‑AACCCGCCTCTT‑CGCCCACTTTGGCCATTGCCCCAGGAAGTAAAAGATCAAAAGGGGTGTCAGCTTTTCCACTAACTCCTCTCTTT‑TTAAATAGGCCACCCGATACTaaaaaaacctataacactctccttaat                                                               >  2:334579/1‑218 (MQ=255)                                                                                                                                                   atgttaataagagacagaaagaaaaCG‑AAGGGAATGGAA‑ATGAAACCCCGGGCGG‑AAACCCCA‑GGATAAAAAGT‑TTAA‑GTTGGATGTTC‑CGTTTCTTTAAAAA‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGGCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGat                                                              <  1:1289281/219‑1 (MQ=255)                                                                                                                                                    aaggggTTGCTTGAATATAAAAGACG‑AATGCGATTGAG‑ATAAAAACCCCGGGGG‑TAATCGCA‑AGATAGAAAGTCTTAA‑G‑CGGTTTTTT‑TGTTCCTTTA‑AGACCACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGatat                                                            <  1:701739/235‑1 (MQ=255)                                                                                                                                                        ggTTGCTGTAAAATAGAAGAAG‑AATGGGTTTGAG‑AT‑AAAACCCCGGGGGTTAACCCCC‑AGATAGAAAGTTTTAA‑G‑CGGATGTTT‑CGTTGCTTTATAGA‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGAt                                                        <  1:542082/239‑1 (MQ=255)                                                                                                                                                           tGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑G‑CGGATGTTT‑CGTTGCTTTATAGA‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGCAGTAAAAGATCAAAAAGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCat                                                     >  1:844605/1‑239 (MQ=255)                                                                                                                                                           tGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGGCTTAA‑G‑CGGATGTTT‑CGTTGCTTTATAGA‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGGAATACTAACAAAGCATTATAGCTgtctcttatacaaatttccc                                                     >  1:334968/1‑220 (MQ=255)                                                                                                                                                                     agaaagcag‑AATGCGAGTGAG‑ATAAACACCCGGGGGG‑TAATCGCA‑AGATAGAAAGTCTTAA‑G‑GGGTTGTTG‑CGTGGTTTTATAGA‑CCCCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTg                                           <  1:195635/231‑1 (MQ=255)                                                                                                                                                                          gaCG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAATCTTAA‑G‑CGGATGTTT‑CGTTGCTTTATAGA‑AACCCGCCTCTT‑CGGCCACTTTGGCCATTTCAGCAGGAAGTAAAAGATCAAAAGGGGTGTCAGCTTTTACACTAACT‑CTCTCTTTTTTAATTCGGCGGCGGGATACTAACAAAACATTTAAAATGAGAAAATGATATAGAACATAAATCCCCTGGCTta                                      >  2:1372602/1‑239 (MQ=255)                                                                                                                                                                           aCG‑AAGGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑G‑CGGATGTTT‑CGTTGCTTTATAGA‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                                                             <  1:388841/215‑1 (MQ=255)                                                                                                                                                                            cG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑G‑CGGATGTTT‑CGTTGCTTTATAGA‑CACCCGCCTCTT‑CGGCCACTTTTGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGGTGTCAGCTTTCCCCCTAACT‑CTCTCTTTATTTAGTCGGCGACGCGGTACTAACAAAGCATTATAAATGAGAAAATTATATAGaaacaaactcccccggcttaaa                                    >  2:1014338/1‑218 (MQ=255)                                                                                                                                                                                     aTTGAG‑ATAAAACCCCGGGGGG‑AAACCGCC‑AGAAAGAAA‑TGTTAAGG‑GGGTTGTTT‑CGTGGTTTTAAAGC‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGt                             <  1:1014338/238‑1 (MQ=255)                                                                                                                                                                                               aaaaCCCCGGGGG‑TAAACCCCAAAATAAAAATTCTTAA‑G‑CGGTTGTTT‑GGTTGTTTTAAAGA‑CCCCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCAt                    <  1:1372602/239‑1 (MQ=255)                                                                                                                                                                                                 aaCACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑G‑CGGATGTTT‑CGTTGCTTTATAGA‑CACCCGCCTCTT‑CGGCCACCTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGGTGTCAGCTTTCACACTAAAC‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCCTTATAGATTAGAAACTGATATAGATCatat                                                   >  1:1138186/1‑205 (MQ=255)                                                                                                                                                                                                      cccggggatagtggaaaagatgAAAGTCTGAA‑G‑CGGTGTTTT‑CTTGGCTTTAGTGA‑CCCCCGCCACGT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGtttt              <  1:290818/218‑1 (MQ=255)                                                                                                                                                                                                                 cGCA‑AGATAGAAAGTCTTAA‑G‑GGGTGTTTT‑CTTTTCTTTATAGC‑CCCCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACg  <  1:770150/239‑1 (MQ=255)                                                                                                                                                                                                                 acccgcagaataaaagtttAA‑G‑GGGGTGTTTGTGTTGTTTTAAAGA‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATa    <  1:1028746/222‑1 (MQ=255)                                                                                                                                                                                                                  cgccagtgatataGTCTGAA‑G‑TGGGTGTTT‑CGTGGTTTGAGAGC‑CCCCCGCCTCGT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACg  <  1:721230/227‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                               GATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCG‑CCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACG‑AATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑G‑CGGATGTTT‑CGTTGCTTTATAGA‑CACCCGCCTCTT‑CGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACT‑CTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACG  >  NC_000913/954152‑954591

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: ‑