Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,631,134 A→G 10.5% N205D (AAC→GAC)  putA → proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,631,1340AG10.5% 98.9 / 3.4 38N205D (AAC→GAC) putAproline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (16/18);  new base G (2/2);  total (18/20)
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.10e-01

AGGAAGGCGTGGCGCTGATGTGCCTGGCCGAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGG  >  NC_002947/5630994‑5631281
                                                                                                                                            |                                                                                                                                                   
aGGAAGGCGTGGCGCTGATGTGCCTGGCCGAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCGACGCCGCCa                                                                                                                                            >  1:416263/1‑150 (MQ=255)
     ggCGTGGCGCTGATGTGCCTGGCCGAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCGACGCCGCCACCTgg                                                                                                                                       <  2:416263/150‑1 (MQ=255)
     ggCGTGGCGCTGATGTGCCTGGCCGAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTgg                                                                                                                                       <  2:481829/150‑1 (MQ=255)
            cgcTGATGTGCCTGGCCGAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGgctgc                                                                                                                                <  2:362040/150‑1 (MQ=255)
                     gCCTGGCCGAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCg                                                                                                                       >  2:534531/1‑150 (MQ=255)
                      ccTGGCCGAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCgg                                                                                                                      >  1:574803/1‑150 (MQ=255)
                            cGAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCt                                                                                                                >  2:3740/1‑150 (MQ=255)
                                gCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTgg                                                                                                                                       >  1:354390/1‑123 (MQ=255)
                                gCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGCGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTgg                                                                                                                                       <  2:354390/123‑1 (MQ=255)
                                  ccTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAg                                                                                                          <  1:544992/150‑1 (MQ=255)
                                       tgcGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCAccc                                                                                                     <  1:576441/150‑1 (MQ=255)
                                          gTATCCCCTACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATa                                                                                                  >  2:468689/1‑150 (MQ=255)
                                           tATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATaa                                                                                                 <  1:534531/150‑1 (MQ=255)
                                            aTCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAAc                                                                                                <  2:59341/150‑1 (MQ=255)
                                               cccGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGaa                                                                                             >  1:489765/1‑150 (MQ=255)
                                                 cGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAAc                                                                                           >  1:89610/1‑150 (MQ=255)
                                                      aaGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGcc                                                                                      >  1:451710/1‑150 (MQ=255)
                                                      aaGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGcc                                                                                      <  1:209903/150‑1 (MQ=255)
                                                                tGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACctcctc                                                                            <  2:451710/150‑1 (MQ=255)
                                                                                  cAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCGACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCgg                                                          >  1:306652/1‑150 (MQ=255)
                                                                                     gATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAAc                                                                                                >  2:175847/1‑109 (MQ=255)
                                                                                     gATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAAc                                                                                                <  1:175847/109‑1 (MQ=255)
                                                                                           tACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAgcgg                                                 >  2:419948/1‑150 (MQ=255)
                                                                                           tACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAgcgg                                                 >  1:24481/1‑150 (MQ=255)
                                                                                                   aCTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCga                                         <  2:24481/150‑1 (MQ=255)
                                                                                                         aGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCATGCTGGTCAGCACCCATAACGAAACCGGCATCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATcc                                   <  2:89610/150‑1 (MQ=255)
                                                                                                          gCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCg                                  >  2:544361/1‑150 (MQ=255)
                                                                                                             ccACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCaa                               <  2:451184/149‑1 (MQ=255)
                                                                                                                 cggggCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCaa                                                       <  1:136512/120‑1 (MQ=255)
                                                                                                                 cTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCaa                                                       >  2:136512/1‑122 (MQ=255)
                                                                                                                    ggCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTc                        <  1:468689/150‑1 (MQ=255)
                                                                                                                       aaCAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGAc                     >  2:274058/1‑150 (MQ=255)
                                                                                                                          aGCCCGTCACTGTTCGTCGACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATg                  <  2:306652/150‑1 (MQ=255)
                                                                                                                          aGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCgatga                                      >  2:474451/1‑130 (MQ=255)
                                                                                                                          aGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCgatga                                      <  1:474451/130‑1 (MQ=255)
                                                                                                                             ccGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGcc               <  2:489765/150‑1 (MQ=255)
                                                                                                                                       tCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGa     <  2:574803/150‑1 (MQ=255)
                                                                                                                                          tCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATgg  >  1:234119/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
AGGAAGGCGTGGCGCTGATGTGCCTGGCCGAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGG  >  NC_002947/5630994‑5631281

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