Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,138,006 T→G 6.2% T344P (ACG→CCG)  gldA ← glycerol dehydrogenase, NAD+ dependent; 1,2‑propanediol:NAD+ oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,138,0060TG 6.2% 232.7 / 7.4 64T344P (ACG→CCG) gldAglycerol dehydrogenase, NAD+ dependent; 1,2‑propanediol:NAD+ oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base G (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

TCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATT  >  NC_000913/4137897‑4138118
                                                                                                             |                                                                                                                
tCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGttt                                                                                     >  1:105551/1‑139 (MQ=255)
tCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGttt                                                                                     <  2:105551/139‑1 (MQ=255)
         cGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCg                                                                   >  4:887348/1‑148 (MQ=255)
         cGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCg                                                                   >  1:153891/1‑148 (MQ=255)
                   gggAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATgg                                                                                        <  2:88469/117‑1 (MQ=255)
                   gggAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATgg                                                                                        >  1:88469/1‑117 (MQ=255)
                   gggAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTc                                                               >  2:784063/1‑142 (MQ=255)
                   gggAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTc                                                               <  1:784063/142‑1 (MQ=255)
                       gTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAAtt                                                     <  3:887348/148‑1 (MQ=255)
                       gTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAAtt                                                     <  2:153891/148‑1 (MQ=255)
                                aGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAAtt                                                     >  4:246871/1‑139 (MQ=255)
                                aGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAAtt                                                     <  3:246871/139‑1 (MQ=255)
                                   ttATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGcc                                         >  2:64875/1‑148 (MQ=255)
                                    tATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCg                                                   >  1:1058756/1‑137 (MQ=255)
                                    tATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCg                                                   <  2:1058756/137‑1 (MQ=255)
                                        cccACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGAc                                    >  1:42910/1‑148 (MQ=255)
                                        cccACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGAc                                    >  2:742742/1‑148 (MQ=255)
                                         ccACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCt                                                                                <  4:961843/103‑1 (MQ=255)
                                         ccACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCt                                                                                >  3:961843/1‑103 (MQ=255)
                                            ctctTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAt                                                                                                   <  4:980675/81‑1 (MQ=255)
                                            ctctTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAt                                                                                                   >  3:980675/1‑81 (MQ=255)
                                              ctTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtgt                                                                                                 >  4:319538/1‑81 (MQ=255)
                                              ctTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtgt                                                                                                 <  3:319538/81‑1 (MQ=255)
                                              ctTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAAtt                                                     >  1:978511/1‑125 (MQ=255)
                                              ctTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAAtt                                                     <  2:978511/125‑1 (MQ=255)
                                              ctTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCAtttt                                             >  1:887229/1‑133 (MQ=255)
                                              ctTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCAtttt                                             <  2:887229/133‑1 (MQ=255)
                                                    ggAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAat                        <  1:742742/148‑1 (MQ=255)
                                                      aaaCGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCggg                                      >  2:627635/1‑132 (MQ=255)
                                                      aaaCGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCggg                                      <  1:627635/132‑1 (MQ=255)
                                                               ccGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAg             >  4:385367/1‑148 (MQ=255)
                                                                 gTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCa                                                                                                    <  4:562978/59‑1 (MQ=255)
                                                                 gTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCa                                                                                                    >  3:562978/1‑59 (MQ=255)
                                                                  tACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTc                                                                                    <  1:1078726/74‑1 (MQ=255)
                                                                  tACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTc                                                                                    >  2:1078726/1‑74 (MQ=255)
                                                                   aCTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCgag         <  3:385367/148‑1 (MQ=255)
                                                                        tCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACAtctt                               <  4:272706/121‑1 (MQ=255)
                                                                        tCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACAtctt                               >  3:272706/1‑121 (MQ=255)
                                                                         cGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTAt   <  2:42910/148‑1 (MQ=255)
                                                                          ggCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTAtt  <  1:64875/148‑1 (MQ=255)
                                                                           gCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCgag         >  2:66872/1‑140 (MQ=255)
                                                                           gCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCgag         <  1:66872/140‑1 (MQ=255)
                                                                                cagcagAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGt                 <  2:469081/127‑1 (MQ=255)
                                                                                cagcagAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGt                 >  1:469081/1‑127 (MQ=255)
                                                                                 agcagAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAAt                                                      <  1:142503/89‑1 (MQ=255)
                                                                                 agcagAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAAt                                                      >  2:142503/1‑89 (MQ=255)
                                                                                   cagAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCgag         >  3:586841/1‑132 (MQ=255)
                                                                                   cagAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCgag         <  4:586841/132‑1 (MQ=255)
                                                                                            gTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATcc                    <  1:413852/112‑1 (MQ=255)
                                                                                            gTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATcc                    >  2:413852/1‑112 (MQ=255)
                                                                                                 ccTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCa                                                 >  3:495998/1‑78 (MQ=255)
                                                                                                 ccTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCa                                                 <  4:495998/78‑1 (MQ=255)
                                                                                                 ccTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTAt   >  3:1066883/1‑124 (MQ=255)
                                                                                                 ccTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTAt   <  4:1066883/124‑1 (MQ=255)
                                                                                                      tCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTg                                                                            <  4:979056/46‑1 (MQ=255)
                                                                                                      tCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATcc                    >  2:111323/1‑102 (MQ=255)
                                                                                                      tCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATcc                    <  1:111323/102‑1 (MQ=255)
                                                                                                      tCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCa                   >  4:35440/1‑103 (MQ=255)
                                                                                                      tCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCa                   <  3:35440/103‑1 (MQ=255)
                                                                                                      tCTGGCGTCGCGCCGCCAGGCATGTTGCGAATGGTTTCACCTTCTg                                                                            >  3:979056/1‑46 (MQ=255)
                                                                                                        tGGCGGCGCGACGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTa                          <  4:108208/94‑1 (MQ=255)
                                                                                                        tGGCGGCGCGACGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTa                          >  3:108208/1‑94 (MQ=255)
                                                                                                        tGGCGGCGCGACGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTAtt  <  2:1059671/118‑1 (MQ=255)
                                                                                                        tGGCGGCGCGACGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTAtt  >  1:1059671/1‑118 (MQ=255)
                                                                                                             |                                                                                                                
TCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATT  >  NC_000913/4137897‑4138118

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: