Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,635,048	A→C	F88C (TTT→TGT)	rob ←	right oriC‑binding transcriptional activator, AraC family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,635,048	0	A	C	94.2%	163.0 / ‑5.8	52	F88C (TTT→TGT)	rob	right oriC‑binding transcriptional activator, AraC family
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base C (23/26);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTC  >  NC_000913/4634900‑4635192                                                                                                                                                       |                                                                                                                                                 cggcgTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACAt                                                                                                                                                <  2:579522/151‑1 (MQ=255)   gcgTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGt                                                                                                                                              <  1:534628/151‑1 (MQ=255)            ccAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTgag                                                                                                                                      <  1:493257/150‑1 (MQ=255)            ccAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTgag                                                                                                                                      <  1:216972/150‑1 (MQ=255)              aGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGGATACCAAAAGGGCTCCCTTCAGGAGAACGGCGGGAAAGTGGAGGAGGTTTGGGCAAA‑TGGTTC‑TTGAATGGGCGGGGGAAAGTCTTTTGggag                                                                                                                                    >  1:205825/1‑150 (MQ=255)                ggTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAAGCGCTCCATTCAGGAGAACGGCGGGGAAAAGGAGGG‑GGCTGGGCAAAATGCCTT‑TTTTAAGCGCGGGGACAAGTCTGttttgagggg                                                                                                                                 >  1:538120/1‑143 (MQ=255)                  tGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAAGCGCACCATTCCGGAGAACGGCGGGTAAAAGCCGGA‑GGCTGGGGAAAACGCTTC‑TTTTATTCCCGGGGACATGTTTGTTgagaagc                                                                                                                                  >  1:590455/1‑145 (MQ=255)                   gACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGaa                                                                                                                               <  1:725094/150‑1 (MQ=255)                   gACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAg                                                                                                                              <  1:831764/151‑1 (MQ=255)                          ttGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTAt                                                                                                                        <  1:209418/150‑1 (MQ=255)                           tgtgCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGGAAAAAGGAGGA‑GGCTGGGGCAAATCCTTCTTTTAATG‑GGGGGGACAAGTCTGTTGGGAGTCGAAACgggaat                                                                                                                       >  1:856952/1‑146 (MQ=255)                             tgCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCCTTCAGGAGAACGGCGGGGAAAGGGAGGG‑GGCTGGGCAAACTGCTTC‑TTTAAAGCGCGGGGACATGTCTGTTGGGAGTCGAAGCGGGAtt                                                                                                                       >  1:664754/1‑148 (MQ=255)                              gCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAAGCGCTCCCTTCAGGAGAACGGCGGGAAAAGGGAGGG‑GGGTGGGGGAAAAGCTTCTTTTAATG‑GCGGGGACAAGTCTGTTTGGAGTCGGAACggggt                                                                                                                        >  1:559688/1‑143 (MQ=255)                                tctGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAAAAACAAAGGCGCTCCATTCAAGGGGACGGCGGGGAAGGGCGAGA‑GGGTGGGGAAAATCGTTC‑TTGTGAGTGGGGGGTAATGTCTGGTTAGAGGCGGAGCgggaatttagc                                                                                                                  >  2:306517/1‑141 (MQ=255)                                 ctGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAgc                                                                                                                  <  2:798048/149‑1 (MQ=255)                                        aGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGGGAACGGCGGGAAAAAGGAGGG‑GGCTGGGGCAAAAGCTTT‑TTTGAAGGGCGGGGAAATGTCTGTTGGGAGGCGGAGGGGGATTTTAgcggggtggtc                                                                                                         >  1:296817/1‑143 (MQ=255)                                        aGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTcc                                                                                                         <  2:90241/151‑1 (MQ=255)                                        aGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTc                                                                                                          <  1:562757/150‑1 (MQ=255)                                        aGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAAGCGCTCCATTCAGGAGAACGGCGGGAAAAAGGAGGG‑GGCTGGGGAAAAAGCCTC‑TTTAAAGCGCGGGGACATGTCTGTTGAGAGGCGAAACGGGATTTCCGCGCGATGTcc                                                                                                         >  1:805263/1‑151 (MQ=255)                                                  cccAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGAc                                                                                               <  2:290658/151‑1 (MQ=255)                                                  cccAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGAc                                                                                               <  2:529413/151‑1 (MQ=255)                                                         gTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGt                                                                                        <  2:707590/151‑1 (MQ=255)                                                           aGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCTGACGCGCAGTCa                                                                                      <  1:107263/151‑1 (MQ=255)                                                             cggcggGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGGGGTGAAGGGCAGGA‑GTCTGGGCAAACCGGTTC‑TTGAAAGCGCGGGGTCATGTTTGTTGGGAGTCGAAGCGGGATTGGAGCGCGGAGTCCAGGATCGGACGCCCAGTCAgg                                                                                    >  2:447749/1‑151 (MQ=255)                                                                   gcgAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTa                                                                                >  2:329664/1‑149 (MQ=255)                                                                     gAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTg                                                                             >  1:299613/1‑150 (MQ=255)                                                                     gAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGc                                                                            <  1:736077/151‑1 (MQ=255)                                                                     gAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGc                                                                            >  1:421278/1‑151 (MQ=255)                                                                      aaTACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCg                                                                           >  2:71384/1‑151 (MQ=255)                                                                      aaTACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCg                                                                           >  2:528678/1‑151 (MQ=255)                                                                        tACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGGGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGAc                                                                         >  2:340903/1‑151 (MQ=255)                                                                              aGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGAccgccg                                                                    >  1:433732/1‑150 (MQ=255)                                                                                        ttCAGGAGAACGGCGGTAAAGTGCAGGA‑GGCTGGGCCAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACaa                                                          >  2:785614/1‑150 (MQ=255)                                                                                          cAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACaaa                                                         >  2:409025/1‑149 (MQ=255)                                                                                           aGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAGGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAac                                                        <  2:62755/149‑1 (MQ=255)                                                                                            ggAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAacgac                                                     <  2:78818/151‑1 (MQ=255)                                                                                            ggAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAacg                                                       <  1:375589/149‑1 (MQ=255)                                                                                              agaACGGGGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAGGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAAcgacga                                                   <  2:862024/151‑1 (MQ=255)                                                                                                 aCGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTGGAAGCGGTATTGCAGCGCGTTGTCCAGAATGGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAg                                                  <  2:549900/149‑1 (MQ=255)                                                                                                    gcggTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGGGGTTACATGTCTGTTGAGAGTCGAAGCGTTATTGAAGCGCGATGTCCAGAATCGGACGCGCAGTAAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCAc                                               <  2:764512/149‑1 (MQ=255)                                                                                                     cggTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGa                                             >  1:20892/1‑150 (MQ=255)                                                                                                           aGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATAcg                                      <  1:56681/151‑1 (MQ=255)                                                                                                                  gA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGCCGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAAt                               >  1:510413/1‑151 (MQ=255)                                                                                                                      ccgatctaaaCTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATa                              <  2:310487/142‑1 (MQ=255)                                                                                                                             aaaCTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATa                              >  1:310521/1‑142 (MQ=255)                                                                                                                                cTGCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCAAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGt                   <  1:171160/150‑1 (MQ=255)                                                                                                                                  gCTTC‑TTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGAc                >  1:520737/1‑151 (MQ=255)                                                                                                                                      c‑ttGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCt            >  1:332642/1‑151 (MQ=255)                                                                                                                                        ttGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCtt           >  1:73633/1‑151 (MQ=255)                                                                                                                                                gcgGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACAtt   >  2:617908/1‑151 (MQ=255)                                                                                                                                                 cgGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACAt    <  1:481199/149‑1 (MQ=255)                                                                                                                                                  gggTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATtc  <  1:251554/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                 CGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGA‑GTCTGGGCAAACTGCTTC‑TTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTC  >  NC_000913/4634900‑4635192

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