Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,830,350	T→A	19.1%	intergenic (+15/‑106)	xanP → / → yicH	xanthine permease/putative inner membrane‑anchored periplasmic AsmA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,830,350	0	T	A	19.1%	120.9 / 5.3	69	intergenic (+15/‑106)	xanP/yicH	xanthine permease/putative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base T (39/13);  new base A (2/12);  total (44/25)
Fisher's exact test for biased strand distribution p-value = 5.79e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CCATCGCCGCCTCCGGTGTGCGTATCGTTTCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGT‑GTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGAGCAGAACATATCAGCGCATGGGTTTCCGAGA  >  NC_000913/3830112‑3830579                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                         ccgtgccgccccgcGGTGTGCGGTTCGTGTCTGGGTAGGCGCGGAACAGTCGGGGGTTTCTGTTTATCGCGTTGTCGCTGGCGGTTGGTCTGGGCTT‑GTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTTCTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGt                                                                                                                                                                                                                                          <  1:819563/226‑1 (MQ=255) ccctCGCCGCCTCGGGTTTGAGTTTTTGTTCTCGGGAGCCGCGGAACCGGCGGGCGTTTGTGATTAGCGCGCTGTCGGTGGCGGTTGGTCGGGGGGT‑GTCTCAGCCGCCGCTGATGTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGt                                                                                                                                                                                                                                          <  1:336874/236‑1 (MQ=255) ccctCGCCCCCTCCGTTGTGCGTATCGTTTCTCGTGAGCCGCTGAACCGGCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGT‑GTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGt                                                                                                                                                                                                                                          <  2:892057/236‑1 (MQ=255) ccATCGCCGCCTCCGGTGTGCGTATCGTTTCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGT‑GTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGt                                                                                                                                                                                                                                          <  2:1193445/239‑1 (MQ=255) ccATCGCCGCCTCCGGTGGGCGTATCGTTTCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGT‑GTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGt                                                                                                                                                                                                                                          <  2:381284/239‑1 (MQ=255)  cATTGCCGCCGGCGGTGAGCGGATCGTTTCGGGTGAGCCGGTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGT‑GTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTa                                                                                                                                                                                                                                         <  2:852267/239‑1 (MQ=255)    tCGCCGCCTCCGGTGTGCGTATCGTTTCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGT‑GTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGGAAAACCTGCTCTCCTCCGGGAT‑CGCCCCGGGCGGGATTACTGCCATCGTGCTGAATCTGATTTTCCCCCC‑CCGAAAACAAGAA‑TTCACTGTCGGCGGTaca                                                                                                                                                                                                                                       >  2:819563/1‑237 (MQ=255)           tttcgtgtggggtttCGTTTCTGGGGAGCCGGTGACCGGTGGGGGTTTTCTGGTTTTCGCGCTGTCGCGGGCGGTGGGTCGGGGGGT‑GTCTCAGCATGCTTTGATTTTGAAGTTGGCCCGGGAAGGGCTGAAAAACCTGCTCCCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATGGTGCTGAATCTGTTTTCCCCCCC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAGGATAATTa                                                                                                                                                                                                                                <  1:158752/225‑1 (MQ=255)                  gtgCGTATCGTTTCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGGGG‑GTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCCCGGGCGGTATTACTGCCATCGGGCTGAATCTGATTTTCCCCCC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGGAAGGAAAATAccgcccc                                                                                                                                                                                                                         >  2:1268434/1‑237 (MQ=255)                         ttgTTTCTCGTTGCCCGGGGAGCCGGGGGGGGATTTTGTTTATCGCGTTGTCTGGTGCGGTTGGT‑TGGGGGTGGTCCCACCACCCGCTGTTTTTGCAGTTGGCCCTGGAATGGTTGAAAACCTTTCCCTCCTCGGGGAC‑CGCGGCGGGCGGTTTTACTGCCAGCGTGCTGAATCTGTTTTTCCCACC‑AGAAAAGCAGTAT‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCa                                                                                                                                                                                                                 <  1:1195294/237‑1 (MQ=255)                             ttCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGT‑GTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGGATTACTGCCATCGTGCTGAATCTGATTTTTCCACC‑AGAAAAACAGTAA‑TTCACTAGCGGCGGTAAGATAATTATCCCCGCC‑GTCTCAcc                                                                                                                                                                                                               >  2:660956/1‑237 (MQ=255)                                          cTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGGCTGGGCGT‑GTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGGATTACTACCATCGTGCGTAATATGATTTTCCCACC‑CGGAAAAACGTAA‑TTCACTGGCGGGCGTAAAGAAATTTCCGCCCCC‑GTCTCCCCCCTTTTAACCATTc                                                                                                                                                                                                 >  2:774995/1‑238 (MQ=255)                                           tctcctgtggggtttttCTGTTTGTCGGGGTGTCGCGGGCGGTGGCCTTGGGTGT‑TTCT‑GGCCGGCCCTGTTTTTGCCGTTTGCCCCTGAAAGGGCGAAAAACCTGTCCCCCCGCGGTATCCGCCGGGGCGGGTATTCCTGCCATCGGGTGGATTTTGTTTTCCCCCCC‑AAAAAAAAATAAT‑TTCCCTGGGGGCGGTAATATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCg                                                                                                                                                                                              <  1:1065506/224‑1 (MQ=255)                                                                    cgcgCTGTCGCTGGCGGTTGGTCTGGGCGT‑GTCTCAGCAGCCGCTGATTTTTCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCCGGAT‑CGCCGCGGGCGGGATTACTGCCATCGTGCTGAATCTGGTTTTCCCACC‑AGAAAAACAGGAA‑TTCACTGGCCGCGGAAAGGAAATTAACCCCGCAGGTTCCACCTC‑TTT‑ACCCTTCCCCCCtctgggatttccccgtaa                                                                                                                                                                          >  2:1204895/1‑217 (MQ=255)                                                                    cgcgCTGTCGCTGGCGGTTGGTCTGGGCGT‑GTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGGAAGATAATTACCGCCGCA‑GTCTCACCCC‑TTTTACCATTCCCGCCCTGAGGGTTTCGCGGaaaacct                                                                                                                                                                      >  2:995033/1‑234 (MQ=255)                                                                    cgcgCTGTCGCTGGCGGTTGGGCTGGGCGT‑GTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCCCCGGGGT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGGATCTGAGTTTCCCACC‑CGACAAACCGTAC‑TTCACTGCCGGCGGGAAGATAATTATCACCGCA‑GTACCACCTC‑TTTACCACTTccccccattacgagacggcctgaaaaagg                                                                                                                                                                     >  1:256432/1‑216 (MQ=255)                                                                                                gttttcccgcggcctgCTTATTTTGAGTTTTCCCCCTGAAGGGTAAAAACCCCTGCTCCCCCCGGGGAT‑CGCCGCGGGGGTTATTTCTGCCATCGTGCTGAATCTGATTTCCCCCCC‑AGAAAAACAGTAA‑TTCACTGGGGGCGGAAAGATAATTCCCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACg                                                                                                                                           <  2:657493/223‑1 (MQ=255)                                                                                                 t‑gtCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGCGGATTGCGCGTAAATCGTGCCTAACTCCCTTTTGTGCATTTCACGGGa                                                                                                                                        >  2:1216814/1‑239 (MQ=255)                                                                                                    tctcAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTCCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTACCGGGATg                                                                                                                                      >  2:1289413/1‑239 (MQ=255)                                                                                                    tctcAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGCCGGATGAATAATTACCGCCGCA‑GACTCACCCT‑TTGAACCATTCCCGCCTTTGGGATTGCCGACAAATCGTGCATACATCCCTTTTGTACCTTTAACggggag                                                                                                                                      >  1:750860/1‑236 (MQ=255)                                                                                                            gcCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATTTGCCTTTCACgg                                                                                                                                          >  2:823935/1‑227 (MQ=255)                                                                                                             ccGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGCTGGAAGACCAt                                                                                                                             >  2:761723/1‑239 (MQ=255)                                                                                                               gCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGa                                                                                                                                 >  2:502638/1‑233 (MQ=255)                                                                                                                             ttGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCC‑GTCTCAACTC‑TTTAACCATTCCCGCCTTGAGGATTTCGCGTAAATCGTGCCTAACTCCCTTTTTGTCCTTTCACGGGTTGGAAGACCATGAAATTTATTGGGAAg                                                                                                             >  1:1116361/1‑239 (MQ=255)                                                                                                                              tGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTgggga                                                                                                              >  2:493195/1‑235 (MQ=255)                                                                                                                                      aaTGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTtc                                                                                                      >  2:718123/1‑237 (MQ=255)                                                                                                                                           cTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTCTGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACAt                                                                                               >  2:815526/1‑239 (MQ=255)                                                                                                                                            tGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGa                                                                                                                                 >  2:581401/1‑204 (MQ=255)                                                                                                                                             gAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGGAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGAGGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACAt                                                                                               >  2:770519/1‑237 (MQ=255)                                                                                                                                              aaaaaCCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGCGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAACTGCTTCtataacttc                                                                                             >  2:1263230/1‑233 (MQ=255)                                                                                                                                                aaaCCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTaaa                                                                                                                                                                          >  2:8324/1‑159 (MQ=255)                                                                                                                                                   cctgcctccccccGGGAC‑CCCCGGGGGGGGTATTCCGGCCCTGGTGCTAAATCGGTTTTCCCCCCA‑AGAAAAACAGTAA‑TTCCTTGGCGGGGGAAAAATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCg                                                                                       <  2:149936/228‑1 (MQ=255)                                                                                                                                                    cgtctctCCCCGGGTAT‑CGCCGGGGGGGGTATTACTGCCATCGTGCTGAATTTGTTTTTCCCACA‑AAAAAAAATGTAA‑TCCGTGGGGGGCGGTAAGATAATTCCCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGc                                                                                      <  2:1051896/236‑1 (MQ=255)                                                                                                                                                    cTGCTCCCCCCGGGTAC‑CCCCGGGGGGGGTTTTACTCCCTTCTGGGTAAATCTGATTTCCCCACC‑AAAAAAACAGTAA‑TTCCTGGGGGGGGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGc                                                                                      <  2:799093/239‑1 (MQ=255)                                                                                                                                                     tGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTTTCTACATTCCCATCGct                                                                                     >  2:984858/1‑239 (MQ=255)                                                                                                                                                     tGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCGTCTTCTCTACATTCTCATCg                                                                                       >  2:683891/1‑237 (MQ=255)                                                                                                                                                          tcctccGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCCCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGaa                                                                                                                                   >  2:625146/1‑188 (MQ=255)                                                                                                                                                            ctccGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCCCATCGCTCTGTTAg                                                                              >  2:348090/1‑239 (MQ=255)                                                                                                                                                               cGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTgg                                                                           >  2:261592/1‑239 (MQ=255)                                                                                                                                                                       gccgcGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTTCGCGGAAATCGTGCAAAACTTCCTTATGTGCATTTCACGGGGTGGAAGACCATGGAATTTATTGGGAAGCTGGTTTTTTACATTATCATCGCTCTTGTTGTGGccgctgc                                                                    >  1:1261610/1‑233 (MQ=255)                                                                                                                                                                       gccgcGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACCTTCTCATCGCTCTGTTAGTGGCGAT‑TGc                                                                   >  2:530748/1‑239 (MQ=255)                                                                                                                                                                        ccggggggggttttccgccattCGTGCTGAATCTGATTTCCCCACC‑AAAAAAACAGTAA‑TCCACGGGGGGCGGAAAGATAATGCCCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCt                                                                  <  2:1176602/218‑1 (MQ=255)                                                                                                                                                                            ggggggTATTCGTGCCTTGGTGCTGATTTTGATTTCCCCCCA‑AAAAAAACAGTAA‑TTCACGGGCGGCGGAAAGATATTTCCCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGc                                                               <  2:792054/238‑1 (MQ=255)                                                                                                                                                                              gcgtTATTACTGCGATCGTGCTGAATCTGTTTTCCCCACA‑AAAAAACAAGAAA‑TTCACGGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCct                                                             <  2:1048059/238‑1 (MQ=255)                                                                                                                                                                                      tcctccctgctgtcaaaTCTGATTTTCCCCCCAAAAAAAATATGAA‑TTCACGGGGGGGGGTAAGATATTTCCCGCCGCA‑GTCTCCCCTT‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTttc                                                     <  1:761723/224‑1 (MQ=255)                                                                                                                                                                                          ccATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGTCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGCGAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTAACTGGCGAA‑CACAGGCTTCTATTttctttac                                                >  1:294302/1‑236 (MQ=255)                                                                                                                                                                                            aTCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGTCGGCGGTAAGATCATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTGAATCTTGCATCACTACCTTAAGTGCATTTCACGGGATGGAAAAACCTGAAATTTTCTGGGAAGCTGCTTCTCTAACTTCTCATCCCCCCGTTAATGGCGGT‑CACCTGCCTCTATTTTTTTTTGCa                                              >  1:901368/1‑239 (MQ=255)                                                                                                                                                                                            aTCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCa                                              >  2:368997/1‑239 (MQ=255)                                                                                                                                                                                            aTCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CCCTGGCCTCTATTTTTTTCTTCa                                              >  2:1128734/1‑239 (MQ=255)                                                                                                                                                                                            aTCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCCCGTAAATCGTGCATAACTCCCTTATTTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCCTCGCTCCGTTAGTGCCGAT‑CCCTTGCCCCTTTTttctttctca                                              >  1:841145/1‑234 (MQ=255)                                                                                                                                                                                              cGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCTTGCCTAACTCCCTTATGTGCATTTCACAGGATGGAAGGCCCTGAAATTTATTGGGAAACTGCTTCTCTACATTCTCATCCCCCTGTTAGTGACGAT‑AGCTGGCCTCTATTTttttctccca                                             >  2:702873/1‑234 (MQ=255)                                                                                                                                                                                              cGTGCTGAATCTGATTTTCCCACA‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTAAGGATTGCCCGTAAATCGTCCATAACTCCCTTATATGCATTTCACAGGATGGCAGACCATGAAATTTATTGTGAATGTGCTTCTCTACATTCTCATCGCCATGGTAGTGGCCAT‑CGCTCGCCCCTATTAttttcttaaa                                             >  1:305953/1‑234 (MQ=255)                                                                                                                                                                                                tgttgattctatttcccccACC‑AAAAAAACATAAA‑TTCACTGGCGGTCGAAAGAAATTTACCGCCGCT‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAAcc                                          <  2:1056286/223‑1 (MQ=255)                                                                                                                                                                                                tgcgtaaTTTGTTTTCCCCACG‑AAAAAAACTTTAT‑TGCATTGGGGGGGTAAAGATAATTCCCGGCCGA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAAcc                                          <  1:995033/234‑1 (MQ=255)                                                                                                                                                                                                     aaTCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGAGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGTAGACCATGAAATTTCTTGGGAAACTGCTTCTCTCAATTTGCTTCGCTCTGTTAGTGGGGATCCGCT‑GACTCTATTTTCTTCTGCAAACCCCCtt                                     >  1:877222/1‑238 (MQ=255)                                                                                                                                                                                                      aTCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGGAGACCATGAAATTTCTTGGGAAGCTGCTTCTCTACCTTCTCATCGCACTGTTAATGGCGTT‑CCATGGCCTCACTTTTCTTCTGCAAACCCGCTgg                                    >  2:580015/1‑239 (MQ=255)                                                                                                                                                                                                       atgtgttttCCCCCCAAAAAAAACAATAA‑TTCAGTGGGGGCTGTAATAAAATGCCCGCCGCG‑TTCTCACCTC‑GTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAACCCGCTg                                     <  2:512538/233‑1 (MQ=255)                                                                                                                                                                                                          taTTTCCCCCCA‑AAAAAAACAGAAA‑TGCACGGGGGGCGGAAAGAAAATTCCCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGAg                                <  2:386157/238‑1 (MQ=255)                                                                                                                                                                                                          gATTTCCCCCCC‑AAAAAAACGTTAT‑TTCCGTGGGGGCGGAAAGATAATTCCCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGAg                                <  2:1116361/239‑1 (MQ=255)                                                                                                                                                                                                           gttttCCCCCCAAGAAAACCAGTAA‑TCCCCTGGCGGGGATAGGATAATTACCGCCGCA‑GTCTCACCTC‑TTTACCCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGAg                                <  2:829832/238‑1 (MQ=255)                                                                                                                                                                                                           gtttcCCCACC‑AAAAAGAAAGTAAGTCCCCGGGCGGCGGTAGGATATTGCCCGCCGCA‑GTCTCACCAC‑TTGAGCCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGTGATGGAAGACCAGGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGAg                                <  2:656751/238‑1 (MQ=255)                                                                                                                                                                                                               ccccACC‑AAAAAAACATTAA‑TTCGTGGGCGGGGTTAAGATAATTCCCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGAGCAGaa                           <  2:474545/238‑1 (MQ=255)                                                                                                                                                                                                                cccACC‑AGAAAAACAGTCA‑TTCACTTGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCACGCCTTGAGGATTGCGCGTAAATCGTGCATCACTCCCTTATGTGCATTTCACGGGATGGAAGACCCTTAAATTTATTGGGAAGCTGTCTCTCTACATTATCATCGCTCTGTTTGTGGCGAT‑CCCTGGGCTCTTTTTTTTTTCACAAAcccacgtgggggaagaca                          >  1:334601/1‑224 (MQ=255)                                                                                                                                                                                                                   acc‑aGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATAGCTCTGTTAGTTGCGAT‑GGCTGGCATCTATTTTTTTCTTCAAAACCCCTGGGGAGCaaaaataa                       >  2:544100/1‑235 (MQ=255)                                                                                                                                                                                                                    cc‑aGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGGGCAGAACata                       >  2:760863/1‑238 (MQ=255)                                                                                                                                                                                                                    cc‑aAAAAAACGTTAT‑TTCAGTGGGGGCGGAAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGAGCAGAACatat                      <  2:22340/239‑1 (MQ=255)                                                                                                                                                                                                                     acaaaaaaaCAGTAA‑TCCCGTGGCGGGGGAAAGAAAATGACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGAGCAGAACatat                      <  2:841145/235‑1 (MQ=255)                                                                                                                                                                                                                           aaatatata‑TTCTCGGGCGGCAGAAAGAAAGTTCCCTCCGCG‑TTCTCTCCTC‑TTTAACCATTTCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGAGCAGAACATATCAgcgc                <  1:815526/231‑1 (MQ=255)                                                                                                                                                                                                                                gTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGGGCAGAACATATCAGCGCCTggg           >  2:615079/1‑239 (MQ=255)                                                                                                                                                                                                                                 taatttcccggggggcggTAAGATATTCCCCGCCGAA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGAGCAGAACATATCAGCGCATggg           <  2:1261610/226‑1 (MQ=255)                                                                                                                                                                                                                                   attccACTGGCGGGGTAAAGATATTTCCCGGCGCA‑GTCTAACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGAGCAGAACATATCAGCGCATGGGtt         <  2:953898/235‑1 (MQ=255)                                                                                                                                                                                                                                       cACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCAGAACTCCCTTATGCGCATTTCACGGGATGGAAGACCATGCAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCCCTCTGTTAGTGACGTT‑CGCTGGCCTCTATTTTCTTCTGCAAACCCACCGGGGGGCAGAACAAAACAGCGAATGGGTTTccc     >  1:946614/1‑238 (MQ=255)                                                                                                                                                                                                                                          tGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATTTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCTCTCTGTTAGTGGCGAT‑CGCTGGCCTCTTTTTTTTTCTTCAAACCCCCTGGGGGGCCGAACATATAAGCCCGTGGGTTTcccaaa  >  1:885426/1‑235 (MQ=255)                                                                                                                                                                                                                                          tGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAACCCGCTggg                                   >  2:1029160/1‑206 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                         CCATCGCCGCCTCCGGTGTGCGTATCGTTTCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGT‑GTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGAT‑CGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACC‑AGAAAAACAGTAA‑TTCACTGGCGGCGGTAAGATAATTACCGCCGCA‑GTCTCACCTC‑TTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGAT‑CGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGAGCAGAACATATCAGCGCATGGGTTTCCGAGA  >  NC_000913/3830112‑3830579

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