Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,254,420 C→A 9.9% G440C (GGC→TGC)  PP_2852 ← sulfatase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,254,4200CA 9.9% 143.1 / 4.9 48G440C (GGC→TGC) PP_2852sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (20/23);  new base A (2/3);  total (22/26)
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.72e-01

GGCCGTAAGTGCCGCGCTTGATGGTTGGCAGTTGGGCCTGCTCCGGCGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGT  >  NC_002947/3254273‑3254558
                                                                                                                                                   |                                                                                                                                          
ggCCGTAAGTGCCGCGCTTGATGGTTGGCAGTTGGGCCTGCTCCGGCGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGccc                                                                                                                                           >  6:65735/1‑149 (MQ=255)
ggCCGTAAGTGCCGCGCTTGATGGTTGGCAGTTGGGCCTGCTCCGGCGCCAATGTAAAGTTGAAGCCCCAGGACGCGGCCAGGCGCACGCCATCGATGCCGTGAGACAAGTCGGAAGTCACCACCACCAGTTTGTCCTTGAGGAAGccc                                                                                                                                           >  1:314263/1‑149 (MQ=255)
ggCCGTAAGTGCCGCGCTTGATGGTTGGCAGTTGGGCCTGCTCCGGCGCCAATGTAAAGTTGAAGCCCCAGGACGAGGCCAGGCGCACGGCATCGATGCCGTGAGACTAGTAGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGccc                                                                                                                                           >  1:319419/1‑149 (MQ=255)
   cGTAAGTGCCGCGCTTGATGGTTGGCAGTTGGGCCTGCTCCGGCGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTgg                                                                                                                                        <  2:46326/149‑1 (MQ=255)
     tAAGGGCCGCGCTTGATGGTTGGCAGTTGGGCCTGCTCCGGCGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGGGCGTCCTCAAGGACACACTggtg                                                                                                                                      <  8:213017/149‑3 (MQ=255)
         tGCCGCGCTTGATGGTTGGCAGTTGGGCCTGCTCCGGCGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTc                                                                                                                                   >  6:304958/1‑148 (MQ=255)
         tGCCGCGCTTGATGGTTGGCAGTTGGGCCTGCTCCGGCGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGGGCGTCCTCAAGGACACACTggtggtgg                                                                                                                                  <  6:414624/149‑7 (MQ=255)
               gcTTGATGGTTGGCAGTTGGGCCTGCTCCGGCGCCAATGTAAGGTTGAAGCCCCAGGACAAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTg                                                                                                                            <  6:291341/149‑1 (MQ=255)
                   gATGGTTGGCAGTTGGGCCTGCTCCGGCGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCa                                                                                                                         <  7:224288/148‑1 (MQ=255)
                       gTTGGCAGTTGGGCCTGCTCCGGCGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGaaa                                                                                                                    <  5:65735/149‑1 (MQ=255)
                            cAGTTGGGCCTGCTCCGGCGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGccc                                                                                                               >  7:189760/1‑149 (MQ=255)
                            cAGTTGGGCCTGCTCCGGCGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGGGCGTCCTCAAGGACACACTggtgg                                                                                                                                     >  6:291622/1‑124 (MQ=255)
                            cAGTTGGGCCTGCTCCGGCGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGGGCGTCCTCAAGGACACACTggtgg                                                                                                                                     <  5:291622/127‑4 (MQ=255)
                              gTTGGGCCTGCTCCGGCGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGcccc                                                                                                              >  5:261817/1‑148 (MQ=255)
                                 gggCCTGCTCCGGCGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACt                                                                                                          >  3:20120/1‑149 (MQ=255)
                                  ggCCTGCTCCGGCGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACt                                                                                                          <  6:372377/148‑1 (MQ=255)
                                        cTCCGGCGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTc                                                                                                     <  8:189760/147‑1 (MQ=255)
                                          ccGGCGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATcc                                                                                                 <  5:143108/149‑1 (MQ=255)
                                             gcgcCAATGTAAGGTTGAAGCCCCTGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCGTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAgg                                                                                              <  5:304958/149‑1 (MQ=255)
                                                ccAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAg                                                                                           <  2:146951/149‑1 (MQ=255)
                                                    tGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGa                                                                                       >  8:306464/1‑149 (MQ=255)
                                                    tGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGa                                                                                       >  2:129905/1‑149 (MQ=255)
                                                        tggTTGAAGCCCCAGTACGACGCCAGGCGCAAGCCATCTATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGc                                                                                   <  5:183782/148‑1 (MQ=255)
                                                         ggTGGAAGCCCCAGGTCGAGGCCAGGCGCACGCCATCGATGGCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCt                                                                                  <  4:264326/149‑1 (MQ=255)
                                                                 ccccAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGttt                                                                          <  6:261817/149‑1 (MQ=255)
                                                                     aGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTATGCGACGGCTGCCTGTTTCGgt                                                                      <  2:319419/149‑1 (MQ=255)
                                                                     aGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGgt                                                                      <  2:314263/149‑1 (MQ=255)
                                                                            ggCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCg                                                                >  2:306255/1‑148 (MQ=255)
                                                                             gCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGaa                                                              <  1:129905/149‑1 (MQ=255)
                                                                              ccAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCTCCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGCTTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAg                                                             <  5:72876/149‑1 (MQ=255)
                                                                                aGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAgcg                                                           >  2:151510/1‑149 (MQ=255)
                                                                                aGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCACTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAgcg                                                           >  3:24476/1‑149 (MQ=255)
                                                                                  gcgcACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGt                                                                                                      >  8:287955/1‑104 (MQ=255)
                                                                                  gcgcACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGt                                                                                                      <  7:287955/104‑1 (MQ=255)
                                                                                    gcACGCCATCGATGCCGTGAGACTCGTCGCAAGTAACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCACGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTc                                                        >  1:389667/1‑148 (MQ=255)
                                                                                      aCGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAg                                                     <  5:206839/149‑1 (MQ=255)
                                                                                           aTCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACt                                                >  2:233302/1‑149 (MQ=255)
                                                                                                     tGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGgcgc                                      >  7:386327/1‑149 (MQ=255)
                                                                                                           tcgtcgGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCATGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGaa                                  >  4:243519/1‑147 (MQ=255)
                                                                                                           tcgtcgGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTcg                                              >  7:196332/1‑135 (MQ=255)
                                                                                                           tcgtcgGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTcg                                              <  8:196332/135‑1 (MQ=255)
                                                                                                             gtcgGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAgg                              >  4:34761/1‑149 (MQ=255)
                                                                                                             gtcgGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAgg                              <  1:306255/149‑1 (MQ=255)
                                                                                                                       accaccaccAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAg                                                     >  7:138945/1‑116 (MQ=255)
                                                                                                                       accaccaccAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAg                                                     <  8:138945/116‑1 (MQ=255)
                                                                                                                            caccaGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCt                           >  2:310234/1‑137 (MQ=255)
                                                                                                                            caccaGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTACCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCt                           <  1:310234/137‑1 (MQ=255)
                                                                                                                              ccaGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGACCACg             >  3:231037/1‑149 (MQ=255)
                                                                                                                               caGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGCTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGACCGGCGCCGAATAGGGCTGGTGGGTGCCCACgg            >  5:15065/1‑149 (MQ=255)
                                                                                                                                  tgtgtCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGa         <  7:306464/149‑1 (MQ=255)
                                                                                                                                         ttGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGt  >  3:31206/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
GGCCGTAAGTGCCGCGCTTGATGGTTGGCAGTTGGGCCTGCTCCGGCGCCAATGTAAGGTTGAAGCCCCAGGACGAGGCCAGGCGCACGCCATCGATGCCGTGAGACTCGTCGGAAGTCACCACCACCAGTGTGTCCTTGAGGACGCCCTGGCGCTCGAGATTGTCAAGAAATGCCCCGACTGCGTCATCCAGGTAGGCGACGGCTGCCTGTTTCGGTGTGTCGAAGCGCTCGAGGTACTCGGCCGGCGCCGAATAGGGCTGGTGGGTGCCCACGGTGAGCAGCGT  >  NC_002947/3254273‑3254558

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