Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,039,273 G→C 11.2% P42R (CCT→CGT)  PP_0897 ← class 1 fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,039,2730GC11.2% 93.2 / 5.6 36P42R (CCT→CGT) PP_0897class 1 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (18/14);  new base C (2/2);  total (20/16)
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.99e-01

GCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGAC  >  NC_002947/1039128‑1039419
                                                                                                                                                  |                                                                                                                                                  
gCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGGGTTGATCAGGTTCTGGGCGATGGAATCGCGTGCAG‑CAGGCGa                                                                                                                                                <  1:339862/150‑1 (MQ=255)
  tCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGAtt                                                                                                                                              >  2:346918/1‑150 (MQ=255)
    aTGGTGGCGCCATCCCAGCGCACGTCAATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGGAG‑CAGGCGATTcg                                                                                                                                            >  2:429650/1‑149 (MQ=255)
     tggtggCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCt                                                                                                                                           <  1:840823/150‑1 (MQ=255)
       gtggCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCCTGCTG‑CACGCGATTccatc                                                                                                                                         <  2:680934/150‑4 (MQ=255)
          gcgcCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGGGCGCCCGTCACGATAGCGGCTAGTTAAGCTCGTGCAT‑CAGGCGATTCCTCACGc                                                                                                                                      >  2:437857/1‑150 (MQ=255)
            gcCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGCCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAg                                                                                                                                    >  1:460583/1‑150 (MQ=255)
                 cccAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCCTGCTG‑CACGCGAtt                                                                                                                                              >  2:592309/1‑135 (MQ=255)
                 cccAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCCTGCTG‑CACGCGAtt                                                                                                                                              <  1:592309/135‑1 (MQ=255)
                   cAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGcc                                                                                                                             >  1:356424/1‑150 (MQ=255)
                   cAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATCACCGGTGGCGCAAATGCGAGAGTGGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCGCGCGGAGAGGcc                                                                                                                             >  1:278011/1‑150 (MQ=255)
                         aCGGCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTAGACGAGGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGccatgcac                                                                                                                       <  1:8913/150‑7 (MQ=255)
                             ccATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATgg                                                                                                                   >  2:172668/1‑150 (MQ=255)
                               aTGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGcc                                                                                                                 >  1:383749/1‑150 (MQ=255)
                                               acacGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACgg                                                                                                 >  2:434599/1‑150 (MQ=255)
                                                            aCCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑AAGGCGATTCCTCAGCCAGGTAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACga                                                                                    <  1:233691/150‑1 (MQ=255)
                                                                      tGGCAGATTGGACGATGGCCTGTGGCGCACATGCGCGAGTTGATCAGAATCGGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATGGc                                                                          >  1:411009/1‑150 (MQ=255)
                                                                      tGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGaaa                                                                                                        >  1:518133/1‑120 (MQ=255)
                                                                      tGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGaaa                                                                                                        <  2:518133/120‑1 (MQ=255)
                                                                      tGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAAtt                                                                            <  2:765545/148‑1 (MQ=255)
                                                                      tGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAAtt                                                                            >  1:765545/1‑148 (MQ=255)
                                                                        gCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACgg                                                                                                 <  1:631331/125‑1 (MQ=255)
                                                                        gCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACgg                                                                                                 >  2:631331/1‑125 (MQ=255)
                                                                                gACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGc                                                                                                                       >  2:301880/1‑95 (MQ=255)
                                                                                gACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGc                                                                                                                       <  1:301880/95‑1 (MQ=255)
                                                                                                   caTGCGCGAGTTGATCAGGATCTGGGCGATGGAATCG‑GCGAATTCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGgtc                                             <  1:437857/150‑1 (MQ=255)
                                                                                                        gcgAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTctt                                        >  2:103619/1‑150 (MQ=255)
                                                                                                           aGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTcttgct                                     <  2:356424/150‑1 (MQ=255)
                                                                                                            gTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTcttgctt                                    >  2:693227/1‑150 (MQ=255)
                                                                                                              tGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGa                                  <  1:346918/150‑1 (MQ=255)
                                                                                                                aTCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCa                                                 >  2:179404/1‑133 (MQ=255)
                                                                                                                aTCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCa                                                 <  1:179404/133‑1 (MQ=255)
                                                                                                                         tGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAg                                                <  2:240988/125‑1 (MQ=255)
                                                                                                                         tGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAg                                                >  1:240988/1‑125 (MQ=255)
                                                                                                                                    aTCGCCTGCTG‑CACGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGctcc            >  1:75752/1‑150 (MQ=255)
                                                                                                                                       gcgTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCctct         <  2:383749/150‑1 (MQ=255)
                                                                                                                                            cag‑cagGCGATTCCTCACGCAGATAGGACTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACg                             >  1:340508/1‑125 (MQ=255)
                                                                                                                                            cag‑cagGCGATTCCTCACGCAGATAGGACTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACg                             <  2:340508/125‑1 (MQ=255)
                                                                                                                                              g‑cagGCGATTCCTCACGCAGATAGGAGTGGTTCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGAc  <  1:429650/150‑1 (MQ=255)
                                                                                                                                                 agGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCctct         >  2:414556/1‑141 (MQ=255)
                                                                                                                                                 agGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCctct         <  1:414556/141‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCTCATGGTGGCGCCATCCCAGCGCACGTCCATGCCCACGCGGATGAACACGGTGACGATACCGGTGTCCTGGCAGATTGGACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAG‑CAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGAC  >  NC_002947/1039128‑1039419

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