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,7550CT98.3% 185.8 / NA 58A80T (GCC→ACC) cysCadenosine 5'‑phosphosulfate kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/29);  minor base G (0/1);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGA  >  NC_000913/2873624‑2873884
                                                                                                                                   |                                                                                                                                 
tcCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTa                                                                                                                            <  2:593972/139‑1 (MQ=255)
     cgcgTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAAcc                                                                                                                       >  1:790714/1‑139 (MQ=255)
       cgTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCt                                                                                                                             >  2:926933/1‑131 (MQ=255)
       cgTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCt                                                                                                                             <  1:926933/131‑1 (MQ=255)
            gcgAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCg                                                                                                                <  1:1546264/139‑1 (MQ=255)
                    tCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATaa                                                                                                        >  2:335939/1‑139 (MQ=255)
                    tCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATaa                                                                                                        >  2:15710/1‑139 (MQ=255)
                                cgcgGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAAc                                                                                            >  2:309931/1‑139 (MQ=255)
                                       tgGCGAGATAAATGCGGTCAGTACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCg                                                                                     <  2:982204/139‑1 (MQ=255)
                                            agaTAAATGCGGTCAGCACCTCCTGTCCGGCTTCAATCATCAAATTCGCCACTTCACCGACGCAACTGATATTCTCTTTACGATCGGTATCGCGAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCagc                                                                                <  1:1664807/139‑3 (MQ=255)
                                               tAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAAt                                                                                                       >  2:268208/1‑113 (MQ=255)
                                               tAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAAt                                                                                                       <  1:268208/113‑1 (MQ=255)
                                               tAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcag                                                                             <  2:1174481/139‑1 (MQ=255)
                                               tAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcag                                                                             <  1:225902/139‑1 (MQ=255)
                                                 aaTGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcagca                                                                           <  1:335939/139‑1 (MQ=255)
                                                   tGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATGCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGGTCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGa                                                                         <  1:680961/139‑1 (MQ=255)
                                                       gTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACg                                                                     >  2:819868/1‑139 (MQ=255)
                                                                accaGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGc                                                            <  1:15710/139‑1 (MQ=255)
                                                                 ccaGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGcc                                                           >  1:1456178/1‑139 (MQ=255)
                                                                 ccaGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGcc                                                           >  1:1559306/1‑139 (MQ=255)
                                                                    gTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATc                                                                                <  2:834588/115‑1 (MQ=255)
                                                                    gTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATc                                                                                >  1:834588/1‑115 (MQ=255)
                                                                     tCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGt                                                       <  2:491739/139‑1 (MQ=255)
                                                                     tCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGt                                                       >  1:491739/1‑139 (MQ=255)
                                                                      ccGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAAt                                                                                                       >  1:1422916/1‑90 (MQ=255)
                                                                      ccGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAAt                                                                                                       <  2:1422916/90‑1 (MQ=255)
                                                                         gCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTAt                                                   >  1:962108/1‑139 (MQ=255)
                                                                               aCCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCga                                                                                                                     <  2:117411/67‑1 (MQ=255)
                                                                               aCCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCga                                                                                                                     >  1:117411/1‑67 (MQ=255)
                                                                               aCCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGt                                                                                                   <  1:27589/85‑1 (MQ=255)
                                                                               aCCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGt                                                                                                   >  2:27589/1‑85 (MQ=255)
                                                                               aCCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACg                                             >  2:1577370/1‑139 (MQ=255)
                                                                                         tCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGGATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAg                                   <  1:313424/139‑1 (MQ=255)
                                                                                          cGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCg                                                                                                                >  1:446461/1‑61 (MQ=255)
                                                                                          cGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCg                                                                                                                <  2:446461/61‑1 (MQ=255)
                                                                                          cGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCg                                                                                     >  2:478961/1‑88 (MQ=255)
                                                                                          cGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCg                                                                                     <  1:478961/88‑1 (MQ=255)
                                                                                          cGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAgc                                  <  1:819868/139‑1 (MQ=255)
                                                                                           gCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAgc                                  <  1:611719/138‑1 (MQ=255)
                                                                                           gCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAgc                                  >  2:611719/1‑138 (MQ=255)
                                                                                                         gcgacgGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGt                    >  2:298226/1‑138 (MQ=255)
                                                                                                         gcgacgGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGt                    <  1:298226/138‑1 (MQ=255)
                                                                                                          cgacgGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGc                                                            <  2:277226/97‑1 (MQ=255)
                                                                                                          cgacgGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGc                                                            >  1:277226/1‑97 (MQ=255)
                                                                                                            acgGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCa                                                                                  <  2:1601178/73‑1 (MQ=255)
                                                                                                            acgGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCa                                                                                  >  1:1601178/1‑73 (MQ=255)
                                                                                                            acgGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTg                  >  1:1441474/1‑137 (MQ=255)
                                                                                                            acgGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTg                  <  2:1441474/137‑1 (MQ=255)
                                                                                                                        tttACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCg    >  2:745899/1‑139 (MQ=255)
                                                                                                                          tACGATTGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCa                                                                                  <  2:1127463/59‑1 (MQ=255)
                                                                                                                          tACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCa                                                                                  >  1:1127463/1‑59 (MQ=255)
                                                                                                                          tACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAgc                                  >  2:750527/1‑107 (MQ=255)
                                                                                                                          tACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAgc                                  <  1:750527/107‑1 (MQ=255)
                                                                                                                          tACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGga  <  1:1577370/139‑1 (MQ=255)
                                                                                                                              aTCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCgg                           >  2:31758/1‑110 (MQ=255)
                                                                                                                              aTCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCgg                           <  1:31758/110‑1 (MQ=255)
                                                                                                                               tCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACg                                                             >  2:277130/1‑75 (MQ=255)
                                                                                                                               tCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACg                                                             <  1:277130/75‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                 
TCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGA  >  NC_000913/2873624‑2873884

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