Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,491,142:1 +A 7.8% coding (805/1902 nt) cysNC → sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,491,1421.A 7.8% 194.2 / 11.5 51coding (805/1902 nt)cysNCsulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase
Reads supporting (aligned to +/- strand):  ref base . (21/26);  new base A (2/2);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATC  >  NC_002947/1490997‑1491285
                                                                                                                                                   |                                                                                                                                               
gTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCg‑tg                                                                                                                                               <  2:128544/148‑1 (MQ=255)
  aCGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCg‑tggtg                                                                                                                                            >  3:205370/1‑149 (MQ=255)
   cGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGc                                                                                                                                           <  1:144449/149‑1 (MQ=255)
   cGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCA‑CGCCACTGGCGATGGTGc                                                                                                                                           <  2:113410/149‑1 (MQ=255)
    gCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGca                                                                                                                                          >  5:100881/1‑149 (MQ=255)
            gACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGt                                                                                                                                   >  7:104948/1‑148 (MQ=255)
               gCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGAc                                                                                                                                >  7:83498/1‑148 (MQ=255)
                    tGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGAtt                                                                                                                          <  2:72338/149‑1 (MQ=255)
                          tCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGtg                                                                                                                    <  8:104948/149‑1 (MQ=255)
                             ccGAAACC‑TCGACGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGtgctg                                                                                                                 >  4:242155/2‑149 (MQ=255)
                             ccGAAACC‑TCGACGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGATGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGtgctg                                                                                                                 >  5:225308/2‑149 (MQ=255)
                              cGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCg‑tg                                                                                                                                               >  5:159844/1‑118 (MQ=255)
                              cGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCg‑tg                                                                                                                                               <  6:159844/118‑1 (MQ=255)
                                   ccGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCg                                                                                                           <  5:119749/149‑1 (MQ=255)
                                      tCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                          <  5:155611/147‑1 (MQ=255)
                                      tCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                          >  6:155611/1‑147 (MQ=255)
                                         acggTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCaa                                                                                                      <  3:242155/147‑1 (MQ=255)
                                           ggTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAaga                                                                                                    <  3:184281/149‑1 (MQ=255)
                                                  gCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCg                                                                                             <  4:205370/149‑1 (MQ=255)
                                                   cTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCgg                                                                                            <  6:269864/149‑1 (MQ=255)
                                                    tGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCggg                                                                                           >  7:35026/1‑149 (MQ=255)
                                                        cGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCaa                                                                                       <  2:121508/149‑1 (MQ=255)
                                                          cAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGt                                                                                     <  6:100881/149‑1 (MQ=255)
                                                           aaCTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTc                                                                                    >  2:62976/1‑149 (MQ=255)
                                                               tCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATc                                                                                >  5:260711/1‑149 (MQ=255)
                                                               tCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCAGTCGGGCAAGAGCAGCCGGGTCAAGTCCATc                                                                                >  2:193205/1‑149 (MQ=255)
                                                                  ccGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTc                                                                             <  1:107144/149‑1 (MQ=255)
                                                                  ccGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTc                                                                             <  4:142017/149‑1 (MQ=255)
                                                                   cGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCa                                                                            <  3:68104/149‑1 (MQ=255)
                                                                   cGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCa                                                                            <  8:35026/149‑1 (MQ=255)
                                                                    gACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCAc                                                                           >  1:59584/1‑149 (MQ=255)
                                                                        tGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                          >  8:164556/1‑114 (MQ=255)
                                                                                 cGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCTGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGaa                                                              <  1:62976/149‑1 (MQ=255)
                                                                                       aGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGa                                                         >  7:37835/1‑148 (MQ=255)
                                                                                        gTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGTACTGGaa                                                        >  6:137538/1‑148 (MQ=255)
                                                                                               aaCCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTgg                                                >  4:86337/1‑149 (MQ=255)
                                                                                               aaCCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTgg                                                <  1:193205/149‑1 (MQ=255)
                                                                                               aaCCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTg                                                 <  8:83498/148‑1 (MQ=255)
                                                                                                    tCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccag                                           <  2:59584/149‑1 (MQ=255)
                                                                                                       cAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccaggc                                         >  2:7492/1‑148 (MQ=255)
                                                                                                           cTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTTCGTAGGTGacga                                                            <  3:128695/125‑4 (MQ=255)
                                                                                                           cTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTTCGTAGGTGacga                                                            >  4:128695/1‑122 (MQ=255)
                                                                                                              aaCTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGt                                 <  1:108088/149‑1 (MQ=255)
                                                                                                              aaCTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTCGAAAACGCTGGCCCAGGCCAGGCCGt                                 <  1:7492/149‑1 (MQ=255)
                                                                                                                cTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGGCAAGTCCAGCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCCGGCCCGGCCGTgc                               >  4:190179/1‑148 (MQ=255)
                                                                                                                       ggCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGAc                        <  5:137538/149‑1 (MQ=255)
                                                                                                                               gggCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATgg                   >  2:49340/1‑146 (MQ=255)
                                                                                                                                 gCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGAc              <  2:162108/149‑1 (MQ=255)
                                                                                                                                    accacGCCACTGGCGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccaggccagg                                     >  7:156606/5‑124 (MQ=255)
                                                                                                                                    accacGCCACTGGCGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccaggccagg                                     <  8:156606/120‑1 (MQ=255)
                                                                                                                                    accacGCCACTGGCGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGAc              >  8:228761/5‑147 (MQ=255)
                                                                                                                                     cATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGAc              >  5:22418/1‑145 (MQ=255)
                                                                                                                                       tCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATc        <  6:225308/149‑1 (MQ=255)
                                                                                                                                          ccAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGAc     <  6:260711/149‑1 (MQ=255)
                                                                                                                                             gTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACAtc  >  7:142253/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATC  >  NC_002947/1490997‑1491285

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