Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,873,755 C→T A80T (GCC→ACC)  cysC ← adenosine 5'‑phosphosulfate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,873,7550CT97.3% 112.8 / NA 37A80T (GCC→ACC) cysCadenosine 5'‑phosphosulfate kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (13/23);  minor base A (1/0);  total (14/23)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCC  >  NC_000913/2873624‑2873890
                                                                                                                                   |                                                                                                                                       
tcCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTa                                                                                                                                  <  2:670427/139‑1 (MQ=255)
       cgTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCga                                                                                                                           <  1:267303/139‑1 (MQ=255)
       cgTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCga                                                                                                                           <  2:642703/139‑1 (MQ=255)
          tCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTTAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGAt                                                                                                                        <  2:533470/139‑1 (MQ=255)
          tCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGAt                                                                                                                        >  2:251505/1‑139 (MQ=255)
            gcgAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCg                                                                                                                      <  1:34248/139‑1 (MQ=255)
                aCCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGc                                                                                                                  <  1:234855/139‑1 (MQ=255)
                    tCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATaa                                                                                                              >  2:34454/1‑139 (MQ=255)
                     cTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTAACCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCAACGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAAt                                                                                                             <  2:550708/139‑1 (MQ=255)
                                              aTAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATc                                                                                      >  2:558886/1‑137 (MQ=255)
                                              aTAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATc                                                                                      <  1:558886/137‑1 (MQ=255)
                                               tAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAATACCGAGATCGCTGCAGAATCCGTGGCGAACATTGTCGCCATCcag                                                                                   <  2:536814/139‑1 (MQ=255)
                                                       gTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTTTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATcc                                                                                                           <  2:659107/107‑1 (MQ=255)
                                                       gTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATcc                                                                                                           >  1:659107/1‑107 (MQ=255)
                                                        tCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGt                                                                          >  1:914255/1‑139 (MQ=255)
                                                          aGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTAc                                                                        <  1:585198/139‑1 (MQ=255)
                                                            caccacCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCg                                                                                                                      <  2:26270/91‑1 (MQ=255)
                                                            caccacCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCg                                                                                                                      >  1:26270/1‑91 (MQ=255)
                                                                accaGACCGGCTTCAACCATCAAATTCGCGACTTCACCGACGCGTCGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGc                                                                  <  2:898040/139‑1 (MQ=255)
                                                                      ccGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGtt                                                            >  1:338098/1‑139 (MQ=255)
                                                                      ccGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACAGACGCCGAGtt                                                            >  1:602984/1‑139 (MQ=255)
                                                                            tCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTa                                                      <  2:658292/139‑1 (MQ=255)
                                                                                    cAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcagcag                                                                                >  1:599247/1‑105 (MQ=255)
                                                                                    cAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcagcag                                                                                <  2:599247/105‑1 (MQ=255)
                                                                                         tCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAg                                         <  1:34454/139‑1 (MQ=255)
                                                                                                          cgacgGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCg                                                                                                     >  2:385857/1‑62 (MQ=255)
                                                                                                          cgacgGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCg                                                                                                     <  1:385857/62‑1 (MQ=255)
                                                                                                          cgacgGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGAc                             <  2:907329/134‑1 (MQ=255)
                                                                                                          cgacgGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGAc                             >  1:907329/1‑134 (MQ=255)
                                                                                                                atatTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTAc                  <  2:602984/139‑1 (MQ=255)
                                                                                                                atatTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTAc                  <  1:900826/139‑1 (MQ=255)
                                                                                                                         ttACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcag                                                                                   <  1:31282/65‑1 (MQ=255)
                                                                                                                         ttACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcag                                                                                   >  2:31282/1‑65 (MQ=255)
                                                                                                                             gATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGagag     <  2:338098/139‑1 (MQ=255)
                                                                                                                               tCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGc   >  1:699607/1‑139 (MQ=255)
                                                                                                                               tCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGc   <  2:669130/139‑1 (MQ=255)
                                                                                                                                cGGAATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTACTTCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGcc  >  2:599355/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                       
TCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCC  >  NC_000913/2873624‑2873890

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