Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,520,504:1	+T	20.0%	intergenic (‑180/+167)	yjhU ← / ← yjhF	putative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,520,504	1	.	T	20.0%	123.0 / 41.0	45	intergenic (‑180/+167)	yjhU/yjhF	putative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand):  ref base . (19/17);  new base T (5/4);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTG  >  NC_000913/4520362‑4520650                                                                                                                                                 |                                                                                                                                                   ttAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCgg                                                                                                                                              >  1:255653/1‑149 (MQ=255)  tAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGc                                                                                                                                             <  2:61243/149‑1 (MQ=255)   aaCCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACTAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCt                                                                                                                                            <  4:462401/149‑1 (MQ=255)      ccTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTa                                                                                                                                         <  1:234258/149‑1 (MQ=255)                cGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACt                                                                                                                               >  3:155860/1‑149 (MQ=255)                       ttGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCttt                                                                                                                                          >  1:290833/1‑131 (MQ=255)                       ttGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCttt                                                                                                                                          <  2:290833/131‑1 (MQ=255)                               gATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAccc                                                                                                                   <  2:98036/146‑1 (MQ=255)                               gATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAccc                                                                                                                   >  1:98036/1‑146 (MQ=255)                                    ccAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGaa                                                                                                                          <  2:6406/134‑1 (MQ=255)                                    ccAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGaa                                                                                                                          >  1:6406/1‑134 (MQ=255)                                     cAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCAt                                                                                                          >  3:802543/1‑149 (MQ=255)                                             tGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAg                                                                                                  <  4:155860/149‑1 (MQ=255)                                               aTTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGt                                                                                                <  2:255653/149‑1 (MQ=255)                                                                     ttttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCt                                                                                                                                            >  2:362300/1‑84 (MQ=25)                                                                     ttttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCt                                                                                                                                            <  1:362300/84‑1 (MQ=25)                                                                     ttttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAg                                                                                                  >  1:64784/1‑125 (MQ=255)                                                                     ttttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAg                                                                                                  <  2:64784/125‑1 (MQ=255)                                                                            tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                                                                   >  3:323580/1‑149 (MQ=255)                                                                                    aCCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCa                                                           <  4:802543/149‑1 (MQ=255)                                                                                        gCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTa                                                                                                                                 >  2:757945/1‑75 (MQ=255)                                                                                        gCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTa                                                                                                                                 <  1:757945/75‑1 (MQ=255)                                                                                                tttGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt                                                                 >  2:416591/1‑131 (MQ=255)                                                                                                tttGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt                                                                 <  1:416591/131‑1 (MQ=255)                                                                                                    aGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGc                                           >  3:153074/1‑149 (MQ=255)                                                                                                          cGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtcac                                     <  2:135819/149‑1 (MQ=255)                                                                                                            cccGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAg                                                                                                  >  4:524939/1‑86 (MQ=255)                                                                                                            cccGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAg                                                                                                  <  3:524939/86‑1 (MQ=255)                                                                                                                       ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGc                                           <  1:98858/131‑1 (MQ=255)                                                                                                                       ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGc                                           >  2:98858/1‑131 (MQ=255)                                                                                                                         ccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTaaa                             >  1:114878/1‑143 (MQ=255)                                                                                                                         ccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTaaa                             <  2:114878/143‑1 (MQ=255)                                                                                                                         ccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATgg                       >  4:769586/1‑149 (MQ=255)                                                                                                                               tGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTAGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACa                 >  1:237630/1‑149 (MQ=255)                                                                                                                                  aaaGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAA‑TATGGGGGGCaagga             >  2:141024/1‑145 (MQ=255)                                                                                                                                   aaGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGa             >  3:562306/1‑149 (MQ=255)                                                                                                                                    aGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATAGACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaa            >  3:696929/1‑149 (MQ=255)                                                                                                                                    aGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaa            >  3:35847/1‑149 (MQ=255)                                                                                                                                    aGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaa            >  1:138476/1‑149 (MQ=255)                                                                                                                                    aGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTCCAATATCGCCCCTTTTTGATCGCCAGCGATTTGTCACCGTGCAGAGGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaa            >  1:221578/1‑149 (MQ=255)                                                                                                                                       cGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaacaa         <  3:769586/149‑1 (MQ=255)                                                                                                                                          agaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAg      <  4:323580/149‑1 (MQ=255)                                                                                                                                              aTG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAgg     <  1:473235/146‑1 (MQ=255)                                                                                                                                              aTG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAgg     >  2:473235/1‑146 (MQ=255)                                                                                                                                              aTG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTg  <  2:221578/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   TTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTT‑CCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTG  >  NC_000913/4520362‑4520650

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