Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,012,247 T→C L190L (TTA→TTG)  metR ← DNA‑binding transcriptional dual regulator MetR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,012,2470TC95.7% 146.2 / ‑1.9 46L190L (TTA→TTG) metRDNA‑binding transcriptional dual regulator MetR
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (19/25);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

CGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGG  >  NC_000913/4012111‑4012383
                                                                                                                                        |                                                                                                                                        
cgGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAAt                                                                                                                                       >  2:988051/1‑140 (MQ=255)
        cccATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGtctc                                                                                                                               >  2:1151172/1‑140 (MQ=255)
            tCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTC‑‑‑‑ACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCga                                                                                                                       >  2:1001108/1‑140 (MQ=255)
               gcgcgGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCg                                                                                                                        <  2:253646/140‑1 (MQ=255)
               gcgcgGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCg                                                                                                                        <  2:807252/140‑1 (MQ=255)
                  cgGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCgaga                                                                                                                     >  2:336188/1‑140 (MQ=255)
                                    acaaTAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTcgc                                                                                                   >  1:944582/1‑140 (MQ=255)
                                      aaTAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTC‑‑‑‑ACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTc                                                                                             <  1:1001108/140‑1 (MQ=255)
                                       aTAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGCTCTTCCGGTGTAATTCGCGtt                                                                                                <  2:829236/140‑1 (MQ=255)
                                         aaGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGtttt                                                                                              >  1:201298/1‑140 (MQ=255)
                                         aaGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGtttt                                                                                              >  2:58606/1‑140 (MQ=255)
                                           ggTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTcg                                                                                            <  1:1151172/140‑1 (MQ=255)
                                               ttATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTcgccgc                                                                                        >  2:36193/1‑140 (MQ=255)
                                               ttATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTcgccgc                                                                                        >  2:945877/1‑140 (MQ=255)
                                                tATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCgccgcc                                                                                       <  2:322576/140‑1 (MQ=255)
                                                 aTCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCa                                                                                      >  2:1119711/1‑140 (MQ=255)
                                                       gCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGAt                                                                                >  2:436126/1‑140 (MQ=255)
                                                       gCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGAt                                                                                >  1:1083920/1‑140 (MQ=255)
                                                          tttCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGt                                                                             <  2:944582/140‑1 (MQ=255)
                                                            tCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCa                                                                           <  1:998761/140‑1 (MQ=255)
                                                               gTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGt                                                                        <  1:436126/140‑1 (MQ=255)
                                                                tGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTg                                                                       >  1:786869/1‑140 (MQ=255)
                                                                tGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTg                                                                       >  2:22991/1‑140 (MQ=255)
                                                                 gACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGc                                                                      >  2:639070/1‑140 (MQ=255)
                                                                     ggCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAac                                                                  <  2:93307/140‑1 (MQ=255)
                                                                        tGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACAcc                                                               <  1:916273/140‑1 (MQ=255)
                                                                             ccTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTTCGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACg                                                          <  2:630965/140‑1 (MQ=255)
                                                                             ccTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACg                                                          <  2:201298/140‑1 (MQ=255)
                                                                             ccTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACg                                                          <  1:639070/140‑1 (MQ=255)
                                                                             ccTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACg                                                          <  1:284616/140‑1 (MQ=255)
                                                                                 ccGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACt                                                      >  1:1063029/1‑140 (MQ=255)
                                                                                    gCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCa                                                   <  1:58606/140‑1 (MQ=255)
                                                                                        aagaaAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCC‑CCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTc                                              >  2:440200/1‑140 (MQ=255)
                                                                                          gaaAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCg                                             <  2:1083920/140‑1 (MQ=255)
                                                                                                 ccgccAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCg                                      <  1:742795/140‑1 (MQ=255)
                                                                                                  cgccAGACATCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCgg                                     >  1:617949/1‑140 (MQ=255)
                                                                                                          aTCCAGTCGACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCC‑CCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATg                            <  1:440200/140‑1 (MQ=255)
                                                                                                                  gACTACGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCa                     >  2:552294/1‑140 (MQ=255)
                                                                                                                       cGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTgcg                <  2:686686/140‑1 (MQ=255)
                                                                                                                       cGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTgcg                <  2:786869/140‑1 (MQ=255)
                                                                                                                       cGCTGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTgcg                <  1:36193/140‑1 (MQ=255)
                                                                                                                          tGCACCGGATAAATCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCgg             <  1:574568/140‑1 (MQ=255)
                                                                                                                                 gATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGGCGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAata      >  2:668233/1‑140 (MQ=255)
                                                                                                                                     aaTTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCgg  <  1:318740/140‑1 (MQ=255)
                                                                                                                                     aaTCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCGGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCgg  <  1:336188/140‑1 (MQ=255)
                                                                                                                                     aaTCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCgg  <  2:617949/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
CGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGG  >  NC_000913/4012111‑4012383

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