Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,631,117 A→T 7.5% K216M (AAG→ATG) ‡ 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,1170AT 7.5% 158.9 / 2.2 53K216M (AAG→ATG) ‡putAproline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (23/26);  new base T (2/2);  total (25/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 = 1.00e+00

CCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCATTGCC  >  NC_002947/5630977‑5631262
                                                                                                                                            |                                                                                                                                                 
ccTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGt                                                                                                                                             >  7:93034/1‑147 (MQ=255)
 ctgctgCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCAGCTTGCCGGTCAg                                                                                                                                          <  6:364304/149‑1 (MQ=255)
   gctgcGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAg                                                                                                                                                  <  2:281011/139‑1 (MQ=255)
   gctgcGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAg                                                                                                                                                  >  1:281011/1‑139 (MQ=255)
           tCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAAc                                                                                                                                >  5:82142/1‑149 (MQ=255)
             cccGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACg                                                                                                                               >  5:84788/1‑148 (MQ=255)
              ccGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGa                                                                                                                              >  2:260375/1‑148 (MQ=255)
               cGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGaaa                                                                                                                            <  8:45229/149‑1 (MQ=255)
                 aCAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAAcc                                                                                                                          >  8:236333/1‑149 (MQ=255)
                    aatgTACCCTTGATGCGCTGATCCGCTACAAGATCAGTACCGGCAACTGTCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGc                                                                                                                       <  8:93034/146‑1 (MQ=255)
                                        aTCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCAcc                                                                                                     <  3:261047/147‑1 (MQ=255)
                                         tCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCATCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCg                                                                                                   >  2:130441/1‑148 (MQ=255)
                                          ccGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGca                                                                                                 >  6:151803/1‑149 (MQ=255)
                                               aCAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCatcatc                                                                                            >  6:121829/1‑149 (MQ=255)
                                               aCAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCatcatc                                                                                            <  6:82142/149‑1 (MQ=255)
                                                  aGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGc                                                                                         <  7:236333/149‑1 (MQ=255)
                                                        gTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAgc                                                                                   <  8:190044/149‑1 (MQ=255)
                                                           ccGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTc                                                                                                                    <  2:252353/113‑1 (MQ=255)
                                                           ccGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTc                                                                                                                    >  1:252353/1‑113 (MQ=255)
                                                            cGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGACACTGTTCCTCAACGCCGCCACCTCGGCCCTCCTGCTGACCCGCAAGCTGCTCAGCACCCATAACGAAACCGGCCTCACCTCCTCCCTCACCCGCATCATCGGCAAGAgcggcg                                                                               >  1:155341/1‑149 (MQ=255)
                                                             ggCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAgcggc                                                                                >  8:181047/1‑147 (MQ=255)
                                                             ggCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGa                                                                              <  4:315047/149‑1 (MQ=255)
                                                             ggCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGa                                                                              <  5:121829/149‑1 (MQ=255)
                                                                         ccGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCg                                                                  >  4:341406/1‑149 (MQ=255)
                                                                          cGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGTGGGCTGCTGCTGACCGGCAAGCGGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGc                                                                 <  4:290794/149‑1 (MQ=255)
                                                                           gCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCgg                                                                                                                        >  1:13346/1‑93 (MQ=255)
                                                                           gCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCgg                                                                                                                        <  2:13346/93‑1 (MQ=255)
                                                                             aCCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAg                                                              <  5:57599/149‑1 (MQ=255)
                                                                             aCCTGGGCAACACCCCGTCACTGTTCGTCAACGCCGCCACCTGGTGGCTGCTGCTGACCGTCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAg                                                              <  7:79005/149‑1 (MQ=255)
                                                                                 gggCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCg                                                          >  1:198501/1‑149 (MQ=255)
                                                                                  ggCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCg                                                          <  2:198501/148‑1 (MQ=255)
                                                                                       cAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAgggatggtgg                                                     >  1:296994/1‑141 (MQ=255)
                                                                                       cAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAgg                                                             <  2:296994/140‑1 (MQ=255)
                                                                                          cccGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATgg                                                 >  3:128920/1‑149 (MQ=255)
                                                                                                    ttCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCgg                                                                                          <  5:159218/98‑1 (MQ=255)
                                                                                                    ttCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCgg                                                                                          >  6:159218/1‑98 (MQ=255)
                                                                                                    ttCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGcc                                        <  7:181047/148‑1 (MQ=255)
                                                                                                        tCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCAt                                                                                                                                   >  5:179117/1‑53 (MQ=255)
                                                                                                        tCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCAt                                                                                                                                   <  6:179117/53‑1 (MQ=255)
                                                                                                          aaCGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATgg                                  <  4:128920/148‑1 (MQ=255)
                                                                                                             gccgccACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCCCCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGt                                >  1:334553/1‑147 (MQ=255)
                                                                                                              ccgccACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCAGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAg                             <  6:84788/149‑1 (MQ=255)
                                                                                                               cgccACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAg                             >  3:8818/1‑148 (MQ=255)
                                                                                                                  cACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGt                         <  3:341406/149‑1 (MQ=255)
                                                                                                                         gggCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTctc                  <  5:151803/149‑3 (MQ=255)
                                                                                                                            tgggtgctgACCAGCTTGCCGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCa                           <  2:123368/133‑1 (MQ=255)
                                                                                                                             tgctgctgACCGGCAAGCTGGTCAGCAACCTTAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCTTCGGCATGAGCGGCGAGCCGATGATCCGCATGGGCGTCGACATGGCCATGCGACTGTTGGGTGAGCTGTTCGTCACCGGc              <  2:155341/149‑1 (MQ=255)
                                                                                                                               ctgctgACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCga            <  1:260375/149‑1 (MQ=255)
                                                                                                                                tgctgcCCAGCTTGCCGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCgag           >  5:81455/1‑149 (MQ=255)
                                                                                                                                tgctgACCAGCTTGCCGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGc              >  3:144403/1‑146 (MQ=255)
                                                                                                                                       cGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATgg                                  >  2:177313/1‑119 (MQ=255)
                                                                                                                                       cGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATgg                                  <  1:177313/119‑1 (MQ=255)
                                                                                                                                          cAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCATTGcc  >  4:6977/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
CCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCATTGCC  >  NC_002947/5630977‑5631262

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