Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,429,296 G→A 10.1% A112T (GCT→ACT) ‡ PP_3923 → phosphoglycerate mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,429,2960GA10.1% 129.0 / 5.8 50A112T (GCT→ACT) ‡PP_3923phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base G (22/22);  new base A (3/2);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

CGCCAACGGCTGCCCGGTGCCGGCGGTTGAGACCGACGCCGCATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTGGCATTCGTTTACGTCACGCGTACAAGGGGGCTTGCAACGCGTGCTGGACA  >  NC_002947/4429163‑4429440
                                                                                                                                       |                                                                                                                                                
cGCCAACGGCTGCCCGGTGCCGGCGGTTGAGACCGACGCCGCATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCg                                                                                                                                     >  4:358534/1‑147 (MQ=255)
 gCCAACGGCTGCCCGGTGCCGGCGGTTGAGACCGACGCCGCATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAg                                                                                                                                  >  2:317913/1‑149 (MQ=255)
   cAACGGCTGCCCGGTGCCGGCGGTTGAGACCGACGCCGCATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGc                                                                                                                                 >  5:233587/1‑148 (MQ=255)
    aaCGGCTGCCCGGTGCCGGCGGTTGAGACCGACGCCGCATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGCCCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGc                                                                                                                                 >  3:83317/1‑147 (MQ=255)
    aaCGGCTGCCCGGGGCCGGCGGTTGAGACCGACGCCGCATTCCATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGa                                                                                                                               <  8:58357/149‑1 (MQ=255)
          tGCCCGGTGCCGGCGGTTGAGACCGACGCCGCATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGCAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCa                                                                                                                         <  5:317815/149‑1 (MQ=255)
            cccGGTGCCGGCGGTTGAGACCGACGCCGCATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAgc                                                                                                                       >  7:312111/1‑149 (MQ=255)
             ccGGTGCCGGCGGTTGAGACCGACGCCGCATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAgcg                                                                                                                      >  8:177055/1‑149 (MQ=255)
              cGGTGCCGGCGGTTGAGACCGACGCCGCATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAgcgc                                                                                                                     <  2:191294/149‑1 (MQ=255)
               ggTGCCGGCGGTTGAGACCGACGCCGCATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGcc                                                                                                                    <  7:2387/149‑1 (MQ=255)
                gTGCCGGCGGTTGAGACCGACGCCGCATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCt                                                                                                                   <  2:365540/149‑1 (MQ=255)
                    cggcggTTGAGACCGACGCCGCATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTc                                                                                                               <  4:158218/149‑1 (MQ=255)
                     ggcggTTGAGACCGACGCCGCATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTcg                                                                                                              >  3:5832/1‑149 (MQ=255)
                                 cGACGCCGCATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGc                                                                                                  <  8:14077/149‑1 (MQ=255)
                                   aCGCCGCATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGCTGGTGCAg                                                                                                >  7:286265/1‑149 (MQ=255)
                                        gcATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTg                                                                                           <  8:312111/149‑1 (MQ=255)
                                         cATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTgg                                                                                          >  1:214824/1‑149 (MQ=255)
                                                 gAGTTCGATGCCGACGGGGTGATCCGCGCCCTCTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATgg                                                                                  >  3:96294/1‑149 (MQ=255)
                                                   gTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATACCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGc                                                                                 >  2:42038/1‑148 (MQ=255)
                                                    ttCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGc                                                                                 >  1:160288/1‑147 (MQ=255)
                                                    ttCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGa                                                                               <  6:233587/149‑1 (MQ=255)
                                                      cGATGCCGACCCGGAGATCCGCGCCCTGTTGCCCGACCTGCAGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGAATCAGCGCCTGTACGCGTTGATGGTGCAGCGCTGGCACGATGGCGAg                                                                              >  4:285176/1‑148 (MQ=255)
                                                          gCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTACTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGc                                                                         >  6:41818/1‑149 (MQ=255)
                                                           ccGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGcc                                                                        <  1:317913/149‑1 (MQ=255)
                                                            cGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGcc                                                                        >  1:368073/1‑148 (MQ=255)
                                                                   gTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCgca                                                                                                                                           >  5:222221/1‑74 (MQ=255)
                                                                   gTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCgca                                                                                                                                           <  6:222221/74‑1 (MQ=255)
                                                                     gATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGcc                                                              >  6:234823/1‑149 (MQ=255)
                                                                                tGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTgg                                                   >  4:147361/1‑149 (MQ=255)
                                                                                     cccGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTGGCATTc                                              >  7:323016/1‑149 (MQ=255)
                                                                                        gACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTGGCATTCGtt                                           <  7:177055/149‑1 (MQ=255)
                                                                                          ccTGCTGCCGGAAGACCGCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTGGCATTCGtt                                           <  1:319877/149‑1 (MQ=255)
                                                                                           ctgctgCCTGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGc                                                               <  8:281537/126‑1 (MQ=255)
                                                                                           ctgctgCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGc                                                               >  7:281537/1‑126 (MQ=255)
                                                                                              ctgcCGGAAGA‑‑GCCCGATGGCCTGCATGTGCTGCGCAATTCTGCGTAGCATCGGAGCGAGTTTCAGCGCCTGTTCTCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTGGCATTCGTTTACg                                       <  3:285176/147‑1 (MQ=255)
                                                                                                  cGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTGGCATTCGTTTACGTCAcgc                                 <  8:323016/149‑1 (MQ=255)
                                                                                                        a‑‑gCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTGGCATTCGTTTACGTCACGCGTAc                             >  8:220967/1‑147 (MQ=255)
                                                                                                            cccGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTGGCATTCGTTTACGTCACGCGTACAAgg                         <  4:96294/149‑1 (MQ=255)
                                                                                                                 tGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTGGCATTCGTTTACGTCACGCGTACAAGGGGGCt                    <  5:260298/149‑1 (MQ=255)
                                                                                                                       tgcgatgctgcGCAGCATTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTGGCATTCGTTTACGTCACGCGTACAAGGGGGCTTGCAAc              >  2:353545/6‑149 (MQ=255)
                                                                                                                        cATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTGGCATTCGTTTACGTCACGCGTACAAGGGGGCTTGCaa               <  2:214824/147‑1 (MQ=255)
                                                                                                                        cATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTGGCATTCGTTTACGTCACGCGTACAAGGGGGCTTGCAAc              <  2:368073/148‑1 (MQ=255)
                                                                                                                          tgtgCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTGGCATTCGTTTACGTCACGCGTACAAGGGGGCTTGCAAcgcg           >  5:225773/1‑149 (MQ=255)
                                                                                                                          gatgctgcGCAGCATTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGG‑CTGGAAACCTGGCATTCGTTTACGTCACGCGTACAAggggg                      <  7:288040/135‑1 (MQ=255)
                                                                                                                          gatgctgcGCAGCATTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGG‑CTGGAAACCTGGCATTCGTTTACGTCACGCGTACAAggggg                      >  8:288040/3‑137 (MQ=255)
                                                                                                                           atgctgcGCAGCATTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTg                                                    >  6:201222/2‑107 (MQ=255)
                                                                                                                           atgctgcGCAGCATTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTg                                                    <  5:201222/106‑1 (MQ=255)
                                                                                                                            tgctgcGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTGGCATTCGTTTACGTCACGCGTACAAGGGGGCTTGCAACGCGTg         <  2:160288/149‑1 (MQ=255)
                                                                                                                              ctgcGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTGGCATTCGTTTACGTCACGCGTACAAGGGGGCTTGCAACGCGTGc        <  6:105194/148‑1 (MQ=255)
                                                                                                                                   cAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTGGCATTCGTTTACGTCACGCGTACAAGGGGGCTTGCAACGCGTGCTGGACa  <  4:5832/149‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                
CGCCAACGGCTGCCCGGTGCCGGCGGTTGAGACCGACGCCGCATTCAATGAGTTCGATGCCGACGGGGTGATCCGCGCCCTGTTGCCCGACCTGCTGCCGGAAGA‑‑GCCCGATGCCCTGCATGTGCTGCGCAATGCTGCGCAGCATCGCAGCGAGTTTCAGCGCCTGTTCGCGTTGATGGTGCAGCGCTGGCACGATGGCGAGCATGCCGATAATGGCCTGGAAACCTGGCATTCGTTTACGTCACGCGTACAAGGGGGCTTGCAACGCGTGCTGGACA  >  NC_002947/4429163‑4429440

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: ‑