Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 508,702 T→C intergenic (+143/+173) ybaQ → / ← copA putative DNA‑binding transcriptional regulator/copper transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913508,7020TC96.0% 44.9 / ‑4.1 25intergenic (+143/+173)ybaQ/copAputative DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (15/9);  minor base A (1/0);  total (16/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAA‑‑CTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGG  >  NC_000913/508608‑508798
                                                                                                |                                                                                                 
catcaggcatctgCGCACGACTGCCGGATGCGACGCTAACGCGTCTTATCAGACCTACAAATTCTCACGCC‑‑ATCCGTAGGCCGGATAAGGCGTTCAcgcc                                                                                              >  1:522051/1‑100 (MQ=14)
                                                           caggtttGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGGCGCGAACGC‑‑CTTATCCGGCCTAcg                                 >  1:537006/6‑99 (MQ=11)
                                                                     ttcattttgTAGGCAGAATAAGACGTTCTCGCAGCATCCGGCATCCGTGCACCATTTCCTGATGCGACGCTTACGCGTCTTATCAGCCCTACAATCCtg                           <  2:2633486/91‑1 (MQ=0)
                                                                         tctcGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCt                      <  1:474169/97‑1 (MQ=2)
                                                                           tcGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGa                    >  1:1949835/2‑100 (MQ=17)
                                                                           tcGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGCGCCTGa                    >  2:2932813/2‑100 (MQ=2)
                                                                            cGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGaa                   >  1:1412672/1‑100 (MQ=18)
                                                                             gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAAc                  <  2:558763/100‑1 (MQ=18)
                                                                             gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAAc                  >  1:2132695/1‑100 (MQ=12)
                                                                              tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAAcc                 <  2:1183507/100‑1 (MQ=18)
                                                                               atgCCGGATGCGGCGTGAACGCCTTATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCg                >  1:28932/3‑100 (MQ=11)
                                                                                ggCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGt               >  2:146946/1‑100 (MQ=18)
                                                                                atctggATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGt               <  2:1485768/96‑1 (MQ=255)
                                                                                 gCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTa              <  1:1989125/100‑1 (MQ=18)
                                                                                  ccGGATAAGGCGTTCACGCCGCATCCGGCATTCGAGCACCATTGCCAGATGCGACGCTTACGCGTCCTACCAGGCCTACAAACCTGTGCCTGAACCGTa              >  1:973830/1‑99 (MQ=11)
                                                                                  ccGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAg             >  2:789728/1‑100 (MQ=18)
                                                                                    ggATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGTCCTACAAACCTGTGCCTGAACCGTAGGc           >  2:2996528/1‑100 (MQ=9)
                                                                                    ggATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAgat           >  1:1485780/1‑98 (MQ=16)
                                                                                    ggATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGc           >  1:3025963/1‑100 (MQ=17)
                                                                                       tAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCgg        <  1:1781705/100‑1 (MQ=18)
                                                                                         aGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGAt      >  1:2680892/1‑100 (MQ=18)
                                                                                           gCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATaa    >  2:331977/1‑100 (MQ=25)
                                                                                            cGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAg   <  2:2187382/100‑1 (MQ=18)
                                                                                            cGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAg   >  2:2505033/1‑100 (MQ=25)
                                                                                             gTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAgg  <  1:1818915/100‑1 (MQ=25)
                                                                                                |                                                                                                 
CATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAA‑‑CTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGG  >  NC_000913/508608‑508798

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