Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,726,876 C→T P73L (CCT→CTT)  nemA → chromate reductase, quinone reductase, FMN‑linked; N‑Ethylmaleimide reductase; old yellow enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,726,8760CT97.9% 133.8 / NA 48P73L (CCT→CTT) nemAchromate reductase, quinone reductase, FMN‑linked; N‑Ethylmaleimide reductase; old yellow enzyme
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (28/18);  minor base G (0/1);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGCAGTATTGAACCGGGTGACATTCCTACCCCGTTGATGGCGGAATACTATCGCCAACGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGC  >  NC_000913/1726744‑1727014
                                                                                                                                    |                                                                                                                                          
tGCGCAGTATTGAACCGGGTGACATTCCTACCCCGTTGATGGCGGAATACTATCGCCAACGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCCTGGCATGCac                                                                                                                                  >  2:695853/1‑142 (MQ=255)
  cgcAGTATTGAACCGGGTGACATTCCTACCCCGTTGATGGCGGAATACTATCGCCAACGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCCACACGATCTGCAGGTGCGCTTGGCCTCCGTAg                                                                                                                                >  1:213794/1‑143 (MQ=255)
         ttGAACCGGGTGACATTCCTACCCCGTTGATGGCGGAATACTATCGCCAACGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGCATATGCAGGTGCGCTAGGCATCCATAGTCCGGAg                                                                                                                         >  2:1006594/1‑143 (MQ=255)
                       aTTCCTACCCCGTTGATGGCGGAATACTATCGCCAACGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTg                                                                                                                                          <  1:181426/112‑1 (MQ=255)
                       aTTCCTACCCCGTTGATGGCGGAATACTATCGCCAACGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTg                                                                                                                                          >  2:181426/1‑112 (MQ=255)
                             aCCCCGTTGATGGCGGAATACTATCGCCAACGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGGGCGCTTGGCATCCATAGGCCGGAGCACGTGGCCGCATGGaaaaa                                                                                                     >  1:309324/1‑143 (MQ=255)
                               cccGTTGATGGCGGAATACTATCGCCAACGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGaaaaaaa                                                                                                   >  2:270746/1‑143 (MQ=255)
                               cccGGTGATGGCGGAACACTATCGCCAACTTGCGAGTGCCGGTTTTATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGTAAAAGGATATGCAGGTGCGCTTGGCATCCGTAGTCCGGAGCAAATTGCCGCATGGaaaaaaa                                                                                                   <  1:995600/143‑1 (MQ=255)
                                      aTGGCGGAATACTATCGCCAACGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGc                                                                                            >  1:682522/1‑143 (MQ=255)
                                                 cTATCGCCAACGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGAATCCATAGGCCGGAGCAAATTGCCGCATGGAAAAACAACACAGCTGGTGTTCATg                                                                                 >  2:631835/1‑143 (MQ=255)
                                                  tgtCGCCAACGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGc                                                                                <  2:309432/141‑1 (MQ=255)
                                                  tATCGCCAACGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATGCATAGTCCGGAGCAAATTGCGGCAGGGAGGAAAATCACCGCTGGCGGTCGTGc                                                                                >  2:1056810/1‑143 (MQ=255)
                                                  tATCGCCAACGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGc                                                                                >  1:1149042/1‑143 (MQ=255)
                                                  tATCGCCAACGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGACCCAAAGGCAACAAAGTTAGAAGTCACCGCTGGCCTTCACGc                                                                                >  1:492301/1‑143 (MQ=255)
                                                  tATCGCCAACGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCCCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGc                                                                                >  1:992740/1‑143 (MQ=255)
                                                   aTCGCCAACGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGGTCATGCt                                                                               >  2:679411/1‑143 (MQ=255)
                                                         aaCGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAAt                                                                         <  1:430590/143‑1 (MQ=255)
                                                         aaCGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCGGGCGCTCATGCGGaacaa                                                                         >  1:729356/1‑140 (MQ=255)
                                                          aCGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATg                                                                        <  1:635955/143‑1 (MQ=255)
                                                           cGTGACAGTGCCGGTTGGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATgg                                                                       <  1:340685/143‑1 (MQ=255)
                                                             tGCCTGTGCCGGGTTGATTATTACTGAAGCCACGCAAATTTCTGCTCAGGCAACAGGATTTGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATAGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTc                                                                     <  2:583552/143‑1 (MQ=255)
                                                              gCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCa                                                                    <  1:624601/143‑1 (MQ=255)
                                                                   tGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCAGAGGAGATAATAGTTTGGTAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATgg                                                               >  1:862591/1‑143 (MQ=255)
                                                                     ccGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCAGGGAAAAAAATCACGGCTGGCGTTCATGCTGAAAATGGTCATATGGcc                                                             <  1:955684/143‑1 (MQ=255)
                                                                     ccGGTTTGATTATTAGTGAAGCCACGCAAAATACTACCCATTCAAAAAACTGACCCGGTGCGCTTTTCATCCATAGTCCGGAGCAAATGGCCGCAGGGAAAAAAATCAACGCTGTCGTGCATGCCGAACAGGGTCAGATGGcc                                                             >  2:1101140/1‑143 (MQ=255)
                                                                        gTTTGATTATGAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGGTATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTg                                                          <  1:44267/143‑1 (MQ=255)
                                                                           tGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAACCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAg                                                       >  2:817902/1‑143 (MQ=255)
                                                                                 ttaGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTGCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCAGGCTGAAAAGGGTCATAGGGCCGTGCAGCTGTgg                                                 >  2:612275/1‑143 (MQ=255)
                                                                                 ttaGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTgg                                                 <  1:838049/143‑1 (MQ=255)
                                                                                       aaGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCAGGAGCAAATTGCCGCATGGAAAAAAATCACCGCACGCGCTCGCGCGGACAATGATCATATGGCCGTGCAGCTGTGGCACAcc                                           >  2:1021532/1‑143 (MQ=255)
                                                                                        aGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATAGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCg                                          >  1:708989/1‑143 (MQ=255)
                                                                                          ccACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGa                                        >  1:1013894/1‑143 (MQ=255)
                                                                                            aCGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACg                                      >  1:636323/1‑143 (MQ=255)
                                                                                             cGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGCTCATATGGCCGTGCAGCTGTGGCACACCGGACGc                                     >  1:365206/1‑143 (MQ=255)
                                                                                                 aaTTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCACATTGCCGCATCGAAAAAACTCACCGCCGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCAttt                                 <  2:106543/143‑1 (MQ=255)
                                                                                                   ttccTGCCCAGGCAAAAGGATATGCACGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAACATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTtct                               <  2:909326/140‑1 (MQ=255)
                                                                                                    ttCTGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTtctc                              <  2:682522/143‑1 (MQ=255)
                                                                                                       tGCCCAGGCAAAAGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATGCCTCACg                           >  2:341727/1‑143 (MQ=255)
                                                                                                             ggCAGAAGGATAGACAGGGGCGCGTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCGGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGcc                     <  2:459679/143‑1 (MQ=255)
                                                                                                               cAAAAGGATATGCAGGTGCGCTTGGCACCCACAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTg                   <  2:992740/143‑1 (MQ=255)
                                                                                                                aaaaGGATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGc                  >  2:612288/1‑143 (MQ=255)
                                                                                                                  aaGGATATGCAGGTGCGCTTGGCCTCCATAGTCAGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGCTCATATGGACGTGCAGCAGTGGCACACCGGACGCATATCCCACGCCAGCCCGccc                >  2:158931/1‑141 (MQ=255)
                                                                                                                  aaGGATATGCAGGTGCGCTTGGCATCCAGAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCaa                <  1:242761/143‑1 (MQ=255)
                                                                                                                    ggATATGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATAGCCGTGCAGCTGTGGCACACCGCAAGCCTATCTCACGCAAGCGTGCAAca              >  2:564865/1‑142 (MQ=255)
                                                                                                                      atatGCAGGTGCGCTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAACATATGTACTAGCGGTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTg            <  1:270746/143‑1 (MQ=255)
                                                                                                                       tatGCAGGTGCGCTTGGCATCCATAGTCCGGAGCACATTGCCGCCTGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAACGGTCATATGGCCGTGCACCTGTGCCACACCGGACGCATTTCTCCCGCCAGCCTGccggctgg           >  2:403204/1‑136 (MQ=255)
                                                                                                                        atGCAGGTGCACTTGACATCCATAGTCCGGATCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGc          <  2:414537/143‑1 (MQ=255)
                                                                                                                                gcgcTTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGc  >  1:1054233/1‑143 (MQ=255)
                                                                                                                                    |                                                                                                                                          
TGCGCAGTATTGAACCGGGTGACATTCCTACCCCGTTGATGGCGGAATACTATCGCCAACGTGCCAGTGCCGGTTTGATTATTAGTGAAGCCACGCAAATTTCTGCCCAGGCAAAAGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGC  >  NC_000913/1726744‑1727014

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