Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,803,836 T→C L112P (CTA→CCA)  nrdF → ribonucleoside‑diphosphate reductase 2, beta subunit, ferritin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,803,8360TC97.4% 136.7 / ‑3.9 39L112P (CTA→CCA) nrdFribonucleoside‑diphosphate reductase 2, beta subunit, ferritin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (17/21);  total (17/22)
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.60e-01

ACGCTGCTCGACACGCTGCAAAATGTTATCGGCGCGCCTTCTCTGATGCCCGATGCACTCACGCCTCATGAAGAAGCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCTATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTT  >  NC_000913/2803691‑2803980
                                                                                                                                                 |                                                                                                                                                
aCGCTGCTCGACACGCTGCAAAATGTTATCGGCGCGCCTTCTCTGATGCCCGATGCACTCACGCCTCATGAAGAAGCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGcc                                                                                                                                             >  1:173979/1‑151 (MQ=255)
     gctCGACACGCTGCAAAATGTTATCGGCGCGCCTTCTCTGATGCCCGATGCACTCACGCCTCATGAAGAAGCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGAcc                                                                                                                                        >  2:516483/1‑151 (MQ=255)
        cGACACGCTGCAAAATGTTATCGGCGCGCCTTCTCTGATGCCCGATGCACTCACGCCTCATGAAGAAGCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCaa                                                                                                                                      <  2:472842/150‑1 (MQ=255)
          acacGCTGCAAAATGTTATCGGCGCGCCTTCTCTGATGCCCGATGCACTCACGCCTCATGAAGAAGCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGa                                                                                                                                   <  2:181207/151‑1 (MQ=255)
             cGCTGCAAAATGTTATCGGCGCGCCTTCTCTGATGCCCGATGCACTCACGCCTCATGAAGAAGCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGt                                                                                                                                >  1:302391/1‑151 (MQ=255)
                  cAAAATGTTATCGGCGCGCCTTCTCTGATGCCCGATGCACTCACGCCTCATGAAGAAGCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATg                                                                                                                           <  2:81337/151‑1 (MQ=255)
                                cgcgcCTTCTCTGATGCCCGATGCACTCACGCCTCATGAAGAAGCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTgg                                                                                                             >  1:77311/1‑151 (MQ=255)
                                      ttCTCTGATGCCCGATGCACTCACGCCTCATGAAGAAGCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTga                                                                                                        >  2:24696/1‑150 (MQ=255)
                                                  cGATGCACTCACGCCTCATGAAGAAGCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCg                                                                                           <  1:150421/151‑1 (MQ=255)
                                                   gATGCACTCACGCCTCATGAAGAAGCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGt                                                                                          >  2:200777/1‑151 (MQ=255)
                                                     tGCACTCTCGCCTCATGAAGAAGCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCTAAGATGTCGATGCCGCCTACGCCTGGAGTGTAGAAAACGCAc                                                                                             <  2:515081/146‑1 (MQ=255)
                                                           cACGCCTCATGAAGAAGCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCg                                                                                   <  2:612726/150‑1 (MQ=255)
                                                                  cATGAAGAAGCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTc                                                                           <  2:96870/151‑1 (MQ=255)
                                                                    tGAAGAAGCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTc                                                                           >  1:399915/1‑149 (MQ=255)
                                                                        gaagCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGatta                                                                     <  1:267592/151‑1 (MQ=255)
                                                                         aagCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGattat                                                                    >  1:515410/1‑151 (MQ=255)
                                                                             ggTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTc                                                                  >  2:298099/1‑149 (MQ=255)
                                                                               tattatCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGc                                                               >  1:583293/1‑150 (MQ=255)
                                                                                attatCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCa                                                              <  1:79888/150‑1 (MQ=255)
                                                                                attatCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCa                                                              >  2:240237/1‑150 (MQ=255)
                                                                                   atCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACAt                                                          >  2:281264/1‑151 (MQ=255)
                                                                                   atCGAATATCAGCTTTACGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCTATGCCAGACCAAAGATGTCGATGCCGACTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACAt                                                          <  2:431031/151‑1 (MQ=255)
                                                                                     cGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACAtt                                                         <  1:22634/150‑1 (MQ=255)
                                                                                        atatCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATcg                                                     <  1:244116/151‑1 (MQ=255)
                                                                                           tCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATCGCgg                                                  >  1:211120/1‑151 (MQ=255)
                                                                                            cAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATCGCGGt                                                 >  1:331342/1‑151 (MQ=255)
                                                                                                 ttATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATCGCGGtgatg                                             <  1:73727/150‑1 (MQ=255)
                                                                                                 ttATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATCGCGGtgatg                                             <  1:266017/150‑1 (MQ=255)
                                                                                                  tATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATCGCGGtgatga                                            <  2:576757/150‑1 (MQ=255)
                                                                                                  tATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATCGCGGtgatga                                            <  2:241199/150‑1 (MQ=255)
                                                                                                   aTGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATCGCGGTGATGATc                                          <  1:139216/151‑1 (MQ=255)
                                                                                                         gCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATCGCGGTGATGATCCGCTg                                     >  1:553635/1‑150 (MQ=255)
                                                                                                             ttCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATCGCGGTGATGATCCGCTGaaaa                                 >  1:381181/1‑150 (MQ=255)
                                                                                                                  gCCCGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATCGCGGTGATGATCCGCTGAAAAAGaaaa                           <  1:525023/151‑1 (MQ=255)
                                                                                                                   cccGCTCTTACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATCGCGGTGATGATCCGCTGAAAAAGAAAAt                          <  1:616450/151‑1 (MQ=255)
                                                                                                                           tACAGTTCGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCa                     <  1:273393/148‑1 (MQ=255)
                                                                                                                                  cGATTTTCTCGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCt           >  1:3191/1‑151 (MQ=255)
                                                                                                                                           cGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCttt  <  1:503669/151‑1 (MQ=255)
                                                                                                                                           cGACGCCATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCttt  <  1:103506/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACGCTGCTCGACACGCTGCAAAATGTTATCGGCGCGCCTTCTCTGATGCCCGATGCACTCACGCCTCATGAAGAAGCGGTATTATCGAATATCAGCTTTATGGAAGCGGTTCATGCCCGCTCTTACAGTTCGATTTTCTCGACGCTATGCCAGACCAAAGATGTCGATGCCGCCTACGCCTGGAGTGAAGAAAACGCACCGTTGCAGCGAAAAGCTCAGATTATTCAGCAACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTT  >  NC_000913/2803691‑2803980

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