BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,446,579	T→C	T299A (ACC→GCC)	aroC ←	chorismate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,446,579	0	T	C	97.5%	140.6 / ‑3.0	40	coding (895/1086 nt)	aroC	chorismate synthase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (22/17); minor base . (1/0); total (23/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.92e-01

TGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCA  >  NC_000913/2446446‑2446713

tGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCt                                                                                                                                  >  1:537562/1‑140 (MQ=255)

 gATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTg                                                                                                                                 >  2:513801/1‑140 (MQ=255)

 gATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTg                                                                                                                                 >  2:171941/1‑140 (MQ=255)

   tCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGa                                                                                                                               >  2:911910/1‑140 (MQ=255)

           aaCGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGtt                                                                                                                       <  2:537562/140‑1 (MQ=255)

               aTCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGCGGttt                                                                                                                          >  2:573191/1‑131 (MQ=255)

               aTCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGCGGttt                                                                                                                          <  1:573191/133‑3 (MQ=255)

               aTCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAg                                                                                                                   >  1:798821/1‑140 (MQ=255)

                tCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAgc                                                                                                                  <  1:71725/140‑1 (MQ=255)

                 cGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAgcg                                                                                                                 >  1:1040439/1‑140 (MQ=255)

                          cGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGttt                                                                                                                      >  2:451024/1‑126 (MQ=255)

                          cGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGttt                                                                                                                      <  1:451024/126‑1 (MQ=255)

                          cGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAg                                                                                                                   <  2:74664/129‑1 (MQ=255)

                          cGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAg                                                                                                                   >  1:74664/1‑129 (MQ=255)

                          cGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATggg                                                                                                        <  2:798821/140‑1 (MQ=255)

                               cTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATg                                                                                                   <  1:791772/140‑1 (MQ=255)

                                        gCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTtgctgc                                                                                          <  1:513801/140‑1 (MQ=255)

                                         cACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGcc                                                                                         <  2:885734/140‑1 (MQ=255)

                                                   aTCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGAt                                                                               <  1:339751/140‑1 (MQ=255)

                                                               ggATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAAt                                                                   >  2:916783/1‑140 (MQ=255)

                                                                  tCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATgcc                                                                <  1:171941/140‑1 (MQ=255)

                                                                    gTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATgccgc                                                              >  2:94647/1‑140 (MQ=255)

                                                                     tGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATgccgcc                                                             >  2:590352/1‑140 (MQ=255)

                                                                       aCGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGccgcccg                                                           <  1:916783/140‑1 (MQ=255)

                                                                           ccTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATg                                                       >  2:145026/1‑140 (MQ=255)

                                                                           ccTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATg                                                       >  1:868590/1‑140 (MQ=255)

                                                                                   gATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGctct                                               <  1:590352/140‑1 (MQ=255)

                                                                                       atcTCAACTTCTTCGCCAAAGCGGTTAATGGTCCGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGaa                                           <  1:1127773/140‑1 (MQ=255)

                                                                                         ctcAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAAc                                         >  2:344160/1‑140 (MQ=255)

                                                                                         ctcAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAAc                                         >  1:1109387/1‑140 (MQ=255)

                                                                                              cttcttCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCt                                    <  2:868590/140‑1 (MQ=255)

                                                                                              cttcttCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCt                                    >  1:598596/1‑140 (MQ=255)

                                                                                                    cGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTg                              >  2:799576/1‑140 (MQ=255)

                                                                                                    cGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTg                              >  1:238811/1‑140 (MQ=255)

                                                                                                         aaGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTc                         >  2:1058872/1‑140 (MQ=255)

                                                                                                             ggTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATcg                     <  1:750265/140‑1 (MQ=255)

                                                                                                                 aaTGGTACGACCCGGCACGG‑CATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGtt                >  1:104179/1‑140 (MQ=255)

                                                                                                                           accGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGcc       >  2:129918/2‑140 (MQ=255)

                                                                                                                             cggcacggcaATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCcgc     >  1:1136793/1‑140 (MQ=255)

                                                                                                                              ggCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCcgcg    <  2:312372/140‑1 (MQ=255)

                                                                                                                              ggCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCcgcg    <  2:1109387/140‑1 (MQ=255)

                                                                                                                                cACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCa  >  2:726874/1‑140 (MQ=255)

TGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCA  >  NC_000913/2446446‑2446713

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`