Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,619,314 T→C T283A (ACG→GCG)  uraA ← uracil permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,619,3140TC91.7% 69.9 / ‑1.5 24T283A (ACG→GCG) uraAuracil permease
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (10/12);  total (11/13)
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.27e-01

CGGCAGCCAGTTTACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTGCAGACCTGGATCGCGCAGCAGATCTTTTTTGACGATATTAGCCGTTACTACCAGGTG  >  NC_000913/2619188‑2619428
                                                                                                                              |                                                                                                                  
cGGCAGCCAGTTTACCGACGCAGGAAAGCAGGACAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAaagaa                                                                                                        >  2:363241/1‑139 (MQ=255)
                                        aaTCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGaa                                                                <  1:363241/139‑1 (MQ=255)
                                        aaTCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGaa                                                                <  1:222720/139‑1 (MQ=255)
                                             ccgccccgccGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTg                                                           <  2:533471/139‑1 (MQ=255)
                                             cccccccgccGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTg                                                           <  2:505966/136‑1 (MQ=255)
                                                                     tGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGaa                                                                <  2:56182/110‑1 (MQ=255)
                                                                     tGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGaa                                                                >  1:56182/1‑110 (MQ=255)
                                                                     tGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCACAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTGCAGACCTGGATCGCGCAGCAGATc                                   >  1:181990/1‑139 (MQ=255)
                                                                      gTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTGCAGACCTGGATCGCGCAGCAGATCt                                  >  1:142232/1‑139 (MQ=255)
                                                                              gggTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTGCAGACCTGGATCGCGCAGCAGATCTTTTTTGAc                          >  2:416698/1‑139 (MQ=255)
                                                                               ggTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTg                                                           >  1:17030/1‑105 (MQ=255)
                                                                               ggTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTg                                                           <  2:17030/105‑1 (MQ=255)
                                                                                   aTCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGaa                                                                <  2:451314/96‑1 (MQ=255)
                                                                                   aTCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGaa                                                                >  1:451314/1‑96 (MQ=255)
                                                                                   aTCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTg                                                           <  1:509777/101‑1 (MQ=255)
                                                                                   aTCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTg                                                           >  2:509777/1‑101 (MQ=255)
                                                                                                    gTTTTCTCCGTAAGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTGCAGACCTGGATCGCGCAGCAGATCTTTTTTGACGATATTAGCCGTTACTACCAgg    <  2:142232/139‑1 (MQ=255)
                                                                                                      tttCTCCGTAAGTAGTATTTGGCACAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTGCAGACCTGGATCGCGCAGCAGATCTTTTTTGACGATATTAGCCGTTACTACCAGGTg  <  2:181990/139‑1 (MQ=255)
                                                                                                            cGTAAGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGaa                                                                <  1:571639/71‑1 (MQ=255)
                                                                                                            cGTAAGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGaa                                                                >  2:571639/1‑71 (MQ=255)
                                                                                                             gTACGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTg                                                           >  2:394760/1‑75 (MQ=255)
                                                                                                             gTAAGTAGTATTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTg                                                           <  1:394760/75‑1 (MQ=255)
                                                                                                                     taTTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTg                                                           >  2:250533/1‑67 (MQ=255)
                                                                                                                     taTTTGGCGCAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTg                                                           <  1:250533/67‑1 (MQ=255)
                                                                                                                              |                                                                                                                  
CGGCAGCCAGTTTACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTGCAGACCTGGATCGCGCAGCAGATCTTTTTTGACGATATTAGCCGTTACTACCAGGTG  >  NC_000913/2619188‑2619428

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: