Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,043,014	A→G	23.5%	G200G (GGA→GGG)	srkA →	Cpx stress response Thr/Ser protein kinase; MazF antagonist protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,043,014	0	A	G	23.5%	137.9 / 9.0	69	G200G (GGA→GGG)	srkA	Cpx stress response Thr/Ser protein kinase; MazF antagonist protein
Reads supporting (aligned to +/- strand):  ref base A (9/41);  new base G (14/3);  total (24/45)
Fisher's exact test for biased strand distribution p-value = 3.35e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATGGAAGCGGTTGGGCGTTATTTAGGGC‑GTATGCA‑CCAGACGGGGCGCAAACAGCTTTTTAT‑CCATCGCCCGACCATCGGTTTGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGC  >  NC_000913/4042776‑4043250                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  gggttgAAGCGGTTGGGCTTTTTTGGGGGC‑GTA‑GCACCCAGCAGGGGCGCAAAACGTTTTTTTC‑CCAGCGCCCTCCCATCGGTTTAAAGGAATATCCC‑ATTGGGCCACGCAAGGTGTTTGGGGACGCTACACTGATCCCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACGGATGAGCTGTTTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGga                                                                                                                                                                                                                                                   <  2:849145/234‑1 (MQ=255)     ggAAGCGGTTGGGCGTTATTTAGGGC‑GTATGCA‑CCAGACGGGGCGCAAACAGCTTTTTAT‑CCATCGCCCGACCATCGGTTTGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGAGGG‑CAGGATTCTCCGTCCTGCGGGTCCATGgtg‑aa                                                                                                                                                                                                                                               >  1:622803/1‑237 (MQ=255)         gCGGTTGGGCGTTATTTAGGGC‑GTATGCA‑CCAGACGGGGCGCAAACAGCTTTTTAT‑CCATCGCCCGACCATCGGTTTGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGAGTTCACCGTTTTGCGGGTACATGgag‑aagtgc                                                                                                                                                                                                                                           >  1:207702/1‑233 (MQ=255)             ttGGGCGTTATTTAGGGT‑GTATCCA‑CCAGACGGGGCCCAAAACGTTTTTTAC‑CCATCCCCCGACCATCGGTTTGAAAGAAGATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTCCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑c                                                                                                                                                                                                                                      <  1:366209/238‑1 (MQ=255)             ttGGGCGTTATTTAGGGC‑GTATGCA‑CCAGACGGGGCGCAAACAGCTTTTTAT‑CCATCGCCCGACCATCGGTTTGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCCCACTGGCGGGAAAGATTTC‑CCGTTCTGCGGCCTCAAGGGG‑ACTTcccagcc                                                                                                                                                                                                                                      >  1:227526/1‑234 (MQ=255)             ttGGGCGTTATTTAGGGC‑GTATGCA‑CCAGACGGGGCGCAAACAGCTTTTTAT‑CCATCGCCCGACCATCGGTTTGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTTACCGTTTTTCGGCTCCAAGGGG‑ACTGCCAcc‑cc                                                                                                                                                                                                                                     >  1:843320/1‑239 (MQ=255)             ttGGGCGTTATTGAGGGC‑GTTTGCA‑CCAGACGGGGCGCAAACAGCTTTTTCT‑CCATCGCCCGACCATCGGTTTGAACGAATATCTT‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑c                                                                                                                                                                                                                                      <  1:381800/238‑1 (MQ=255)              ggggCGTTA‑TTAGGGCGGTATGCA‑ACAGTCGGGGCTCAAAAAGCTTTTTAC‑CCATCGCCAGACAATGGGTTTGAACGAATTGCTC‑ATTGAGCCAGGCAAGCTGTTTGAGGACGCTACACTGATCCCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACATATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCg                                                                                                                                                                                                                                    <  1:1045981/238‑1 (MQ=255)                 gCGTTATTTAGGGC‑GTATGCA‑CCAGACGGGGCGCAAACAGCTTTTTAT‑CCATCGCCCGACCATCGGTTTGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑GAGAGTTTACCGTTTTGCGGCTCCACAGGG‑AGTGCCACG‑CCcggg                                                                                                                                                                                                                                 >  1:256538/1‑238 (MQ=255)                            ggC‑GTATGCA‑CCAGACGGGGCGCAAACAGCTTTTTAT‑CCATCGCCCGACCATCGGTTTGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCCCACTGGCGGG‑CAGAATTCAACGTGCTGCGGCAACCAGGAG‑GAACCCACG‑CAGGGAATAAtcttcgg                                                                                                                                                                                                                      >  2:15743/1‑235 (MQ=255)                                 tATGCA‑CCAGACGGGGCGCAAACAGCTTTTTAT‑CCATCGCCCGACCATCGGTTTGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTAACCGTTCTCGGGGCACATGGGG‑AGTGTCACG‑CCAGGAGAATTTTCTGCcgcg                                                                                                                                                                                                                  >  1:220150/1‑239 (MQ=255)                                  tgtcA‑CCAGGCGGGGGGCCAACCG‑TTTTTATACCTGCGCCCGACCGT‑GGGTTGAACAGAAATCCCTATTGGGCAACGAAGGCTTTTGGAGGACGCTACACTGTTCCCTACCGGTTTGAAAGCGGAATTCCGAAAAGG‑GAAAGATGGGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGcgcg                                                                                                                                                                                                                  <  2:757850/235‑1 (MQ=255)                                  aTGCA‑CCAGATGGGTCGCAAACTGCTTTTTAT‑CCAGCGCCCGACCTTCGGTTTGAACGAATATCTC‑ATTGAGCCCCGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGa                                                                                                                                                                                                                 <  1:496735/239‑1 (MQ=255)                                  aTGCA‑CCAGACGGGGCGCAAACAGCTTTTTAT‑CCATCGCCCGACCATCGGTTTGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTTTGCGGCCACACGGGG‑ACTGCAACA‑CCCGGGGAATTCTCTgggtcga                                                                                                                                                                                                                 >  1:70152/1‑234 (MQ=255)                                  aTGCA‑CCAGACGGGGCGCAAACAGCTTTTTAT‑CCATCGCCCGACCATCGGTTTGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCAACCTTCTGCGGCTACATTGAG‑AATGCCCCC‑CCCGGGAGATTCTCTGGcgcgc                                                                                                                                                                                                                 >  1:1085401/1‑238 (MQ=255)                                   tGCA‑CCAGACGGGGCGCAAACAGCTTTTTAT‑CCATCGCCCGACCATCGGTTTGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTTACCGGTCTGCGGCTACATGGGG‑ACTGCCCCG‑CCGGGAATATTTTCTGGcgcgcg                                                                                                                                                                                                                >  1:1010910/1‑237 (MQ=255)                                         cAGACGGGGCGCAAACAGCTTTTTAT‑CCATCGCCCGACCATCGGTTTGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAAGATTT‑AACGGTCTGCGGGTACCAGGGG‑ACTGCCCCC‑CCCGGGATATTTTCTgggcgccggggcc                                                                                                                                                                                                           >  1:853655/1‑228 (MQ=255)                                          cgACGGGGCCAAAATCGTTTTTTCT‑CCATCGCCCGCACATCGGTTTGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATCCCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCa                                                                                                                                                                                                          <  2:256538/238‑1 (MQ=255)                                               gggCCCAACCGGTTTTTTTC‑CCCTCCCCCGGCCTTGGGTTTGAACGAATTCCCC‑TTTGGCCCACGCAAGCTGTTTGGGGACCCTATACTGATCCCTTCCGGGTTGAAAGCGTCATTCCGGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑t                                                                                                                                                                                                       <  2:207702/236‑1 (MQ=255)                                                 gcgcAAACAGCTTTTTAT‑CCATCGCCCGACCATCGGTTTGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑GAAAATTTACCCTTCTGCGGCTACATGGGG‑ACTGCCACG‑CCCGGGGAATTTTCCGGCGCGCTGGTTCAA‑tggttg                                                                                                                                                                                                  >  1:589702/1‑239 (MQ=255)                                                    caaaacgCTTTTTAC‑CCATCCCCCGGCCATCGGTTTGAACGAATATCCC‑ATTGGGCCCCGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑tgtttgtt                                                                                                                                                                                                <  1:1025971/232‑1 (MQ=255)                                                          gCTTTTTAT‑CCATCGCCCGACCATCGGTTTGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAAATTTCACCGCTCTGCGGCTACATGGGG‑AATGCCACC‑CCCGGGCAATTTCCTGGGCCGACGGTCCAC‑TGTtttttgtctcgt                                                                                                                                                                                         >  1:319523/1‑233 (MQ=255)                                                                  acccAGCGCCCGACATGCGGTTTGACCGAATATCTC‑ATTGAGCCACGAAAGCTGTTTGAGGACGCTACACGGATCCCTTCGGGTTTGAAAGCGGCATTCCTGAAAGG‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGc                                                                                                                                                                                  <  1:728063/237‑1 (MQ=255)                                                                        cGCCCGACCTTCGTTTTGAACGAATATCTC‑TTGGAGCCACCAAGGTTTTTTGAGGCCGCTACCCGGATGCCTTCCGGGTGGAAAGCGGCATTCCTGAAGGC‑GACAGATGAGCTGATGGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTaa                                                                                                                                                                            <  2:319523/239‑1 (MQ=255)                                                                           ccGACCATCGGTTTGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTGCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGa                                                                                                                                                                                                                                 <  1:463800/184‑1 (MQ=255)                                                                                cctCGTTTTGAACGAATATCTC‑ATTGGCCCACGGAAGTTTGTGGAGGACGCGACATTGATCCCTCCTGGGTTGAAGGCGGCATGCCGAAAAGC‑GACAGTGTAGCTGATTGCCGCCGTTACAGCACCCTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAg                                                                                                                                                                    <  2:220150/237‑1 (MQ=255)                                                                                      tttGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTCCACTGATCCCTTCGGGGTTGAAGGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTc                                                                                                                                                              <  1:446374/239‑1 (MQ=255)                                                                                       ttGAACGAAATCCCC‑ATGGGGCCACGAAAGCTGTTGGAGGACGCTCCACTGATCCCTTCCGGGTTGAAAGGGGCATTCCTGAAAGC‑GACAGAGGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTc                                                                                                                                                              <  1:965985/238‑1 (MQ=255)                                                                                         gAACGAATATCCA‑ATGGGGCCACGCAAGCTGTTTGGGGACGCTACACTGATCCCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑gg                                                                                                                                                          <  1:382113/239‑1 (MQ=255)                                                                                               ttttctc‑TTTGACCCCCGCAAGCTGTTGTAGGACCCTACCTTGATTCCTTCGGGTTTGAAGGCGGCTTTCCGAAAGGC‑GACAGTTGAGTTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTtgt                                                                                                                                                    <  2:853655/236‑1 (MQ=255)                                                                                                     taaTTGAGCCCGGCAAGGTGTTTGAGGACCCGAGACTGATGCCTTCCGGGTTGAAAGCGGCATTCCGGAAGGC‑GACAGATGAGCTGATGGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATg                                                                                                                                               <  1:409899/237‑1 (MQ=255)                                                                                                         tGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTACTGAAAGC‑GACAGATGAACTGATTGCCGCCGTTACAGCACCCGGGCGGG‑AAACATTTAACGCTATGCGGGTACAGAGAGAACTAGCACA‑CGGGGAAAATAC‑CTGCTGGCGTGGGCAAA‑TGTTTGTTGGAATGGATGAAGCACGTAATGGGGCCGCCGGCCC‑GGAGTTGTgggggtggg                                                                                                                                           >  2:185596/1‑232 (MQ=255)                                                                                                           aGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑GAGATTTTACCGTTCTGCGGCTAAATTGGG‑ACTGCCCCG‑CCGGGGATATTTTCTGGGGCGCTGGGCCCA‑TGTTTGTTTATCTGGATGATGGACGGAATGGTCCAGCCGTTCA‑GGATTTTTGGGTG‑TTGTTc                                                                                                                                        >  1:167459/1‑239 (MQ=255)                                                                                                                          ttGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑GAGATTTTACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTTTCTGGCGCGATTGTTCAA‑TGTTTGTTTGTCTGGGTGATGCACGGAATGGTCCAGCCGTTCA‑GGGTTTTTTGATG‑TTGTTCAATGGGGTTAAAAcc                                                                                                                         >  1:829308/1‑239 (MQ=255)                                                                                                                               gACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTTACCGTTCTGCGGCTACATGGGG‑ACTGCCACG‑CCGGGGATATTCTCTGGCGCGAGGGGCCAA‑TGTTTGTTGATCTGGATGATGCACGGAATGGTCCAGCCGTTCA‑GGATTTTTGGATG‑TTGTTCAATGGGGGTAAAAcccggca                                                                                                                    >  1:129724/1‑234 (MQ=255)                                                                                                                               gACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGGATATTCTCTGGCGCGCTGGTCCCA‑TGTTTGTTTATCTGGATGATGCACGTAATGGTCCAGCCGTGCC‑GGGTTTTTGGGTT‑TTGTTCAAAGGCGCTAAAGCCCagca                                                                                                                    >  1:217957/1‑239 (MQ=255)                                                                                                                                aCGCTACACTGATCCCTGCCGGGTTGAATGCGTCATTCCTGAAGGC‑GACAGAGGTGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGagca                                                                                                                    <  1:948900/238‑1 (MQ=255)                                                                                                                                 cGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAAAGTTTACCCTTCTTCGGCTCCCTCGAG‑ACTGCCCCGCCCGGGAATA‑TCTCTTTCCCGCTTGTCCAC‑AGTTTGTTGGTTTTGGTGATGGTCGCAATGGTGCAGCCCTTCA‑GGATTTTTGGGTG‑TTTCTCAAAGGCGATAAAAACGagccgc                                                                                                                  >  2:211421/1‑239 (MQ=255)                                                                                                                                   cTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGAATTCTCCGTTCTGCGGGCCCATGGGG‑AATGCCCCC‑CCGGGGAAATTTTCTGTCGCGAGGGTCCCA‑TGGTTGTTGTTCTTGGGGAGGCACGGAATGGTCCCGCCCGCCAGGGTTTTGTTGGTG‑TGTCTCAA‑GGCGGTAAAGcagcgaaaggg                                                                                                                >  2:306249/1‑229 (MQ=255)                                                                                                                                    ttcgcggATGCCTTCCGGTTTAGATGCGGCATTCCGGAAGGC‑GAATGATGAGCTGGTTGCCGCCGTTACAGCCACCGGGCGGG‑AAGATTTCCCCGTTCGGGGGCTACATGGAG‑CCTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGa                                                                                                               <  2:217957/233‑1 (MQ=255)                                                                                                                                    tACACTGATCCCTTCCGGGTTGAAAGCGGAATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTaa                                                                                                                                                                            <  1:1038332/180‑1 (MQ=255)                                                                                                                                      cacgtaTACCTTCCGGGTTAAAAGCGGCATTCCTGAAAGC‑GACAGATGACCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATg                                                                                                             <  1:367201/234‑1 (MQ=255)                                                                                                                                      cacTGCTCCCCTCCGGGTTAAAAGCGGCATCCCGTAAAGC‑AAGAGTGGGGCTGCTTGCCGGCGTTACAGCACCCTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATg                                                                                                             <  2:843320/239‑1 (MQ=255)                                                                                                                                      cacTGATCCCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGa                                                                                                                                                                                                                 <  1:378560/142‑1 (MQ=255)                                                                                                                                         tGATACCTTCCGGGTTGAAAGCGGCATGCCGGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCa                                                                                                                                                                                 <  1:1068091/170‑1 (MQ=255)                                                                                                                                            tACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGCG‑ACTGCCACG‑CCGGGGATATTCTCTGGCGCGATGGTGCAA‑tg                                                                                                                                                                                                      >  1:265452/1‑146 (MQ=255)                                                                                                                                             cccTTGCGGTTGGAAAGGGGCA‑TCCAGAAAGCGGAAAGATGAGCTTAATGCCGACGTTACAACACACTGGCGGG‑AAGATTTCTCCGTTCTGCGTCTACATGGCG‑ACTGCCAAG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTgg                                                                                                      <  2:1091746/238‑1 (MQ=255)                                                                                                                                              ctttCGGGGTTGAACGCGGATTTCCTGAAAGG‑GACAGTGGGGGTGTTGGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                                     <  2:129724/237‑1 (MQ=255)                                                                                                                                              ctgtCCGGTTTAAAAGCGGAATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                                     <  1:829951/236‑1 (MQ=255)                                                                                                                                              ccTTCGGGTTTGAAGGCGAGTTTCCTAAAGGC‑GACAGATGAGCTGATTGCCGCCGTTCCAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                                     <  1:686/239‑1 (MQ=255)                                                                                                                                              ccTTCCGTGTTGAAAGCGGCATTCCTGAAACC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACTCTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTGAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                                     <  1:93704/239‑1 (MQ=255)                                                                                                                                              ccTTCCGGGTTGAAAGCGGAATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                                     <  1:306249/239‑1 (MQ=255)                                                                                                                                               ctgcggtgTTGAACGAGGCATCCATTAAGGC‑GACAGATGTGTTGTTTCGCGCGGTTGCAGCCCCCTGGCGGG‑AAGATTTCCCGTTTGTGTGGCTAAATGGGG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa                                                                                                     <  2:70152/231‑1 (MQ=255)                                                                                                                                                   cGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGGTTTCACCCGTCCGCTGCTCAATGGGG‑AATGCCACG‑CCGGGGAAATTTTCTGGCCGGATTGTCCAA‑TGTTTGTTTATCTGGATGATGCACGTAATTGGCCCGCCGGTCT‑GGGTTTTTGGGTGTTTGCTC‑ATGGCGATAAAGACGCGCAGGGGGAGCAACAGGGAAAta                                                                                                >  2:6305/1‑239 (MQ=255)                                                                                                                                                          aaaGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCa                                                                                                                                                                                 <  1:154353/153‑1 (MQ=255)                                                                                                                                                          aaaGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACC‑CCGGGAATATTTTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTTGATGATGCACGTAATGGTCCAGCCGTTCA‑GGGTTTTTGGATG‑TTGCTCAACGGGGGTAAAACCGAGCAGCGGATGCAACTGGAAACTattattta                                                                                         >  1:2573/1‑237 (MQ=255)                                                                                                                                                          aaaGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAAATTTCCCCCCTTCGCGGGTCCAAGGGA‑AATGCCACC‑CCGGGAATATTTCCTGGGCGGACGGGCCAACTGTTTTTTGCTC‑GGCGGGTGCACGAAATGGGCCAGCCGGTCG‑GGGATTGGGGGGT‑TTGTTCAAAGGGGCTAAAACCCAGCCGCGGAGGCAAATGGAAAAAATTATTGa                                                                                         >  2:831286/1‑239 (MQ=255)                                                                                                                                                          aaaGCCGCTTTCTGGAAAGA‑GACAGGTGAGCTGTTGCCCGCGGTTACAGCACCCTGGCGGG‑AAGATTTCCCCGTTCGGCGGCTACATGGAG‑GCTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGa                                                                                         <  2:1085401/239‑1 (MQ=255)                                                                                                                                                                   ttCCGGAAGGG‑GACAGGTGGGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCCCCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATga                                                                                <  1:211421/239‑1 (MQ=255)                                                                                                                                                                                            aTTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGc                                                        <  1:15743/239‑1 (MQ=255)                                                                                                                                                                                                      ttACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGa                                              <  1:185596/239‑1 (MQ=255)                                                                                                                                                                                                                 tGGCGGG‑AAGATTTCACCGTTCTGCGTCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGc                                                                                      <  1:921675/188‑1 (MQ=255)                                                                                                                                                                                                                            aTTTCCCCGTTCGGCGGCTCCATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGa                                              <  1:42284/218‑1 (MQ=255)                                                                                                                                                                                                                                   cGTTTTGGGGCTACATGGGG‑GTGGCCACG‑CCGGGAATTTTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGtt                  <  2:829308/239‑1 (MQ=255)                                                                                                                                                                                                                                   cGTTCTGCGGCTACATGTAG‑CCTGCCACG‑CCGGGAATATTCTGTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCCCGTAAGTGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGtt                  <  1:948163/239‑1 (MQ=255)                                                                                                                                                                                                                                   cGTTCTGCGGCTACATGGGG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGtt                  <  1:466430/239‑1 (MQ=255)                                                                                                                                                                                                                                   cGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGtt                  <  1:605695/239‑1 (MQ=255)                                                                                                                                                                                                                                       cTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCCGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGt              >  1:563373/1‑239 (MQ=255)                                                                                                                                                                                                                                          cGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGATTATGAAGAATTTAGCGAGTTCGAACCCCCTGAAATCGGAATGATTGAAACTTTACCCCCGCTGCGTTTGGtttt           >  2:1152/1‑238 (MQ=255)                                                                                                                                                                                                                                                   ggAG‑ACTGCCACG‑CCGGGAATTTTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGc  <  1:445251/239‑1 (MQ=255)                                                                                                                                                                                                                                                   ggAG‑ACTGCCACG‑CCGGGAATATTCTCGGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGc  <  1:85173/239‑1 (MQ=255)                                                                                                                                                                                                                                                   ggAG‑ACGGCCACG‑CCGGGAATATTCTCTGGCTCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGc  <  1:28870/239‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  AGATGGAAGCGGTTGGGCGTTATTTAGGGC‑GTATGCA‑CCAGACGGGGCGCAAACAGCTTTTTAT‑CCATCGCCCGACCATCGGTTTGAACGAATATCTC‑ATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGC‑GACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑AAGATTTCACCGTTCTGCGGCTACATGGAG‑ACTGCCACG‑CCGGGAATATTCTCTGGCGCGATGGTCCAA‑TGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCA‑GGATTTGTGGATG‑TTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGC  >  NC_000913/4042776‑4043250

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