Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	414,575	C→T	12.4%	E393E (GAG→GAA)	sbcC ←	exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	414,575	0	C	T	12.4%	165.7 / 2.0	67	E393E (GAG→GAA)	sbcC	exonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (9/47);  new base T (6/2);  total (17/50)
Fisher's exact test for biased strand distribution p-value = 1.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACT  >  NC_000913/414338‑414809                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               ttctttATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCt                                                                                                                                                                                                                                               <  2:15849/239‑1 (MQ=255)     ttATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑c                                                                                                                                                                                                                                          <  2:988957/239‑1 (MQ=255)     ttATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑c                                                                                                                                                                                                                                          <  1:1027628/239‑1 (MQ=255)      tataACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑Ca                                                                                                                                                                                                                                         <  1:929570/239‑1 (MQ=255)      tataACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑Ca                                                                                                                                                                                                                                         <  1:172578/239‑1 (MQ=255)      tataACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑Ca                                                                                                                                                                                                                                         <  1:1252533/239‑1 (MQ=255)      tataACGCTGGCGCATTTCGGTAAGTGCGGCGTTACGGTGCGTCGGTTCTTGCGTGCCATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑Ca                                                                                                                                                                                                                                         <  2:357339/239‑1 (MQ=255)              tGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTTACGTCAACGTGCTCGCCGCAAGCG‑GCGTCA‑GTTTTTTTTCCG‑CATAGGGTTa                                                                                                                                                                                                                                 >  2:59624/1‑239 (MQ=255)              tGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGGTCGGCGCCAACG‑GATTAATGTTTTTG‑TTAG‑CAT‑GGGttta                                                                                                                                                                                                                                >  1:828182/1‑237 (MQ=255)              tGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTaa                                                                                                                                                                                                                                <  1:1209987/239‑1 (MQ=255)                         gTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGcc                                                                                                                                                                                                                     <  1:635751/239‑1 (MQ=255)                           tAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAATGGACACGCGATCTCCGCAAGGGCCATTAG‑GTTTTT‑CTTAG‑CAT‑GGGtttaatttttcgtcc                                                                                                                                                                                                                     >  1:515419/1‑224 (MQ=255)                                         cGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACCTGAATGCCGCCACCG‑CATTTA‑CGTTTTTTTTAG‑CCT‑GGGGTTACTGTttcttcgccttccccagctggtt                                                                                                                                                                                                     >  1:1240026/1‑220 (MQ=255)                                              cGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTTATGTCAAACTCATCGCCCCCAGCG‑CACTAA‑GGTTTTTGTCAG‑GCT‑GGGTTTACTCTTTCTGCCCTTGCCCCCGAtttccgcga                                                                                                                                                                                                 >  1:523819/1‑230 (MQ=255)                                                 cTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCg                                                                                                                                                                                              <  2:1232691/238‑1 (MQ=255)                                                  tGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑Ct                                                                                                                                                                                           <  2:336542/239‑1 (MQ=255)                                                  tGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑Ct                                                                                                                                                                                           <  2:438069/239‑1 (MQ=255)                                                   gTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAAAGTCAACGGGATCGCCGCCAGCG‑CACTAA‑GTTTTTTTTCAG‑CAT‑GGGGTTAATGTTGCTGGCATTGCCGCCCATGCTCTCGagctccg                                                                                                                                                                                           >  1:1165163/1‑233 (MQ=255)                                                     tCTTGCTTTACATTCTGGATAGCTACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGt                                                                                                                                                                                        <  2:1015725/239‑1 (MQ=255)                                                     tCTTGCGGTACATTCTGGCTAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGt                                                                                                                                                                                        <  2:1143803/239‑1 (MQ=255)                                                       ttGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGttt                                                                                                                                                                                      <  1:315324/239‑1 (MQ=255)                                                              aCATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAAAGACCACGGGTTCGCCGCAAGCG‑CATTTA‑GGTTTTTCTCTG‑CAT‑GGGTTAACTGTTTCGGCCCTTTCCGCAGAGGCGAGCGGTCGCCTGTTTTTTTggg                                                                                                                                                                                >  1:708421/1‑239 (MQ=255)                                                                       aTAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTgcg                                                                                                                                                                       <  1:785399/238‑1 (MQ=255)                                                                                  tagcTGCGTCAGGCG‑TTTGGTTGGGTGAACAATCTGTCCATGCAGCCCGCCCAGGTGCTGCCGCGGTGGGCTTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTCCATCGGCGGTTAACGTCATCGTGATCGCCGCAAGCG‑CATTAA‑TTTTTTGCTCAG‑CAT‑GGGTTAACTGGTGCTGCCATTGCCGCAGATGCTCGCTAGCG‑CTGGTTTGTTGGGAGAACTGCGCACGCAAACccg                                                                                                                                                           <  2:464511/236‑1 (MQ=255)                                                                                                  ttttGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCa                                                                                                                                           <  1:1236284/239‑1 (MQ=255)                                                                                                                     tttCCATGCGGCGCGACCAGGTGCTGCCGCAGTGGGCGTTGCCCGGCATGTTGCGCCAGGGGGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCgtg                                                                                                                         <  2:1165163/237‑1 (MQ=255)                                                                                                                     tgcccATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTGGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCgtg                                                                                                                         <  1:331143/236‑1 (MQ=255)                                                                                                                     tGTCCATGCAGCGCGACCAGGGGATGCCGCAGTGGGCTTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCACCCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCgtg                                                                                                                         <  2:377410/239‑1 (MQ=255)                                                                                                                     tGTCCAGGCCGCGCGCCCAGGTGCTGAGGCAGTGGGCGTTTCTCAGCATTTTGCGCCGGGGCGGTAGCAATTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAG‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCgtg                                                                                                                         <  2:798329/239‑1 (MQ=255)                                                                                                                     ggTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCgtg                                                                                                                         <  1:59624/238‑1 (MQ=255)                                                                                                                       tCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTTACCTCAACCGGGTCGCCGCCCGCG‑CCTTAA‑GTTTTTTCTCTG‑CCC‑GGGTTAACTATTTTTGGCCTTGCCCCAGATGCTCTCGATCG‑CTGGGTTGTTTGGGGAACTGTGCCCGCCCACCCGCCCGTTCGGTGTTCCCCTGGCGGGAGCGGGCGTGtt                                                                                                                       >  2:263566/1‑239 (MQ=255)                                                                                                                              aGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGACCTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑Ta                                                                                                                <  1:1104015/238‑1 (MQ=255)                                                                                                                               gcgcGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTTAAGGCAACCTGATCGCCCCCAGCG‑CCCTAA‑ATTTTTTTTTTG‑CAT‑GGGGTAAATGTTTGTTCCCTTTCCCCCGGTGCCTCCGATCG‑CCGGGTTTTTGGGGGGAATGCGGCCCCCCCCCCCCCCGGTCGGTGTTTCACTGTCCGAGGCGGGCGGGGTccctgta                                                                                                                >  2:331143/1‑233 (MQ=255)                                                                                                                                 gcgACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGGTTACGTCAACGAGGTCGCCCCCAGCG‑CCGTAA‑ATTTTTTTTCTGTAAG‑GGGTTACTTTTTTCCCCCCTTTCC‑CACAAGGTTCCGGTCT‑CTGGCTTGTTTGGTGAACCGCCCACGCCAACCCCCCCCGTCGTTGTTTCACTGACGGGAGCgggggcttttctggaaccc                                                                                                             >  2:1314059/1‑222 (MQ=255)                                                                                                                                         gTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAg                                                                                                                                 <  1:523639/211‑1 (MQ=255)                                                                                                                                          tcgtgcagCAGTGGGTTTTGCCTCGCATGTTGCGCCAGGGCGGTAGCAACTGCATCGGCGTTAAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATGGCGGGAGATGCTCGCGAGCG‑CTGGTTTGTTGGGAGAACTGCGCTCGCCACCCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTc                                                                                                     <  2:158536/230‑1 (MQ=255)                                                                                                                                          ggCTGACGCGGTCGGCGTTGCTCACCATGTTGCGCCAGGGCGGTAGCAACTTCATGGGCGTTTAACGACAACTTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTGGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCa                                                                                                    <  2:1098577/237‑1 (MQ=255)                                                                                                                                           gCTGACGCAGTGGGCGTTTCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGGGGTTAACGTCAACGTGAGCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAgg                                                                                                  <  1:231383/239‑1 (MQ=255)                                                                                                                                                 gCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTttgt                                                                                            <  1:603950/239‑1 (MQ=255)                                                                                                                                                   aGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGGTTACGGCAAAGTGATCGCCCCCAGCG‑CCCTTA‑ATTTTTTTTCAT‑CCC‑GGGTTTAATGTTTCTGGCCTTGCCGCCGATTCTCCCGCTCG‑CTGGTTTTTTTGGGGGACTTCGGCCGCCCACCCCCCCGGTCGGTGTTTCCCCGACGGGAGGGGGCGGGGTTCTG‑tctcttattcacatccttctctg                                                                                           >  2:523639/1‑216 (MQ=255)                                                                                                                                                      ggggtttGCTCGGCATGTTGCGCCAGGGCGGTAGCAACTTCAGCGGCGGTTAACGTCAACGGGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTtgctg                                                                                       <  2:1117790/234‑1 (MQ=255)                                                                                                                                                        gCGTTGCTCAGCAGGGTCCGCCAGGGGGGTAGAAAGTTCATGGGCGGTTAACGTAAACGTGATCGCCGCAAGCG‑AATTAA‑GTTTTTGCGCAG‑TAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTtgctgc                                                                                      <  2:646251/238‑1 (MQ=255)                                                                                                                                                            ggCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCgg                                                                                                                              <  1:674579/194‑1 (MQ=255)                                                                                                                                                               tCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTaa                                                                              <  1:622080/239‑1 (MQ=255)                                                                                                                                                                        tGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTTACGCCAACGTGAACCCCGCCAGCG‑CGTTAA‑GATTTTGCTCCG‑CCT‑TGGTTTACTTTTGTTGCCACTTCCGCCCACACCTCCCCTCG‑CTGGTTTTTTGGGTGGAAGCCGCACGCCCACCCCCCCGGTCGTTGTTTCTCTTCTGGGAGGGGGCGGGTGTCTGTTAACAAGGT‑TTCTGGGTTTGTtgcttcttgtgcgattcttctt                                                                      >  2:611835/1‑223 (MQ=255)                                                                                                                                                                             cacgacagagggaatttCACCGGCGTTTAACGGCACCGTGTGCGCGGCAAGGC‑CATTAG‑GTTTTTGACGAG‑GAT‑GGGTTACCTGTTGCAGCAGTTGCCGCAGATGCTCGCGGTCG‑CTGGTTTGTTGGGGGAACGGCCCGCGCGAACCCGCCGGTTCGTTGTTCCACTGACGGACGCGGTCGTGTTCCTG‑TAGCCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCGGCTGACTGCt                                                                <  2:76409/224‑1 (MQ=255)                                                                                                                                                                                    gTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTcgccg                                                          <  1:192353/238‑1 (MQ=255)                                                                                                                                                                                                         aaCGTCAACTTGATCGCCGCAAGCG‑CATGAA‑TTTTTGTTTCCT‑CCG‑GGTTTACCGGTTGCTGCCTTTGCCGCAGATGCCCGCGATCG‑CTGGTTTGTTGGGAAAACTGCGCCCGCCACCCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTcgcgcg                                    <  2:559376/239‑1 (MQ=255)                                                                                                                                                                                                            gTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCCCGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAg                                 <  1:1155514/239‑1 (MQ=255)                                                                                                                                                                                                            gTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAg                                 <  1:399914/239‑1 (MQ=255)                                                                                                                                                                                                            gTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CGGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAg                                 <  1:1336336/239‑1 (MQ=255)                                                                                                                                                                                                            cacagcGTGATGGCCGAGAGCG‑GATTAA‑GTGTTTCCGTGG‑CGT‑GGGGTAGCTGTTGTTGCCAGTGCCGGAGATGCTCGGGATCG‑CTGGTTTTTTGGGAGACCGGCGCGCCCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAg                                 <  2:410389/234‑1 (MQ=255)                                                                                                                                                                                                             gcaaggtgtgcgccgcAAGCG‑CATTAA‑GTTTTTGTCCGG‑CTG‑GGGTAAACGGTGGCTGCCATTGCCGAGGTTGCTGGCGAGCG‑CTGGTTGGTTGGGAGAACTGCGCCCGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAAGGCTCGCGCGAAg                                 <  2:820461/228‑1 (MQ=255)                                                                                                                                                                                                              cAACGTGATCGCCCCGAGCG‑CATTAA‑GTTTTTTCTCAG‑CCT‑GGGGTAACTGTTGCTGCCATTTCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGGGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGGATTCAGGCTTTGTTGCTGATGGTGTAATTCTTGTGACTGCTTCCCCCCGTGGTGGCGGAAGCTCCCGCGagga                                >  2:600532/1‑235 (MQ=255)                                                                                                                                                                                                                      tCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAAcc                                                                                                             <  1:1199004/153‑1 (MQ=255)                                                                                                                                                                                                                      tCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAAcc                                                                                                             >  2:1199004/1‑153 (MQ=255)                                                                                                                                                                                                                      tCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATtgtg                       <  1:646415/239‑1 (MQ=255)                                                                                                                                                                                                                      tCGCCGCAAGCG‑CATTAA‑GTTTTGGCTCAG‑CGT‑GGGTTAACTGTTGCTGCCATGGCCGCAGATGCGCGCGAGCG‑CTGGTTGGTGGGGAGAACGGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATtgtg                       <  1:470304/239‑1 (MQ=255)                                                                                                                                                                                                                      tCGCCCCAAGCG‑CATTAA‑GTTTTTGCTCGG‑CTT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATtgtg                       <  2:1240026/239‑1 (MQ=255)                                                                                                                                                                                                                      ccgccgcAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAATTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATtgtg                       <  1:73196/238‑1 (MQ=255)                                                                                                                                                                                                                        gccgcAAGCT‑CATTAA‑GTTTTTGCCCAG‑CAT‑GGGTTATCTGTGGCTGCCATTGCCGCAGATGCGCGCGATCG‑CTGGTGTGTTGGGAGAACTGCCCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGct                     <  2:1146446/239‑1 (MQ=255)                                                                                                                                                                                                                         ccgcAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTTGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGCATGCTCGCGCGAAGCGCCATtgtgg                      >  2:904462/1‑236 (MQ=255)                                                                                                                                                                                                                            cAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGt                 >  1:343760/1‑239 (MQ=255)                                                                                                                                                                                                                                 g‑cATTAA‑GTTTTTGCTGAG‑CAG‑GGGTTAACTGTTGCTGCCATTGCGGCGGATGCTCGCGATCG‑CTGTTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAg              <  1:1145363/237‑1 (MQ=255)                                                                                                                                                                                                                                 g‑cATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAg              <  1:1062325/237‑1 (MQ=255)                                                                                                                                                                                                                                   gcttaaggTTTTTCCTCGG‑CTG‑GGTTTAACTGTTGCGTCCTTTCCCCCGGATTCCCGCGTTGC‑CGGGTTTTTTGGGAGACCTCGGCACCCCACCCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGc             <  2:1140687/231‑1 (MQ=255)                                                                                                                                                                                                                                               tGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTAct  >  2:143146/1‑237 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               TTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGACATTCTGGATAGCGACCTGTAACTGCGCCAGACGTTTTTGTTGGG‑GAACAATCTGTCCATGCAGCGCGACCAGGTGCTGACGCAGTGGGCGTTGCTCAGCATGTTGCGCCAGGGCGGTAGCAACTTCATCGGCGGTTAACGTCAACGTGATCGCCGCAAGCG‑CATTAA‑GTTTTTGCTCAG‑CAT‑GGGTTAACTGTTGCTGCCATTGCCGCAGATGCTCGCGATCG‑CTGGTTTGTTGGGAGAACTGCGCACGCCAACCCGCCGGTTCGTTGTTCCACTGACGGAAGCGGTCGTGTTCCTG‑TAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACT  >  NC_000913/414338‑414809

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